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Preface 



The purpose and motivation of these lectures can be summarized in the following 
two questions: 

• What is the ground state (and its properties) of dense matter? 

• What is the matter composition of a compact star? 

The two questions are, of course, strongly coupled to each other Depending on your 
point of view, you can either consider the first as the main question and the second 
as a consequence or application of the first, or vice versa. 

If you are interested in fundamental questions in particle physics you may take 
the former point of view: you ask the question what happens to matter if you squeeze 
it more and more. This leads to fundamental questions because at some level of suf- 
ficient squeezing you expect to reach the point where the fundamental degrees of 
freedom and their interactions become important. That is, at some point you will 
reach a form of matter where not molecules or atoms, but the constituents of an 
atom, namely neutrons, protons, and electrons, are the relevant degrees of freedom. 
This form of matter, and its variants, constitute one important topic of these lec- 
tures and is termed nuclear matter. If you squeeze further, you might reach a level 
where the constituents of neutrons and protons, namely quarks and gluons, become 
relevant degrees of freedom. This form of matter, termed quark matter or strange 
quark matter, is the second important subject we shall discuss. By studying dense 
matter, we shall thus learn a lot about the fundamental theories and interactions of 
elementary particles. When trying to understand this kind of dense matter, we would 
like to perform experiments and check whether our fundamental theories work or 
whether there are new phenomena, or maybe even new theories, that we have not 
included into our description. Unfortunately, there are currently no experiments on 
earth which can produce matter at such ultra-high densities we are talking about. 
However, this does not mean that this kind of matter does not exist in nature. On 
the contrary, we are pretty sure that we have observed objects that contain matter 
at ultra-high density, namely compact stars. We may thus consider compact stars 
as our "laboratory". Thinking about the first question has therefore led us to the 
second. 
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If you are primarily interested in phenomenology, or if you are an astrophysicist, 
you may start from the second question: you observe a compact star in nature and 
would like to understand its properties. In this case you start from observations like 
the rotation frequency, the temperature of the star etc. and ask, why does the star 
rotate so slow/so fast, why does it cool down so slow/so fast? And these questions 
will inevitably lead you to the attempt to figure out the microscopic structure of the 
star, although you have started from macroscopic observables. You need to know 
whether the star contains nuclear matter or quark matter or both, in which phase the 
respective matter is present, and which properties these phases have. It is thus very 
natural, also from the astrophysicist's point of view, to study the first question. 

In any case, we see that both questions are closely related and we don't have 
to decide which of the two points of view we take. If I have to characterize what 
awaits you in these lectures I would nevertheless say that we shall lean a bit more 
towards the fundamental aspects. In other words, we shall neglect many complica- 
tions that arise from considering a realistic compact star. A star is a finite system, it 
is inhomogeneous, it underlies the laws of general relativity etc. Although our dis- 
cussions are always motivated by the astrophysical application, we mostly discuss 
infinite, homogeneous systems and do not elaborate on general relativistic effects. 
Only in discussing the consequences of our microscopic calculations we shall, on a 
qualitative level, discuss the more realistic setting. 

So what kind of physics will we discuss and which theoretical tools do we need? 
Since our focus is on nuclear and quark matter, the dominant interaction that gov- 
erns the states of matter we are interested in is the strong interaction. The underlying 
theory for this interaction is Quantum Chromodynamics (QCD). Although this the- 
ory is uniquely determined by very simple symmetry principles, it is extremely hard 
to solve for most applications. Unfortunately (or fortunately, because this makes it 
interesting and challenging) matter at compact star densities eludes rigorous first- 
principle calculations. Therefore, we often have to retreat to simple phenomenolog- 
ical models or have to perform rigorous QCD calculations at asymptotically large 
densities and then extrapolate the results down to the density regime we are inter- 
ested in. 

In the physics of compact stars also the weak interaction plays an important role. 
We shall see that it is responsible for the chemical equilibration of the system, i.e., 
it fixes the various chemical potentials. It is also important for the understanding of 
cooling mechanisms of the star or for transport properties of nuclear and quark mat- 
ter. Furthermore, our (mostly field-theoretical) treatment always includes nonzero 
chemical potentials and sometimes nonzero temperature (for many applications the 
zero-temperature approximation is sufficient). In this sense it goes beyond the stan- 
dard vacuum field theory formalism. Basic elements of thermal quantum field theory 
at finite chemical potential are therefore explained in the appendix. 

All this may sound exciting on the one hand, because it shows that the physics 
of compact stars is extremely rich (due to the diversity of involved physics I found 
it helpful to include a glossary of important terms at the end of these lecture notes). 
But on the other hand it may also sound like a big challenge for you if you are not fa- 
miliar with advanced field theory. Nevertheless, these lecture notes are not primarily 
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intended as a review for researchers (although they might find it useful too) but as 
a pedagogical introduction for graduate students and advanced undergraduate stu- 
dents. For some of our discussions all you need as a prerequisite is some knowledge 
in thermodynamics and statistical physics, for instance in chapter |2] which deals 
almost exclusively with noninteracting systems. Some other sections, for instance 
the calculation of the neutrino emissivity in chapter|5]indeed makes use of advanced 
field-theoretical methods at finite temperature. It is not the intention of these lectures 
to develop the theoretical tools in all details before we use them. More importantly, 
all calculations are physically motivated, thus by understanding the physics behind 
the results and calculations, these lectures aim at making you familiar with the the- 
ories and technicalities via "learning by doing". So at the end of these lectures you 
will have heard about the basic phenomena and possible microscopic explanations 
of the physics of compact stars, but also will be prepared to start theoretical re- 
search in this exciting field yourself, to possibly contribute to the answers to the two 
questions we have started with. 

These lectures are based on a course given in the summer semester 2009 at the 
Vienna University of Technology. I thank all participants of this course for the lively 
discussions and the numerous questions and comments that helped improve these 
lecture notes. I am also grateful to M. Alford, P. Jaikumar, P. van Nieuwenhuizen, 
F. Preis, A. Rebhan, T. Schafer, I. Shovkovy, S. Strieker, N.-O. Walliser, Q. Wang, 
and F. Weber for helpful comments and discussions. 



Vienna, January 2010 



Andreas Schmitt 
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Chapter 1 
Introduction 



1.1 What is dense matter? 

The QCD phase diagram collects the equilibrium phases of QCD in the plane of 
quark (or baryon) chemical potential jj. and temperature T. We show a sketch of 
this phase diagram in Fig. 11.11 In this introduction, we are not concerned with the 
details of this diagram. We observe that compact stars, on the scales of this diagram, 
live in the region of small temperatures and intermediate densities. They may live 
in the region where quarks are confined, i.e., in the hadronic phase. This would 
imply that they are neutron stars. They may also live in the deconfined region which 
would make them quark stars. A compact star may also contain both deconfined 
and confined quark matter because the star actually has a density profile rather than 
a homogeneous density. In the interior, we expect the density to be larger than at 
the surface. Therefore, the third possibility is a hybrid star with a quark core and a 
nuclear mantle. 

We do currently not know the exact location of most of the phase transition lines 
in Fig. II. II Therefore, we do not know the ground state of strongly-interacting quark 
(or nuclear) matter at the relevant density. As a consequence, we can to some extent 
only speculate about the matter composition of the star The reason is, simply speak- 
ing, that QCD is notoriously hard to solve for temperatures and densities present in 
a compact star With the help of the phase diagram we can put this statement in a 
wider context: QCD is asymptotically free, which means that for large momentum 
exchange the interaction becomes weak. Hence, at sufficiently large temperatures 
and/or densities, we deal with weakly interacting quarks and gluons. In the case of 
large densities (or large chemical potentials) this can be understood from the uncer- 
tainty principle which relates small distances (the interacting particles are very close 
to each other) to large momenta. As a result of asymptotic freedom, regions in the 
phase diagram where jj. and/or T are sufficiently large can be understood from rigor- 
ous first-principle calculations. These regions, although theoretically under control, 
are far from being experimentally (even astrophysically) accessible. 
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Fig. 1.1 Conjectured phase diagram of QCD in the plane of quark chemical potential fl and tem- 
perature T. While matter at low density and high temperature is probed in heavy-ion collisions, 
cold and dense matter can only be found in neutron stars (compact stars). We may find (superfluid) 
nuclear matter and/or deconfined quark matter inside a star. Deconfined quark matter is, at high 
temperatures, termed quark-gluon plasma (QGP) and is, at low temperatures, expected to be in a 
color-superconducting state, here labelled by CFL (color-flavor locking), discussed in Sec. 14.21 and 
non-CFL (some color superconductor other than CFL). 



If we now go to lower temperatures and densities we have to cross a large un- 
known territory. Only at small temperatures and densities, when we are deep in the 
hadronic phase we have reached an area which again is under control, at least to 
some extent. Theoretically, it is more complicated than the perturbatively treatable 
asymptotic regions. After all, hadrons are quite complicated objects once we try to 
describe them in terms of their constituents. However, we can use effective descrip- 
tions which can be supported, confirmed, and improved by experiments in the lab. 
Furthermore, at least for vanishing chemical potentials, we can perform brute-force 
QCD calculations on the computer which gives us solid theoretical knowledge for 
certain quantities (at nonvanishing chemical potentials these calculations are prob- 
lematic due to the so-called sign problem). 

We thus see that compact stars (as well as the quark-gluon plasma created in 
heavy-ion collisions) reside in a region of the phase diagram which is hard to ac- 
cess. More positively speaking, this region is interesting and challenging because 
exciting and unknown physics may be discovered and new theoretical tools may 
need to be developed. Or, in other words, the cold and dense matter we talk about 
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in these lectures is interesting from the theoretical point of view because, on the 
characteristic scale of QCD, it is only moderately, not extremely, dense. 

The theoretical tools used in current research to describe cold and dense matter 
are based on the above observations: if we describe quark matter we may use per- 
turbative methods which are valid at asymptotically large densities and extrapolate 
the results down to intermediate densities. We shall do so for instance in chapter|2] 
where we treat quarks as noninteracting or in Sec. l4.3l where we calculate the color- 
superconducting gap within perturbative QCD. However, we have to be aware that 
the extrapolation of the results pushes the calculations out of their range of validity 
by many orders of magnitude. On the other hand, we may use models for nuclear 
matter which are established at low densities by experimental data. We do so for 
instance in chapter |3] This time we have to extrapolate to larger densities. Again, 
the extrapolation is in principle uncontrolled. 

These theoretical challenges emphasize the significance of astrophysical obser- 
vations: we do not simply like to confirm the results of our calculations by using 
astrophysical data, we need astrophysical input to understand the theory which we 
believe to be the underlying theory of strongly interacting matter, namely QCD. 
Therefore, the connection between astrophysical observables and microscopic cal- 
culations is one of the main subjects of these lectures. 



1.2 What is a compact star? 

Only beaten by black holes, compact stars are the second-densest objects in nature. 
They have masses of the order of the mass of the sun, M ~ 1.4M0, but radii of 
only about ten kilometers, R ^ 10 km. Thus the mass of the sun Mq = 1 .989 • 10^^ g 
is concentrated in a sphere with a radius 10^ times smaller than that of the sun, 
Rq = 6.96 • 10^ km. We thus estimate the average mass density in a compact star to 
be 

p~7•10"^gcm-^ (1.1) 

This is a few times larger than the density present in heavy nuclei, the nuclear ground 
state density 

Po = 2.5-10'4gcm-\ (1.2) 

which corresponds to a baryon number density of «o — O.lSfm^^. Mass and radius 
of the star are determined by the equation(s) of state of the matter phase(s) inside 
the star This is the subject of chapter|2] 

In the traditional picture of a compact star, the star is made out of neutron-rich 
nuclear matter. Hence the traditional name is actually neutron star. This is some- 
times still the preferred term, even if one talks about a star that contains a quark 
matter core (which then might be called "exotic neutron star"). Here we shall al- 
ways use the more general term compact star to include the possibilities of more 
exotic matter; after all, a significant part of these lectures is about this exotic matter. 
The term compact star will thus be used in these lectures for an object with charac- 
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teristic mass, radius etc. as given in this subsection. It can either be made of nuclear 
matter or variants thereof {neutron star), of a quark matter core with a surround- 
ing mantle of nuclear matter {hybrid star) or exclusively of (strange) quark matter 
{quark star or strange irflr)Q Here are a few more basic properties of compact stars: 

• Compact stars are born in a supernova, a spectacular explosion of a giant or 
supergiant star due to the gravitational collapse of its core. Supernovae are very 
complex, nonequilibrium processes which astrophysicists try to understand with 
hydrodynamic simulations. We shall not be concerned with supernovae in these 
lectures but should keep in mind that some properties of the star may be a result 
of these violent explosions. A possible example is the high velocity with which 
some of the compact stars travel through space. 

• Compact stars are not only extreme with respect to their density. Some of them 
also rotate very fast with rotation periods in the millisecond regime, such that 
their frequencies are 

v<lms"'. (1.3) 

To see that this is really fast, notice that a point on the equator has a velocity 
of 2;r/?/lms ~ 0.2c, i.e., it moves with 20% of the speed of light. The current 
record holder is the star PSR J1748-2446ad0 rotating with a period of 1.39 ms. 
Several observations are related to the rotation frequency. First of all, compact 
stars have been discovered as pulsars, by observing pulsating radio signals, for 
the first time in 1967. These periodic signals are due to the lighthouse effect, i.e., 
radio emission is aligned in a beam along the magnetic axis of the pulsar which 
spins around the rotation axis, crossing the earth's telescopes periodically. More 
interestingly for our purpose, the pure fact that the rotation of some compact stars 
can be so fast requires some explanation. From the microscopic point of view, this 
is related to transport properties such as viscosity of the matter inside the star, see 
Sec. 16.21 Also pulsar glitches, sudden jumps in the rotation frequency, must find 
an explanation in the properties of dense matter. 

• Compact stars also have huge magnetic fields, 

B-IO'^G. (1.4) 

Even larger surface magnetic fields of the order of B lO'^ G have been observed 
(the magnetic field in the core of the star possibly being even higher). Such highly 
magnetized stars are also termed magnetars. Compare these magnetic fields for 
instance to the earth's magnetic field, B ^ 0.6 G, a common hand-held magnet. 



' The term compact star is in general also used to include white dwarfs, stars which are less dense 
than neutron stars, hybrid stars, or quark stars, and sometimes even to include black holes. Since 
we are not concerned with either of these objects here, we can reserve the term compact star as 
explained in the text. 

^ The label of the star says that it is a "Pulsating Source of Radio emission" (PSR) located on the 
celestial sphere at right ascension 17 h 48 min with —24° 46' declination; the 'J' indicates the use 
of the J2000 coordinate system, the suffix 'ad' is used to distinguish the object from other pulsars 
in the same globular cluster Terzan 5. 
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B ^ lOOG, or the strongest steady magnetic fields in the laboratory, B ^ 4.5 • 
lO^G. 

• Compact stars are cold. This may sound odd, given their temperatures which, 
right after they are born in a supernova explosion, can be as high as T lO'^ K. 
This corresponds, in units where the Boltzmann constant is one, A:^ = 1, to 
T ^ lOMeV. During the evolution of the star, the temperature decreases down 
to temperatures in the keV range. The dominant cooling mechanism is neutrino 
emission which we discuss in chapter |5] There are two reasons why in our con- 
text it is appropriate to call compact stars cold, in spite of the apparently large 
temperatures. First, temperatures in the keV range are small compared to the 
scale set by QCD, for instance the deconfinement transition at vanishing quark 
chemical potential of about ~ 170MeV. This means compact stars are located 
basically on the horizontal axis in the QCD phase diagram in Fig. 11.11 Second, 
temperatures in compact stars are small compared to the quark (or baryon) chem- 
ical potential, T <C /i. This is important for our calculations since it implies that 
r = is a good approximation in many cases. 



1.3 Further reading 

Before we start with the main part, let's mention some literature. Extensive text- 
books about compact stars are Refs. ||Tl|2][3l . A shorter introduction to compact stars 
and dense matter can be found in the review article Q. Similar reviews are Refs. 
|l5]|6l, with emphasis on quark matter, and Ref. ||71, with emphasis on astrophys- 
ical observations. A more theoretical review about quark matter (more precisely, 
about color-superconducting quark matter), with a section about compact star appli- 
cations is Ref. fSi. For an introduction to thermal field theory see the textbooks ||9l 
and ifTol . on which the appendix of these lecture notes is partially based. As became 
clear above, in this course we shall deal with questions which are under debate in 
current research. Therefore, some of the material included here has so far only been 
available in research papers. The respective references will be given in the various 
chapters. I will not try to be exhaustive in the reference list but rather point out se- 
lected references which are useful for a deeper understanding of what we discuss 
in these lectures. If you are interested in more references you can find plenty in the 
quoted papers and textbooks. 
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Chapter 2 

Mass and radius of the star 



In this chapter, we will discuss the most basic properties of a compact star, its mass 
and radius. We have already given typical values for these quantities above. Below 
we shall connect them with microscopic properties of nuclear and quark matter This 
connection is made by the equation of state from which, in particular, an estimate for 
the maximum mass of the star can be obtained. Let us begin with a simple estimate 
of mass and radius from general relativity. For the stability of the star we need R> Rs 
where R is the radius of the star, and R^ = 2MG the Schwarzschild radius, with the 
mass of the star M and the gravitational constant G = 6.672 • 10^" m^kg^'s^^ = 
6.707 • lO^^^GeV^^. (We shall mostly use units common in nuclear and particle 
physics, Ti = c = ks = I, although astrophysicists often use different units.) For 
R < Rs the star becomes unstable with respect to the collapse into a black hole. Let 
us build a simple star made out of a number of nucleons A with mass m ~ 939 MeV 
and a distance fq ~ 0.5 • 10^'^ cm (that's where the nucleon interaction becomes 
repulsive). We thus cover a volume ^ r^A and thus have a radius 

R^roA'/\ (2.1) 
(for our rough estimate we are not interested in factors of n), and a mass 

M - mA . (2.2) 



Now from the limit R = 2MG we obtain 

\ ZmG 



2.6 -10^^ (2.3) 



This is the number of nucleons up to which we can fill our star before it gets unsta- 
ble. In other words, the Schwarzschild radius is proportional to the mass of the star 
and thus increases linearly in the number of nucleons A, while the radius increases 
with A'/^; therefore, for A smaller than the limit A ^ 2.6 • 10^^ the star is stable, 
while it collapses into a black hole for nucleon numbers larger than this limit. We 
can plug the limit value for A back into the radius and mass of the star to obtain 
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2 Mass and radius of the star 



R^Tkm, M-2.3M0. (2.4) 

Adding more nucleons would make the star too heavy in relation to its radius. We 
see that these values are not too far from the observed ones given in Sec. II. 21 

Besides giving an estimate for the baryon number in the star, we see from 
this simple exercise that general relativistic effects will be important because the 
Schwarzschild radius will be a significant fraction of the radius of the star We can 
also estimate the gravitational energy of the star. To this end, we need the differential 
mass of the star at a given radius (i.e., the mass of a thin spherical layer) 

dm^p{r)dV, (2.5) 

where dV = Anr^dr is the volume of the thin spherical layer at radius r. For a rough 
estimate let us (unrealistically) assume a constant density p (r) = p such that the 
mass m{r) of the star up to a radius r < is given by m{r) ~ ^r'p. Then we 
obtain 

f'^ Gm{r)dm{r) 3 GM^ 

£grav^/ — ^~-^~0.12M, (2.6) 

Jo r 5 R 

where we have used the above realistic values M ~ 1 .4M© and R ~ 10km. We thus 
see that the gravitational energy figrav is more than 10% of the mass of the star. 
This suggests that for the mass-radius relation we need an equation that incorpo- 
rates effects from general relativity. For simplicity, let us first derive the equation 
that relates mass and radius without general relativistic effects and include them 
afterwards. We are looking for an equation that describes equilibrium between the 
gravitational force, seeking to compress the star, and the opposing force coming 
from the pressure of the matter inside the star. In the case of a compact star, this 
pressure is the Fermi pressure plus the pressure coming from the strong interactions 
of the nuclear or quark matter inside. The differential pressure dP at a given radius 
r is related to the gravitational force dF via 

dF 

with 

The equation for the differential mass ( 12.5b . together with Eq. (12.71 ) (into which we 
insert Eqs. (12.5b and ( 12.8b ). yields the two coupled differential equations, 

dm 7 , , 

— ^4-Kr^e{r), (2.9a) 
dr 

dP Ge{r)m{r) 

— = 2 ■ (2.9b) 

dr r^ 



where we have expressed the mass density through the energy density e(r) ~ p(r) 
(in units where c — 1). The second equation, which is easy to understand from el- 
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ementary Newtonian physics, receives corrections from general relativity. It is be- 
yond the scope of these lectures to derive these corrections. We simply quote the 
resulting equation. 



dP 



Ge{r)m{r) 



P{r) 



e{r) 



m(r) 



2Gm(r) 



(2.10) 



This equation is called Tolman-Oppenheimer- Volkov (TOV) equation and the deriva- 
tion can be found for instance in Ref. |1J. In order to solve it, one first needs the 
energy density for a given pressure. Only then do we have a closed system of equa- 
tions. This input is given from the microscopic physics which yields an equation 
of state in the form P(e). We have thus found a first example how the microscopic 
physics can potentially be "observed" from astrophysical data, namely from mass 
and radius of the star. We shall encounter many more of these examples. The equa- 
tions of state for noninteracting nuclear and quark matter will be discussed in the 
subsequent sections. 

For a given equation of state one needs two boundary conditions for the TOV 
equation. The first is obviously m{r = 0) = 0, the second is a boundary value for the 
pressure in the center of the star, P{r = 0) = Pq- Then, the solution of the equations 
will produce a mass and pressure profile m{r), P{r) with the pressure going to zero 
at some point r ^ R, giving the radius of the star The mass of the star is then read 
off at this point, M = m{R). Doing this for varying initial pressures Pq yields a curve 
M{R) in the mass-radius plane, parametrized by Pq. This curve depends strongly on 
the chosen equation of state. 



2.1 Noninteracting nuclear matter 

We start with a very simple system where we neglect all interactions. In this case, 
all we need is basic statistical physics and thermodynamics. The thermodynamic 
potential for the grand-canonical ensemble is given by 

Q=E-iiN-TS, (2.11) 

with the energy E, chemical potential jx, particle number A^, temperature T and 
entropy S. The pressure is then 

Q 

P^-y^-e + lin + Ts, (2.12) 

where V is the volume of the system. Number density n = N/V, energy density 
e = E/V, and entropy density s = S/V are, for a fermionic system, given by 
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n = 2 



{my 



[(l-A)ln(l-/,)+/,lnM 



(2.13a) 
(2.13b) 
(2.13c) 



We shall first be interested in a system of neutrons («), protons (p), and electrons 
(e), each giving a contribution to the pressure according to Eqs. ( 12.121 ) and ( 12.131 ). 
Since they are spin-j fermions, we have included a factor 2 for the two degenerate 
spin states. The Fermi distribution function is denoted by fj^. 



fk- 



1 



and the single-particle energy is 



(2.14) 



(2.15) 



Inserting number density, energy density, and entropy density into the pressure 
dznt yields 



2T 



{my 



■In 



1- 



(2.16) 



This corresponds to the result obtained from field-theoretical methods in appendix 
IA.2I see Eq. jAJlj . There also the antiparticle contribution is included, which can 
here, due to the large positive chemical potential, safely be neglected. One can easily 
check that number density and entropy are obtained from the pressure ( 12.161 ) via the 
usual thermodynamic relations, i.e., by taking the derivatives with respect to jJ. and 
T. For the following we now take the limit T = 0. This is a good approximation 
since the temperature of a compact star is typically in the keV range and thus much 
smaller than the chemical potentials and masses of the nucleons. 

For r = the Fermi distribution is a step function, fit ~ ©{kp — k), and thus the 
k integrals will be cut off at the Fermi momentum kp, i.e.. 



1 



dkk^ 



3n^ 



(2.17a) 



1 /'^ 

e = — r / dkk^\/k^+m 

1 



Jo 



87r2 



{2k'p +m^kF ) J kp+m^-m^ In - 



kp 



(2.17b) 



Then, with pL = ^kp+rrp-, the pressure is 
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P = ^ I ' dkk^in-^/W+m2) 
n-^ Jo 



(24 - im^kp )\ kl+m^ + 3m^ In - 



kp + x kl + m^ 



(2.18) 



This can either be obtained by inserting Eqs. (I2.17ab and (I2.17bb into Eq. (12.12b or, 
equivalently, by taking the T = limit in Eq. ( 12. 16b . For the latter one makes use of 
Um7'^orin(l+e^/^) =x0(x). 

For n, p, e matter, the total pressure is 

P=\ f'' dkk^i^i-^Jk^ + mf). (2.19) 

The Fermi momenta can be thought of as variational parameters which have to be 
determined from maximizing the pressure, i.e., from the conditions 



0, i = n,p,e, (2.20) 



dP 

dkp^i 

which implies 

II, = ^kli + mj. (2.21) 

We have additional constraints on the Fermi momenta from the following two condi- 
tions. Firstly, we have to require the star to be electrically neutralfl, i.e., the densities 



^ In fact, a compact star has to be electrically neutral to a very high accuracy, as one can see from 
the following simple estimate. Suppose the star has an overall charge of Z times the elementary 
charge, Ze, and we consider the Coulomb repulsion of a test particle, say a proton, with mass m and 
charge e (e having the same sign as the hypothetical overall charge of the star Ze). The Coulomb 
force, seeking to expel the test particle, has to be smaller than the gravitational force, seeking to 
keep the test particle within the star. This gives the condition 

(Z£)£^GM^ 

with the mass M and radius R of the star. Even if we are generous with the limit on the right-hand 
side by assigning the upper limit M < Am to the mass of the star (if the star contains A nucleons, 
its total mass will be less than Am due to the gravitational binding energy), we will get a very 
restrictive limit on the overall charge. Namely, we find 

{Ze)e GAin- m- 

^ ' < — ^ Z<G^A. (2.23) 



R2 r2 e 



With the proton mass m ~ lO-'MeV, the elementary charge ~ 10 ' (remember a = e- /{4k) ; 
1/137), and the gravitational constant G ~ 7 • 10"^^^ GeV" " 



(2.24) 



i.e., the average charge per nucleon has to be extremely small in order to ensure the stability of the 
star. Since we have found such an extremely small number, it is irrelevant for the argument whether 
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of protons and electrons has to be equal, 

iif, = Hp . (2.25) 

With Eq. ( I2.17al ) this means 

kF,e = kF,p ■ (2.26) 
Secondly, we require chemical equilibrium with respect to the weak processes 

n-^ p + e + Ve, (2.27a) 
p + e ^ n + Ve. (2.27b) 

The first of these processes is the usual P -decay, the second is sometimes called 
inverse {5 -decay or electron capture. We shall assume that the neutrino chemical 
potential vanishes, /iy^ = 0. This is equivalent to assuming that neutrinos and an- 
tineutrinos, once created by the above processes, simply leave the system without 
further interaction. This assumption is justified for compact stars since the neutrino 
mean free path is of the order of the size of the star or larger (except for the very 
early stages in the life of the star). Consequently, j3 -equilibrium, i.e., equilibrium 
with respect to the processes (12.27b . translates into the following constraint for the 
chemical potentials, 

ll„=lip + lle. (2.28) 
Inserting Eq. ( 12.211 ) into this constraint yields 



kin = \jk\,, + '«p + V + ■ (2.29) 

We can eliminate the electron Fermi momentum with the help of Eq. ( 12.26b and 
solve this equation to obtain the proton Fermi momentum as a function of the neu- 
tron Fermi momentum. 

To illustrate the physical meaning of this relation, let us consider some limit cases. 
First assume a vanishing proton contribution, kf p = 0. Then the equation gives 
(which is most easily seen from Eq. ( 12. 29b ) 

kl,„^{mp + mef~ml<0. (2.31) 

This expression is negative because the neutron is slightly heavier than the elec- 
tron and the proton together, nip ~ 938.3 MeV, m„ ~ 939.6MeV, ~ O.SllMeV. 
Therefore, kf p = is impossible and there always has to be at least a small fraction 



we use a proton or an electron as a test particle. The essence of this argument is the weakness of 
gravitation compared to the electromagnetic interactions: a tiny electric charge per unit volume, 
distributed over the star, is sufficient to overcome the stability from gravity. 
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of protons. Now let's assume u = 0, which leads to 



2m„ 



-ml - 1.4 Me 



(2.32) 



This is the threshold below which there are no neutrons and the charge neutral sys- 
tem in j3 -equilibrium contains only protons and electrons of equal number density. 

In general, we may consider a given baryon density «b = «n + iip to express the 
neutron Fermi momentum as 



kF.n = {3n^nB-klp 



,1/3 



(2.33) 



Inserting this into Eq. ( |2.30t yields an equation for kf p as a function of the baryon 
density. In the ultrarelativistic Umit, i.e., neglecting all masses, Eq. ( |2.30t obviously 
yields kp^ = kf,„/2 and thus Wp = «„/8 or 



11 
9 



(2.34) 



By solving Eq. ( 12.301 1 numerically one can check that this limit is approached from 
below for large baryon densities, i.e., in a compact star containing nuclear matter 
we deal with neutron-rich matter, which justifies the term neutron star. 

As a crude approximation we may thus consider the simple case of pure neutron 
matter. We also consider the nonrelativistic limit, m„ ^ A:/r_„. In this case, the energy 
density ( I2.17bl) and the pressure ( 12.181 ) become 



(2.35a) 
(2.35b) 



(To see this, note that the In term cancels the term linear in kp in the case of e, and 
the linear and cubic terms in kp in the case of P.) Consequently, keeping the terms 
to lowest order in kf^„/m„, 



4 


k^ 






mi 




4 








ml 





P(e): 



3k' 



5/3 ./3 



15m„7r2 



(2.36) 



We have thus found a particularly simple equation of state, where the pressure is 
given by a power of the energy density. The general (numerical) discussion of the 
equation of state, including protons and electrons, is left to the reader, see problem 

The next step to obtain the mass-radius relation of the star is to insert the equation 
of state into the TOV equation. The simplest case is a power-law behavior as in Eq. 
(12.361 1. The general form of such a so-called "polytropic" equation of state is 
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P{e)=Ke^. (2.37) 

Using the Newtonian form of the mass-radius equations, Eqs. ( 12. 9b , this yields 

dm 47Z , I /v/ N 
^ = ^^^^'/^«, 

dP _ G P^ly{r)m{r) 
'd? ^ 

Even in this simplest example, we need to solve the equations numerically, see prob- 
lem [Z2] The results of this problem show that the maximum mass reached within 
this model is about M < 0.7M© which is well below observed neutron star masses. 
(See also Refs. 121 El HI for a pedagogical introduction into the equation of state and 
mass-radius relation from solving the TOV equation.) This small maximum mass is 
a consequence of the assumption of noninteracting nucleons. Taking into account 
interactions will increase the maximum mass significantly. 



(2.38a) 
(2.38b) 



2.2 Noninteracting quark matter 

Whenever we talk about quark matter in these lectures we ignore the charm (c), 
bottom (b), and top (?) quarks. The quark chemical potential inside the star is at 
most of the order of 500 MeV and thus much too small to create a population of 
these states. Therefore, we only consider at most three quark flavors, namely up (m), 
down (d), and strange (s). We shall mostly neglect the masses of the u and d quarks; 
their current masses are m„ ~ m^i ~ 5 MeV <C /X ~ (300 — 500) MeV. The mass 
of the strange quark, however, is not negligible. The current strange quark mass 
is ms ~ 90 MeV, and the density-dependent constituent mass can be significantly 
larger, making it non-negligible compared to the quark chemical potential. 

If we consider free quarks, the energy density e, the number density «, and the 
pressure P, are obtained in the same way as demonstrated for nucleons in the previ- 
ous subsection. We only have to remember that there are three colors for each quark 
flavor, Nc = 3, i.e., the degeneracy factor for a single quark flavor is 2Nc = 6, where 
the factor 2 counts the spin degrees of freedom. Consequently, for each quark flavor 
f — u,d,s we have at zero temperature (cf. Eqs. (I2.17l i. (12.18b ). 

nf = ^, (2.39a) 



£/ ^ ^ ' dkk^^k^+mj, (2.39b) 



Jo 



'F.f 



-2 ' dkk^ijXf- Jk^ +mj) . (2.39c) 
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Again, we need to impose equilibrium conditions with respect to the weak interac- 
tions. In the case of three-flavor quark matter, the relevant processes are the leptonic 
processes (including a neutrino or an antineutrino) 

d^u + e + Ve, i^-M + e + Vf, (2.40a) 
u + e^d+Ve, it + e^s + Ve^ (2.40b) 

and the non-leptonic process 

s + u^d + u. (2.41) 

These processes yield the following conditions for the quark and electron chemical 
potentials, 

= lie + fiu , fis^lle+fiu- (2.42) 

(This automatically implies jU^/ ~ jUj .) The charge neutrality condition can be written 
in a general way as 

L q/nf-n^^Q, (2.43) 

f=ii,d,s 

with the electric quark charges 

2 1 

9h = ^ , qd = = - 2 ' (2-44) 

and the electron density rig. 



2.2.1 Strange quark matter hypothesis 

Before computing the equation of state, we discuss the strange quark matter hypoth- 
esis within the so-called bag model. The bag model is a very crude phenomenolog- 
ical way to incorporate confinement into the description of quark matter. The effect 
of confinement is needed in particular if we compare quark matter with nuclear 
matter (which is ultimately what we want to do in this section). Put another way, 
although we speak of noninteracting quarks, we need to account for a specific - in 
general very complicated - aspect of the interaction, namely confinement. 

To understand how the bag constant accounts for confinement, we compare the 
pressure of a noninteracting gas of massless pions with the pressure of a noninteract- 
ing gas of quarks and gluons at finite temperature and zero chemical potential. The 
pressure of a single bosonic degree of freedom at = and at large temperatures 
compared to the mass of the boson is 

r c/^k / ,.,r\ n'^T^ 

This is derived in appendix lA. 1 1 within thermal field theory, see Eq. (IA.37I I. Anal- 
ogously, a single fermionic degree of freedom gives (see Eq. (IA.72b of appendix 
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.,_.r/^i„(i+.-/^).Z5!l!. (..46) 

Therefore, since there are three types of pions, their pressure is 

Pn = 3 . (2.47) 

90 

This is a simple approximation for the pressure of the confined phase. In the decon- 
fined phase, the degrees of freedom are gluons (8 x 2) and quarks {4-NcNf = 24). 
Thus with 2 X 8 + 7/8 X 24 = 37 we have 



90 



^'?.s = 37^^-B, (2.48) 



where the bag constant B has been subtracted for the following reason. If B were 
zero, the deconfined phase would have the larger pressure and thus would be pre- 
ferred for all temperatures. We know however, that at sufficiently small tempera- 
tures, the confined phase (that's the world we live in) must be preferred. This is 
phenomenologically accounted for by the bag constant B which acts like an energy 
penalty for the deconfined phase. Without this penalty, at least in this very sim- 
ply model description, the deconfined phase would be "too favorable" compared 
to what we observe. As a consequence, by including the bag constant there is cer- 
tain critical temperature Tc below which the confined phase is preferred, Pt^ > Pq g, 
and above which the deconfined phase is preferred, Pc|,^ > Ptz. This is indeed what 
one expects from QCD, where the deconfinement transition temperature is expected 
to be Tc ~ 170 MeV. (As can be seen in the QCD phase diagram in Fig. 11.11 this 
deconfinement transition is rather a crossover than a phase transition in the strict 
sense.) 

In the context of compact stars we are not interested in such large temperatures. 
In this case, the chemical potential is large and the temperature practically zero. 
Nevertheless we compare nuclear (confined) with quark (deconfined) matter and 
thus have to include the bag constant in the pressure and the free energy of quark 
matter, 

P + B = Y,Pf, (2.49a) 

/ 

e = ^e/+fi. (2.49b) 

/ 

This phenomenological model of confinement is called the bag model ||5]|6l because 
the quarks are imagined to be confined in a bag. One can view the microscopic 
pressure Y.fPf of the quarks to be counterbalanced by the pressure of the bag B and 
an external pressure P. 

Equipped with the bag model, we can now explain the strange quark matter hy- 
pothesis. For simplicity we consider massless quarks. A nonzero strange quark mass 
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will slightly change the results but is not important for the qualitative argument. We 
will also ignore electrons. They are not present in three-flavor massless quark mat- 
ter at zero temperature. They are however required in two-flavor quark matter to 
achieve electric neutrality. But also in this case their population is small enough to 
render their effect unimportant for the following argument. 
With m/ = we simply have 

which in particular implies 

P/ = f. (2.51) 

For the strange quark matter hypothesis we consider the energy E per nucleon num- 
ber A, 

E £ 

-T = -, (2.52) 

A llB 

where is the baryon number density, given in terms of the quark number densities 
as ^ 

nB = -Y^nf, (2.53) 

^ f 

because a baryon contains A^^ = 3 quarks. At zero pressure, P — 0, Eqs. ( |2.49l l and 
( 12.51b imply e = 4B and thus 

E 4B 

T = — ■ (2-54) 

A llB 

We now apply this formula first to three-flavor quark matter ("strange quark mat- 
ter"), then to two-flavor quark matter of only up and down quarks. For strange quark 
matter, the neutrality constraint ( 12.431 ) becomes 

2«j, — rij — ris = . (2.55) 

Together with the conditions from chemical equilibrium ( 12.421 ) this implies 

jU„ = pLci = Us = At • (2.56) 

We see that strange quark matter is particularly symmetric. The reason is that the 
electric charges of an up, down, and strange quark happen to add up to zero. Now 
with ng = p? /k^ and 

^-1^/ = ^ (2.57) 



(still everything at P = 0) we have 



= (47r2)^/433/4B'/4^ 5 714^1/4^ MgVfll/j. (2.58) 

Nf=3 
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We have expressed ^^/^ in units of 145 MeV, b\{1 = ^^/^/(USMeV). 

For two-flavor quark matter (neglecting the contribution of electrons), the charge 
neutrality condition is 

nu = 2n„ . (2.59) 



Hence, 



Then, with «b = /iu/Tf^ and 



fid 



2^/3 



B = LPf 

f 



(2.60) 



(2.61) 



we find 



(47r2)i/4 (1 +24/3)3/4fii/4 ^ 6.441 b'/^ ~ 934MeVB;/^ . (2.62) 



By comparing this to Eq. ( 12.58b we see that two-flavor quark matter has a larger en- 
ergy per baryon number than three-flavor quark matter. This is a direct consequence 
of the Pauli principle: adding one particle species (and keeping the total number of 
particles fixed) means opening a set of new available low-energy states that can be 
filled, thus lowering the total energy of the system. 

We can now compare the results (12.581 1 and ( I2.62l i with the energy per nucleon in 
nuclear matter. For pure neutron matter, it is simply given by the neutron mass. 



= m„ =939.6MeV. 



(2.63) 



For iron, ^^Fe, it is 



"Fe 



56mA, -56 -8.8 MeV 
56 



930 MeV, 



(2.64) 



with the nucleon mass mn = 938.9 MeV and the binding energy per nucleon in iron 
of 8.8 MeV. Since we observe iron rather than deconfined quark matter, we know 
that 

E 
< — 

Nf=2 



b1/4 > 144.4 MeV. 



(2.65) 



We have thus found a lower limit for the bag constant from the stability of iron with 
respect to two-flavor quark matter. Now what if the bag constant were only slightly 
larger than this lower limit? What if it were small enough for three-flavor quark 
matter to have lower energy than iron? The condition for this would be 



E 



B''^ < 162.8 MeV. 



56Fe 



(2.66) 
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This would imply that strange quark matter is absolutely stable (stable at P = 0), 
while nuclear matter is metastable. This possibility, which would be realized by a 
bag constant in the window 145 Me V < B'/"* < 162MeV, is called strange quark 
matter hypothesis, suggested by Bodmer Q and Witten ||8l, see also ||9l- 

Note that the existence of ordinary nuclei does not rule out the strange quark 
matter hypothesis. The conversion of an ordinary nucleus into strange quark matter 
requires the simultaneous conversion of many (roughly speaking A many) u and 
d quarks into s quarks. Since this has to happen via the weak interaction, it is 
practically impossible. In other words, there is a huge energy barrier between the 
metastable (if the hypothesis is true) state of nuclear matter and absolutely stable 
strange quark matter. This means that strange quark matter has to be created in an- 
other way ("going around" the barrier), by directly forming a quark-gluon plasma. 
This can for instance happen in a heavy-ion collision. Or, more importantly in our 
context, it may happen in the universe, giving rise to stars made entirely out of quark 
matter, so-called strange stars. 

Small "nuggets" of strange quark matter are called strangelets (a strange star 
would then in some sense simply be a huge strangelet). If a strangelet is injected 
into an ordinary compact star (a neutron star), it would, assuming the strange quark 
matter hypothesis to be true, be able to "eat up" the nuclear matter, converting the 
neutron star into a strange star. Note the difference between this transition and the 
above described impossible transition from ordinary nuclear matter to strange quark 
matter: once there is a sufficiently large absolutely stable strangelet, successive con- 
version of up and down quarks into strange quarks increase the size of the strangelet; 
the energy barrier originating from the simultaneous creation of a large number of 
strange quarks now cannot cause the system to relax back into its original nuclear 
(metastable) state. This argument has important consequences. If there exist enough 
sizable strangelets in the universe to hit neutron stars, every neutron star would be 
converted into a strange star. In other words, the observation of a single ordinary 
neutron star would rule out the strange quark matter hypothesis. Therefore, it is 
important to understand whether there are enough strangelets around. It has been 
discussed recently in the literature that there may not be enough strangelets ifTOl . in 
contrast to what was assumed before. 



2.2.2 Equation of state 

Next we derive the equation of state for strange quark matter. We include the effect 
of the strange quark mass to lowest order and also include electrons. It is convenient 
to express the quark chemical potentials in terms of an average quark chemical 
potential jU = (/i„ + /Xj + /Xj)/3 and the electi-on chemical potential He, 
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2 

jXu^ 11- -lie , (2.67a) 

IJ.,1 = H + ^lie, (2.67b) 

jj., = Ix + ^He. (2.67c) 

Written in this form, the conditions from )3 -equilibrium ( I2.42l i are automatically 
fulfilled. Taking into account the strange quark mass, the Fermi momenta for the 
approximately massless up and down quark and the massive strange quark are given 
by 

kfM = III, , (2.68a) 
kF,d = lid , (2.68b) 



kF,, = ^llj-ml (2.68c) 
The energy density and the pressure are 

L ^' ^rr + T^ + A (2.69a) 

y P^^p^ + ^ + I^ dkk^ ( PL, - Jk^+mi] + j(2.69b) 

where we have neglected the electron mass. The neutrality condition can now be 
written as 

d ^ 2 11 

(Note that fie is defined as the chemical potential for negative electric charge.) Solv- 
ing this equation to lowest order in the strange quark mass yields 

2 

ill 

He^^. (2.71) 

Consequently, the quark Fermi momenta become 

2 

kF.u - M - ? , (2.72a) 

Ofl 

+ (2.72b) 

kF,s - M - • (2.72c) 
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Fig. 2.1 Illustration of the Fermi momenta for neutral, unpaired quark matter in j8 -equilibrium 
with quark chemical potential /i . The splitting of the Fermi momenta is due to the strange quark 
mass ojj which is assumed to be small compared to jj,. 



We see that the Fermi momenta are split by an equal distance of m^/{4-jj.), and 
kF,s < kp^ii < kf j, see Fig. 12.11 The splitting and the order of the Fermi momenta 
can be understood from the following physical picture: start from the symmetric 
situation nis = /ie = 0. In this case, all quark flavors fill their Fermi spheres to a 
common Fermi momentum given by ji, and the system is neutral. Now switch on 
the strange quark mass. This lowers the Fermi momentum of the strange quark ac- 
cording to Eq. ( I2.68cb . Consequently, there are fewer strange quarks in the system 
and thus there is a lack of negative charge. To counterbalance this missing negative 
charge, the system responds by switching on a chemical potential jie- Because of 
/3 -equilibrium, the Fermi momenta of all quark flavors are rigidly coupled to this 
change. Electric neutrality is regained by lowering the up quark Fermi momentum 
and raising the down and strange quark Fermi momenta (which makes for the catchy 
phrase "the Fermi momentum of the down goes up"). Since the strange quark Fermi 
momentum was already lowered by the finite mass, it is clear that the resulting or- 
der is kp^s < kF,u < kp^- The electron contribution to the negative charge density is 
negligibly low, °^ ji^ °^ m^/ii^, while the contribution of the quarks due to the 
strange quark mass is proportional to Hin^. The splitting of the Fermi momenta due 
to the effects of the strange quark mass, j3 -equilibrium, and electric neutrality is 
very important in the context of color superconductivity. Since color superconduc- 
tivity is usually based on Cooper pairing of quarks of different flavor, a mismatch 
in Fermi surfaces tends to disfavor this pairing. We shall discuss superconductivity 
in quark and nuclear matter in chapter |4] and give a brief qualitative discussion of 
the consequences of Fermi surface splitting for color superconductivity at the end 
of that chapter 
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Here we continue with unpaired quark matter and insert the result for ji^ ( |2.71| i 
back into the energy density and the pressure. Again keeping only terms to lowest 
order in the strange quark mass yields 



(2.73a) 



LP'^T^.-^^- (2.73b) 



Consequently, 



Is,--31:p, + ^. (2.74) 
With Eq. ( 12. 49 al l the pressure, including the bag constant, becomes 

An^ An^ 

and, expressing P in terms of the energy density, we obtain with the help of Eq. 

-^-^^ 

This is the equation of state of noninteracting, unpaired strange quark matter within 
the bag model with strange quark mass corrections to lowest order. 



2.3 Mass-radius relation including interactions 

Let us briefly discuss the results for the mass-radius relation of a compact star for 
given equations of state for nuclear and quark matter Since the underlying calcula- 
tions in general are complicated and have to be done on a computer, we only quote 
some results to illustrate the physical conclusions. So far we have only discussed 
the simplest cases of noninteracting matter. Interactions have a significant effect on 
both the equation of state and the mass-radius relation. We now discuss these effects 
briefly, only in the subsequent chapters shall we study the nature and details of these 
interactions (and discuss their relevance to other observables than the mass and the 
radius of the star). 

The maximum mass of a star for noninteracting nuclear matter is ^ O.IMq (see 
for instance Ref. IH or solve problem l2!2b : including interactions increases the mass 
to values well above IM^j. The significance of the equation of state and interactions 
for the maximum mass is easy to understand; if the pressure P{e) for a given en- 
ergy density e is large, the system is able to sustain a large gravitational force that 
seeks to compress it. Comparing two equations of state over a given energy density 
range, the one with the larger pressure (for all energy densities in the given range) is 
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thus termed stiff, the one with the smaller pressure is termed soft. Soft equations of 
state can sustain less gravitational force and thus lead to stars with lower maximum 
masses. In the case of noninteracting nuclear matter, it is only the Fermi pressure 
from the Pauli exclusion principle that prevents the star from the collapse. Inter- 
actions increase this pressure because the dominant effect in the case of nuclear 
matter at the relevant densities is the short-range repulsion between the nucleons. 
Therefore, the maximum mass is significantly larger in this case. 

In Figs. I2.2l and l2.3l several models for the nuclear equation of state are applied 
to obtain maximum masses up to 2.4Mm. For the case of quark matter, we can 
understand some of the corrections through interactions in the following simple 
way. A generalization of the pressure (12.75b is 



This equation contains two corrections compared to Eq. (|2.75l l. One is included in 
the coefficient c and originates from the (leading order) correction of the Fermi 
momentum due to the QCD coupling as. 



resulting in a correction of the ji'^ term in the pressure with c = laj n. (This mod- 
ification of the Fermi momentum will also become important in the context of neu- 
trino emissivity in chapter[5]) Higher order calculations suggest c > 0.3 at densities 
relevant for compact stars. However, the exact value of c is unknown because per- 
turbative calculations are not valid in the relevant density regime, cf. discussion in 
Sec. II. II Therefore, c can only be treated as a parameter with values around 0.3, as 
done for example in Fig. 12.21 To get an idea about perturbative calculations beyond 
leading order in a^, you may consult the recent Ref. ifTTl . 

The second correction in Eq. ( 12.771 ) is the quantity A . This is the energy gap aris- 
ing from color superconductivity whose microscopic origin we discuss in chapter]?] 
It gives a correction to the jx^ term in the pressure. One might think that this correc- 
tion is negligible compared to the jx^ term and the bag constant. However, it turns 
out that for reasonable values of the bag constant these two terms largely cancel 
each other and the /i^ term becomes important. However, the effect of superconduc- 
tivity is still hard to determine. Firstly, it would require a precise knowledge of the 
strange quark mass. Secondly, it turns out that the maximum mass of a hybrid star 
is not very sensitive to the value of — 44^ |[T2l . 

As a result of this discussion and the results in Figs. l2.2l and l2.3l two points are 
important for the further contents of these lectures. Firstly, we should now be moti- 
vated to learn more about the nature and the consequences of interactions in nuclear 
and quark matter. Secondly, we have learned that, given our ignorance of the pre- 
cise quantitative effects of the strong interaction and the uncertainty in astrophysical 
observations, the mass and the radius of the star are not sufficient to distinguish be- 
tween a neutron star, a hybrid star, and possibly a quark star. Therefore, we also have 




(2.77) 




(2.78) 
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APR + Phenomenological QM EoS 
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Fig. 2.2 Mass-radius plot from Ref. 1121 which shows the dependence of the mass-radius curve 
on the (uncertain) parameters of the quark matter equation of state in a hybrid star. We see that 
reasonable choices of the parameters lead to similar curves as for nuclear matter (here with the 
APR equation of state). In this plot, the transition density p^. (in units of the nuclear ground state 
density no) between quark matter and nuclear matter has been used as a parameter, rather than 
the bag constant. From our discussion it is clear that one can be translated into the other. The 
coefficient c describes QCD corrections to the quark Fermi momentum and thus to the term in 
the pressure, see Eq. ilJll . 



' ' Neuiron star 
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Fig. 2.3 Mass-radius plot from Ref. 1131 . A comparison of a neutron star, different hybrid stars, 
and a quark star is shown, using several nuclear equations of state (DBHF, APR, HHJ) and several 
quark phases (CFL, 2SC). For more details and explanations of the various abbreviations, see Ref. 

fS). 



to take into account other observables which are linked to the microscopic physics. 
While the equation of state is a bulk property, i.e., it is determined by the whole 
Fermi sea, there are other phenomena which are only sensitive to the low-energy 
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excitations at the Fermi surface. One class of such phenomena is given by transport 
properties. They can possibly be related to observables which are more restrictive 
than mass and radius for the question of the matter composition of the star. We 
shall discuss such observables in chapter |5] where we relate the cooling of the star 
to neutrino emissivity, and in chapter |6] where we qualitatively discuss other such 
observables. 



Problems 



2.1. Equation of state for noninteracting nuclear matter 

Find the full equation of state for noninteracting n, p, e matter at T = numerically 
by plotting P versus e. You should see the onset of neutrons and identify a region 
where the equation of state is well approximated by the power-law behavior of pure 
neutron matter in the nonrelativistic limit, Eq. ( 12.36b . 

2.2. Mass-radius relation 

(a) Solve Eqs. ( 12.38b numerically (for nonrelativistic pure neutron matter, i.e., y = 
5/3) and plot m{r), P{r) for a given value of the pressure Pq = P{r = 0). 

(b) Use Pq as a parameter to find the mass-radius relation M{R). To this end, you 
need to do (a) for several values of Pq and find for each Pq the radius R at which 
P{R) ~ and the corresponding mass M{R). 

(c) You may incorporate general relativistic effects from the TOV equation (12.10b 
and/or the full (numerical) equation of state for noninteracting nuclear matter from 
Problem|2l] 
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Chapter 3 

Basic models and properties of dense nuclear 
matter 



There are numerous models to describe cold and dense interacting nuclear matter. 
Some of them have been used to obtain the curves in Figs. l2.2l and l2.3l From these 
curves we see that the models may differ significantly in their predictions of the 
properties of neutron stars and hybrid stars. The reason is that they all are extrapo- 
lated into a regime where there is little theoretical control. In other words, for densi- 
ties below the nuclear ground state density there are experimental data for instance 
from atomic nuclei or neutron scattering which serve to fit the parameters of the nu- 
clear models unambiguously. However, it is very challenging to construct a model 
which reliably predicts the properties of nuclear matter for larger densities. Put an- 
other way, currently the only "experiments" in this density regime are astrophysical 
observations which themselves are naturally less controlled than experiments in the 
laboratory. Therefore, the state of the art in describing interacting nuclear matter at 
high densities is a competition between several models which all are prone to uncer- 
tainties. In these lectures we do not attempt to give an overview over these models. 
We rather focus on two basic models and discuss them in detail. The first is the 
Walecka model and its extensions. The second is chiral perturbation theory, which 
is an effective model based on chiral symmetry of QCD and spontaneous breaking 
thereof in nuclear matter. We shall use it to discuss kaon condensation in nuclear 
matter. 

To put the following in the perspective of understanding QCD, we should keep 
in mind that nucleons are ultimately built of quarks and gluons which are the fun- 
damental degrees of freedom of the strong interactions. It is a highly nontrivial task 
to describe even the mass of a nucleon from quarks and gluons, let alone nuclear 
interactions. An important tool for such a case is an effective theory which has non- 
fundamental degrees of freedom, baryons and mesons instead of quarks and gluons. 
An effective theory can in principle be obtained by taking the low-energy limit of 
the underlying fundamental theory, in this case QCD. However, this procedure may 
turn out to be very difficult. Therefore, one tries to "guess" an effective theory, for 
instance guided by symmetry principles. One obtains a theory with some unknown 
parameters which have to be fit, for instance to experimental results. Once the pa- 
rameters are fitted, one may extrapolate the theory beyond the regime where the fit 
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has been done. In our case, this will be the high-density region for which we have no 
experiments in the laboratory. There is of course no guarantee that this extrapolation 
works. Models for interacting nuclear matter at high densities have to be understood 
in this spirit. Of course, an upper density limit for the validity is the deconfinement 
phase transition to a phase where quarks and gluons are the relevant degrees of free- 
dom. This limit density is not precisely known but may well be reached in compact 
stars. 



3.1 The Walecka model 

The Walecka model contains nucleons which interact via the exchange of the scalar 
a meson and the vector © meson. The Lagrangian is 

= J^N + -^a ,m + -^I , (3.1) 

Here, the free nucleon Lagrangian is 

^N^w{i'fdi,-mN + ^lf)w. (3.2) 

where = ™d W = [ ^" ] with the neutron and proton spinors y/,, and Xj/p- 

For a basic discussion of the field-theoretical treatment of noninteracting fermions, 
in particular the roles of finite temperature and chemical potential, see appendix 
IA.2I The free mesonic Lagrangian is 

^a.a = ^ {d^adf^a-mla^-) -^(O^yO^" + ^ml(0^(0^ , (3.3) 

where (O^y = d^COy ~ dyCO^i, and the interaction Lagrangian with Yukawa interac- 
tions between the nucleons and the mesons is 

Ji'i ^ gaWa\if + g(oWr^co^W- (3.4) 

We shall consider isospin-symmetric matter, i.e., the masses and chemical poten- 
tials of protons and neutrons are assumed to be identical. In general, /i is a matrix 
jj. = diag(/i„,/ip) = diag(/iB + flijfiB — M/) with the baryon and isospin chemical 
potentials /ig and jj.]. Thus, in other words, we assume the isospin chemical poten- 
tial to vanish. We can then simply denote jj. = jxe = fi„ = fip. Also the interactions 
between the nucleons are assumed to be symmetric, i.e., the ««, pp, and np in- 
teractions are identical. An isospin asymmetry in the interactions can be included 
by adding p-meson exchange. We will briefly discuss this in Sec. 13.21 Also kaon 
condensation induces an asymmetry, discussed in Sec. 13.31 

The parameters of the model are the masses and the coupling constants. The 
masses are 
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mA, = 939MeV, mffl = 783MeV, = (500 - 600) MeV . (3.5) 

The (7 meson is in fact a broad resonance and thus we can only approximately assign 
a mass to this meson. Below we shall use nia = 550MeV. The additional parameters 
are the coupling constants ga, gm- We shall discuss below how they are fixed. 

In order to compute the equation of state, we need to consider the partition func- 
tion 

Z = / ^VA^v^^a^O) exp / if, (3.6) 



X 



where we abbreviated 



= I dz (3.7) 

X Jo J 

with the inverse temperature j3 = 1 /r. We shall allow for vacuum expectation values 
of the mesons. To this end, we write the meson fields as a sum of the condensate 
and fluctuations, 

a^a + a, (3.8a) 
co^ ^ OoSo^ + co^ , (3.8b) 

as explained in appendix lA.ll for a general bosonic field. Now the simplest approxi- 
mation is to neglect the fluctuations. This corresponds to the mean-field approxima- 
tion. In this case the interaction between the nucleons and the mesons is simplified 
to a mesonic background, or mesonic mean field, which is seen by the nucleons. 
We can then simply drop all derivative terms of the mesons. As a consequence, 
the meson mean fields merely act as corrections to the nucleon mass and chemical 
potential, and we obtain the Lagrangian 

if = V/ {iy'^d^ - m*f^ + n*Yo) xj/ - ^"4^^ + ^'"i^^ > (3-9) 



with 

mil = niN-gaO , (3.10a) 
PL* = PL- g 0,6)0. (3.10b) 

It is important to keep in mind that the actual chemical potential, associated with nu- 
cleon number, is /i, not pL* . This becomes important for the correct thermodynamic 
relations, see footnote before Eqs. ( 13.26b . The new effective "chemical potential" 
jX* nevertheless has physical meaning since it determines the Fermi energy as we 
shall see below. 

The partition function now becomes 

Z = eT(-i'"'^^'+i"4'36) [ ^xf,^xf,exp [ \tf{i'fd^-m% + ix*Yo)¥- (3.11) 



X 



The evaluation of the free fermionic part (with modified mass and chemical poten- 
tial) is now straightforward and is done in detail in appendix IA.2I Here we only 
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repeat the most important steps. One first introduces the Fourier transforms 

K K 



12) 



Our conventions are K = (— ia)„,k), X = (— /T,x), and K ■ X = k^XQ — k • x = 
— (a)„T + k • x), with the fermionic Matsubara frequencies o,, = (2« + l)nT. Thus, 
after performing the X integral in the exponent one obtains 



with the inverse nucleon propagator 

G-\K) = -tK^-yQpi* 



'AT- 



(3.13) 



(3.14) 



Now using the standard formula for the functional integral over Grassmann variables 
one obtains 



T 



(3.15) 



where the determinant is taken over momentum space, Dkac space, and the (here 
trivial) neutron-proton space. Consequently, 



InZ 



+4y 



V 



^ .2 =,2 . ^ 2 ,-,2 



2 

T 



+ ln (l +in (l 



, (3.16) 



where we have performed the Matsubara sum and taken the thermodynamic limit, 
and where we have defined the single-nucleon energy 



The pressure then becomes 



r 1 1 

P = - InZ = - i^mla^ + -rn^^o" + , 
with the nucleon pressure (after subtracting the vacuum part) 



AT 



{2nf 



(3.17) 



(3.18) 



(3.19) 



We have thus derived the fermionic pressure already used in chapter |2] see Eq. 
(I2.I6I 1, from thermal field theory. The factor 4 counts the two spin degrees of free- 
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dom and the two baryon degrees of freedom (proton and neutron). We also have 
obtained the contribution of antiparticles, for which jj.* — jU*. 

The meson condensates have to be determined by maximizing the pressure. We 
obtain 

dp , dPN 

0^^^-«4^-..^, (3.20a) 

In terms of the baryon and scalar densities 

dPN SPn .V- f d^^ 1 



".H*V'>^-|^^4E/^=i,,-±_-. ,3.21b, 



we can write the equations for the condensates ( |3.20| i as 

ga 



yj - ^ris, (3.22a) 
mi 

OQ^^riB. (3.22b) 



It is useful to rewrite the first of these equations as an equation for the corrected 
mass m*fj rather than for the condensate a. 



mfj=mN ^Hj, (3.23) 

ml 

where we have used Eq. (I3.10al ). We now take the zero-temperature limit, T <C 
mN,l-L, which is justified since the temperatures of interest are at most of the or- 
der of 10 MeV, while the baryon chemical potentials are above 1 GeV. The Fermi 
distribution function then becomes a step function. In particular, all antiparticle con- 
tributions vanish. We obtairQ 



' One has to be careful with the thermodynamic relations in deriving the energy density l |3.26bt : 
remember that the actual chemical potential associated with baryon number ng is /x, not /i*. This 
means that the pressure at zero temperature can be written as P = — e + flnB. The last term of the 
pressure (term in square brackets on the right-hand side of Eq. )3.26at ) comes from a term of the 
stmcture — eo + fl*nB, cf. for instance Eq. l |2.18| l. With ^* = jl — g a, Wo and the expression for Wq 
in Eq. j3.22bt we can write this as 



P = -eo + ^^*nB+ ^ J 



1 gl 2 1 gl 



2 mi, 2 mi 



^ gi 2 , ^ gl 2 
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1^ 2_ 2 

2 nim 2. mt 



+ 



1 



2 mi 



-kl-{m*f,fkF]E*p + {m*f^f\n 



kp+Et 



„2 "B+2,n2"-^ 



+ 



1 



(2kl + {m*i^fkF)E*p-lm*j^Y\n 



4 kf+Ep 



where we have defined the Fermi energy 



and where the zero-temperature densities are 
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kfEp - (mj^)^ln 



kp+Ei 



(3.26a) 

(3.26b) 
(3.27) 

(3.28a) 
(3.28b) 



Pressure and energy density in Eqs. (13.261 1 define the equation of state which has to 
be determined numerically. We may discuss the limits of small [kp — > 0) and large 
(kp — > °°) density analytically. For small density we find 



24 
37r2 



(3.29) 



neglecting terms of the order of 4 / ("^]v)^ ^'^'^ higher Therefore, from Eq. 

( 13.231 1 we conclude 

OT^~mAr, (3.30) 

where we have suppressed terms of the order of kp/m^ ^ mn- The pressure and 
the energy density are, within this small-density approximation, dominated by the 
nucleonic pressure Pn, 
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15a;2 



niN 



2mNkp 
37r2 



(3.31) 



Comparing with Eqs. ( 12.351 ) we see that we have exactly reproduced the noninteract- 
ing limit. This is no surprise because the only effect of the interactions in the present 



from which we can read off the energy density 

e = eo 

which yields Eq. j3.26bt . 



1 



1 Sl .2 



" 11- 

2 mi, 2 inf, 



(3.25) 
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Fig. 3.1 Density-dependent effective nucleon mass at T = in the Walecka model in units of 
the zero-density mass mjv and as a function of the baryon density ng. Solid line: full numerical 
result. Dashed line: high-density approximation from Eq. <3.33t . 



approach is the modification of /i and niN. In the small-density limit these effects 
are negligible and we are back to the noninteracting result, where the equation of 
state has the form P °^ e^^^. 

For large fcf , on the other hand, we have 



We see that the effective nucleon mass goes to zero for large densities. For general 
values of the Fermi momentum, the effective mass has to be computed numerically 
from Eqs. (3l3\) and dllSbl ), see Fig.im 

At large densities, the nucleonic pressure fV as well as the pressure from the 
scalar meson (which is proportional to n^/ni^) behave like kp. Therefore, the total 
pressure is dominated by the vector meson contribution which is proportional to 
rig/ml and thus behaves like kp. 




(3.32) 



and thus 




(3.33) 




(3.34) 



Consequently, the speed of sound approaches the speed of light at large densities. 
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Fig. 3.2 Binding energy per nucleon at zero temperature in the Walecka model as a function of 
baryon density, obtained from computing the energy density with the density-modified nucleon 
mass. The binding energy has a minimum at which the pressure P is zero, i.e., at this point nuclear 
matter is self-bound, and the corresponding density is called saturation density. The two parameters 
of the model, namely the coupling constants and g^, are fixed such that the binding energy per 
nucleon is Eq — — 16.3MeV at the saturation density hq — 0.153fm^"'. 



So far, our model cannot be used quantitatively since we have not yet fixed the 
numerical values of the coupling constants. To do so one requires the model to 
reproduce the saturation density no and the binding energy per nucleon at saturation 



Note the difference between the binding energy in (infinite) nuclear matter and the 
binding energy in finite nuclei. The latter is — 8.8MeV for iron, see Eq. ( 12.64b . 

We leave it as an exercise to compute the coupling constants from the values 
( 13361 ). see problemITT] One obtains gi/{An) = 14.717 and gl/{An) = 9.537. The 
result for the density-dependent binding energy with these values for the coupling 
constants is shown in Fig. 13.21 This figure shows that there is a finite density «o 
where the binding energy is minimal. This is a basic feature of nuclear matter which 
has to be reproduced by any physically meaningful model, reflecting the properties 
of the nuclear forces. It says in particular that if you add nucleons to a large nucleus 
the density will stay approximately constant because there is a preferred distance 
between the nucleons that minimizes the energy. We have implicitly made use of 
this fact in our estimate of the nucleon number in a neutron star at the beginning of 
chapter|2] In the limit of infinite symmetric nuclear matter and ignoring the Coulomb 
forces, this density at which the binding energy is minimal is «o = 

0.153fm-^ It 

is called saturation density. At the saturation density nuclear matter is self-bound, 
i.e., it is stable at zero pressure. We have indicated in Fig. 13.21 that the minimum 
of the binding energy divides the stable density regime with positive pressure from 





(3.36) 
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the unstable regime with negative pressure. The behavior of the pressure follows 
from the thermodynamic relation P = — where E = eV is the energy and V the 
volume, which implies 

P = ,4^(|M. (3.37) 

Consequently, at the minimum of e/iig as a function of the system has zero 
pressure. Moreover, we see that a decrease in the binding energy per baryon number 
upon increasing the baryon number leads to a negative pressure. (At very small 
densities, barely visible in the plot, the energy also increases with density, i.e., P > 0. 
This is the regime where the nucleons are too far apart to feel their attraction; the 
increasing energy is then a consequence of the increasing kinetic energy.) 

In our context of compact stars, the self-boundedness of nuclear matter implies 
that nuclear matter can exist at the surface of the star where the pressure vanishes. In 
the interior, the gravitational pressure compresses the matter to densities larger than 
«o- As we see from the figure, this compressed matter, in turn, has itself positive 
pressure to counterbalance the pressure from gravity. This is the reason why the 
high-density part of the curve in Fig. l3.2l is relevant for astrophysical applications. 
We shall see in the next subsection why the Walecka model in the simple form 
discussed here cannot be trusted for densities much larger than «o and how the 
model can be improved to yield predictions for the high-density regime. 



3.1.1 Including scalar interactions 

The Walecka model accommodates important aspects of nuclear matter such as the 
existence of a saturation density whose realistic value is reproduced upon fitting the 
parameters of the model. We have already discussed on general grounds that ex- 
trapolations to high densities are uncontrolled, and thus the Walecka model (and all 
similar models of this kind) have to be improved in an interplay with experimen- 
tal observations, for example astrophysical data. But there is even a more obvious 
shortcoming of the simple version of the Walecka model discussed so far Even at 
the saturation density it fails in its prediction for the incompressibility of nuclear 
matter which is defined as 

This quantity is a measure for the stiffness of nuclear matter. In some literature, K 
is also called compression modulus or, somewhat misleadingly, "compressibility". 
To see that a large value of K corresponds to "stiff" matter, start from the usual 
thermodynamic definition for the compressibility 
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This definition says that easily compressible ("soft") matter has a small change in 
pressure upon changing the density. For the second equality we have used Eq. ( 13.371 ). 
On the other hand, from the definition ( |3.38t we obtain 



= 9nB-^ + 18 — -— , (3.40) 
on-g \nB diiB J 

where dns/dkp ~ linB/kf (see Eq. ( I3.28ab ) has been used. Now recall that in equi- 
librium, i.e., at the saturation density where the pressure vanishes, e/w^ as a function 
of «B has a minimum. 



^_ d{e/nB) 



driB 



1 / £ <9e \ 
= \--^] ■ (3.41) 



Consequently, the second term on the right-hand side of Eq. ( |3.40t vanishes at ng = 
«o and the relation between % and K becomes at saturation 

- = ^, (3.42) 

i.e., a large compressibility % corresponds to a small incompressibility K, as it 
should be. 

The calculation of the incompressibility in the given model yields K ~ 560 Me V. 
This is more than twice as much as the experimentally inferred value. Also the nu- 
cleon mass itself can be determined experimentally and compared to the prediction 
of the model. In total, there are thus four values which the model should reproduce. 
To improve the model, we add cubic and quartic scalar self-interactions of the form 

^La^-^mMigaOf-^^igaO)" (3.43) 

to the Lagrangian (13. 11 1. Besides the phenomenological need of these terms, there is 
also a theoretical reason for their presence: the model becomes renormalizable. With 
the self-interactions we have introduced two new dimensionless constants b and c 
which can be used, together with the two couplings g^, g(o to fit four experimental 
values. Namely, the two from Eq. ( 13.36b plus the incompressibility and the Landau 
mass 

K ~ 250MeV , otl = 0.83miv . (3.44) 
The Landau mass is defined as 

niL = — , (3.45) 

Vf 



where 



(3.46) 

k=kF 
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Fig. 3.3 Binding energy per nucleon as a function of density in the Walecka model, including cubic 
and quartic scalar self-interactions (solid line). The four parameters of the model are fixed to the 
saturation density, the binding energy per nucleon at the saturation density, the incompressibility, 
and the Landau mass. For comparison, the dashed line shows the result from Fig. 13.21 i.e., without 
scalar interactions. The scalar interactions account for a much softer equation of state. 



is the Fermi velocity. It is plausible that the Landau mass is experimentally more 
accessible than the mass parameter mj^ since it is an effective mass for fermions at 
the Fermi surface where all low-energy excitations are located. 

In the mean field approximation, it is easy to include the effect of the scalar 
self-interactions. The pressure becomes 

P = -\'nla^ - ^MNigaof - ^(gaO)^ + ^micb^ + , (3.47) 

with Pn defined in Eq. (|3.19t . The implicit equation for the effective nucleon mass 
(13.23b now receives contributions from the additional terms and becomes 

2 2 

= tiiN + — ^ \btnN{mN - nipjf + c{mN - 'wjv)^] ■ (3.48) 

To fit the four above mentioned values, one has to choose g^^/{A%) = 6.003, 
^ffl/(4?r) = 5.948, b = 7.950 • 10-^ and c = 6.952 • lO^^^. The numerical evalua- 
tion of the binding energy is left as an exercise. The result is plotted in Fig. 13.31 
and shows that the behavior at large densities has changed significantly compared 
to the case without scalar interactions. In particular, the lower value of the incom- 
pressibility goes along with a softer equation of state at large densities. In other 
words, if you choose a fixed binding energy on the vertical axis you find a larger 
baryon density after taking into account the scalar interactions. The matter has thus 
become easier to compress in the high-density regime, in accordance with a lower 
incompressibility. (See also discussion about stiff and soft equations of state in Sec. 

B) 
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P 1 n 


A 


Z+ 1 Z" 1 £- 




m (MeV) 


939 


1115 


1190 


1315 




1/2 


-1/2 





1 





-1 


1/2 


-1/2 


Q 


1 








1 





-1 





-1 


s 





-1 


-2 


J 


1/2 


quark content 


uud 1 udd 1 uds \ uus \ uds \ dds \ uss | dss 



Table 3.1 Mass, isospin, electric charge, strangeness, spin, and quark content for the spin- 1/2 
baryon octet. 



3.2 Hyperons 

In the interior of a compact star, densities can be as high as several times nuclear sat- 
uration density. Therefore, baryons with strangeness, hyperons, may occur (as well 
as muons). The lightest of these states are given by the baryon octet, see Table 13.1] 
It is rather straightforward to incorporate hyperons in the kind of model discussed 
above. Of course, the evaluation becomes more laborious, and the model has many 
more parameters. Let us therefore briefly discuss the model with the hyperon octet 
without going into too much detail. 

The interaction between the baryons is now extended by interactions mediated 
by the and p vector mesons. (The meson has quark content ss\ the p meson has 
the same quark content as a pion, i.e., it can be considered as an excited state of the 
pion.) The Lagrangian is 

j 

-\prp';^v+\'4Pap;- (3.49) 

Here, j runs over all eight baryons and Ta are the isospin generators. In a compact 
star, we have to require chemical equilibrium with respect to the weak interactions. 
In the case of hyperons, the conditions are 







-Me, 


Ma = Mn 


(3.50a) 






-Mo 


MxO = Mn 


(3.50b) 






+ Me, 


MsO = Mn 


(3.50c) 




= M« 


+ Mo 




(3.50d) 
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Fig. 3.4 Density fractions of baryons and leptons (and quarks, for a bag constant B = 
250MeV/fm^-^) as a function of the baryon nuinber density. The figure is taken from Ref. (T] 
where inore details about the underlying calculation can be found. For sufficiently large densities, 
hyperons and muons appear, and at densities of a few times nuclear ground state densities their 
density fractions becomes comparable to the fractions of nucleons and electrons. There is a region 
of coexistence of deconfined quark matter and baryonic matter, see Sec. l3.4l for a discussion of 
these inixed phases. The curves shown here depend on the chosen models for nuclear and quark 
matter and the value of the bag constant. 



and, including muons, iie = ji^. The conditions ( 13.501 1 all come from weak processes 
which we have already discussed, see Eqs. (I2.40l i. For example the process n ^• 
+ e + Ve, which gives rise to the condition = /i„ — /ig, can be understood 
from the elementary processes as 



We show the result of baryon and lepton density fractions in a model similar to the 
one discussed here in Fig. 13.41 

As a result of this rough discussion and the curves in the figure we learn that hy- 
perons can be included in a rather straightforward extension of the simple Walecka 
model and that hyperons do appear for sufficiently large densities. The physical 




udd — > uus + e + Ve- 



(3.51) 



Electric neutrality is given by the constraint 



(3.52) 
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reasons are that (/) they can appear because the baryon chemical potential is large 
enough to provide energies larger than their mass, (//) they do appear because (a) 
the systems seeks to acquire neutrality and does so with electrons at low densities; if 
hyperons are available, electrons in high-energy states can be replaced by hyperons 
in low-energy states and (b) the system seeks to become isospin symmetric; at low 
densities it is highly isospin asymmetric, and hyperons with nonzero isospin number 
provide a means to symmetrize the system. 



3.3 Kaon condensation 

Another possible variant of dense nuclear matter, besides the occurrence of hyper- 
ons, is the condensation of mesons. Originally, pion condensation was suggested 
|l2|. Only many years later, it was realized that kaon condensation is possible in 
compact stars [,3]. This is somewhat surprising since kaons are much heavier than 
pions and thus pion condensation seems more likely. However, in the medium, the 
effective kaon mass becomes sufficiently small to allow for a kaon condensate. 

Kaon condensation is of interest for these lectures for several reasons. Besides 
being a variant of dense matter and thus relevant for the physics of compact stars, its 
discussion requires the introduction of several important concepts in the theory of 
the strong interaction. It is thus also interesting from a fundamental point of view. 
Moreover, we shall encounter kaon condensation again later in these lectures, when 
we discuss the quark-matter relatives of the kaon, see Sec. 14.2.1] 

To explain kaon condensation, we will first have to say what a kaon is and will 
do so with the help of chiral symmetry and spontaneous breaking thereof. Then, we 
will discuss chiral perturbation theory. This is one possible method to study kaon 
condensation and has been used in the original work [[3]. For another approach, us- 
ing models similar to the above discussed Walecka model, see for instance Ref. ||4l 
and references therein. The evaluation of the chiral model has to be done numeri- 
cally, so we will more or less only be concerned with setting up and understanding 
the basic equations. As a modest goal, we will try to understand the onset of kaon 
condensation, i.e., we will show how to compute the critical baryon density at which 
there is a second-order phase transition to the kaon-condensed phase. 



3.3.1 Chiral symmetry of QCD 

Kaon condensation can be discussed in a low-energy effective theory, here chiral 
perturbation theory. This theory should describe the fundamental theory, QCD, in 
the low-energy limit. In order to construct the theory, we need to understand the 
underlying symmetries. The QCD Lagrangian is 



^QCD ^\lf{ir''D^+llYo-M)Y + ^gluons , 



(3.53) 



3.3 Kaon condensation 



41 



with the quark spinor \j/ in color, flavor, and Dirac space, the mass matrix in flavor 
space 

/m„ \ 
M= I m,/ , (3.54) 
V mJ 

and the covariant derivative = d^ — igTaA'^, where r„ = A„/2 (a = 1 , ... 8) are 
the generators of the color gauge group SU{3)c with the Gell-Mann matrices A„, 
are the corresponding gauge fields, and g is the strong coupling constant. The chem- 
ical potential /i is a diagonal matrix in flavor space. Without taking into account the 
weak interactions, each flavor is conserved and there are three independent chem- 
ical potentials. We have already seen in the previous sections that after taking into 
account weak interactions there are only two chemical potentials, one for quark 
(baryon) number, and one for electric charge. 

The purely gluonic contribution to the Lagrangian is given by 

■2gluons = ^jiv ' (3.55) 

where G"^y = d^A" - dyA'^ +gf"'"'A''^Ay with the SU{3)c structure constants is the 
gluon field strength tensor. Here we are not interested in this gluonic part, since 
we focus on the transformations of the fermion fields and the resulting symmetries 
of the Lagrangian. Also later, when we shall use QCD for explicit calculations, 
the gluonic part is negligible because we always work at very small temperatures 
compared to the quark (or baryon) chemical potential. The interactions of the quarks 
via gluon exchange, included in the covariant derivative, is of course important; in 
Sec. l4.3l this interaction will be used on the microscopic level. 
We now introduce the chirality projectors 

p 1 + 75 p 1-75 

Pr = 2 ' " 2 ■ (3.5b) 

They obey the identities 

Pr/l^Pr/l, pI/l^Pr/l. PrPl^O, Pr+Pl=1, (3.57) 

i.e., they form a complete set of orthogonal projectors. (These identities are obvious 
with = 1 and 75 = 75 .) For a physical picture, remember that, for massless quarks, 
chirality eigenstates are also eigenstates of helicity. Therefore, in this case, there is 
a one-to-one correspondence between chirality and the projection of the fermion 
momentum onto its spin. We define left- and right-handed quark spinors by 

Vr/l^Pr/lV, (3.58) 
such that \(/ = Pr\(/ + Pi\j/ = \(/r + \j/L. Then, using 



{75,7m} = 0> 



(3.59) 
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we can write the Lagrangian as 

- \(fRM\(/L - v/lMv/r + ifgiuons • (3.60) 

Let us first discuss the massless case, M = 0. In this case, separate rotations of left- 
and right-handed spinors leave the Lagrangian invariant, 

Yr ^ Xj/R, xi/L^ Wl ■ (3.61) 

Since we are interested in three quark flavors, ?„ are the nine generators of the fla- 
vor group U{3), tQ ~ 1 and = r„ (a = 1 , . . . 8). Consequently, the Lagrangian is 
invariant under U{3)l'x U {3)r. The corresponding Noether currents are 

They can be rewritten in terms of vector and axial-vector currents 

Ja.V = J^.R+Ja.L = wft^xi/, (3.63a) 
JaA ^ Ja.R-Ja.L = Wtt.YsW ■ (3.63b) 

To see this, note that PrJs — Pr and PlJs = —Pl- In QCD the singlet axial-vector 
current is in general not conserved, 

Vo'^ = -f3^MvGr, (3.64) 

where G^^ = \e^^°PGap is the dual field strength tensor This is referred to as the 
axial anomaly. We are left with the symmetry group 5?/ (3)« x SU{3)l xU{l)v- The 
vector symmetry U{l)v corresponds to baryon number conservation and is therefore 
also denoted as U{1)b- The flavor symmetry group SU{3)r x SU{3)l is referred to 
as chiral symmetry. As we can see from Eq. ( 13.60b . nonzero masses break the chiral 
symmetry ejcp/ic/fZy. They do not break the t/(I)y symmetry, and for the special case 
nil, =md — nis the subgroup SU{3)r+l of simultaneous R and L rotations remains a 
symmetry of the Lagrangian. 

Spontaneous breaking of chiral symmetry is realized by a chiral condensate of 
the form (i/ZlI/Zr). This is analogous to spontaneous symmetry breaking in simple 
models such as theory (see for instance the discussion of Bose-Einstein conden- 
sation in appendix lA.lb . or in a superconductor, or in the Higgs mechanism. The 
chiral condensate is only invariant under simultaneous right- and left-handed rota- 
tions, i.e., the symmetry breaking pattern is 

G = SU{3)r X SU{3)l -^H = SU{3)r+l. (3.65) 

As a comparison, in 0'* theory with a complex scalar field (/), we have G = ?7(1), H = 
1, which gives rise to the familiar "Mexican hat" potential with a negative quadratic 
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and a positive quartic term in 1 1 . Spontaneous breaking of a global symmetry goes 
along with massless Goldstone bosons. In the Mexican hat, there is one massless 
excitation along the bottom of the Mexican hat, given by the angular component of 
the order parameter (while the radial component corresponds to a massive mode). 
Here, the bottom of the Mexican hat is not just a one-dimensional line. It is rather 
given by the coset space G/H (which is simply t/(l) in 0^ theory). This space has 
dimG — dim// generators. Consequently, with dimG = 8 + 8 = 16 and dim// = 8, 
there are 8 Goldstone modes. They are described by the SU (3) matrix 



U = , 



(3.66) 



with the pion decay constant fji ~ 93 MeV. The meson fields Oa of the Goldstone 
octet are usually reparametrized as 



\/2 ^/6 



tO 



V 



K 



7t" rj 



K+ 

/ro 



\ 



(3.67) 



3V 



Since the rows (columns) of this matrix carry left-handed flavor (right-handed anti- 
flavor) labels, it is easy to read off the quark content of the various mesons, e.g., 
~ SU, ~ du etc. According to its chiral structure, the chiral matrix transforms 
under a transformation g = (gL,gR) G G as 



U^gLUgl 



(3.68) 



3.3.2 Chiral Lagrangian 

In the (unrealistic) case of vanishing quark masses, the chiral symmetry is an exact 
symmetry and the Goldstone bosons are exactly massless. Exploiting the analogy 
to the Mexican hat potential, this means that the bottom of the Mexican hat is truly 
flat. Quark masses break the chiral symmetry explicitly. However, if the masses 
are small compared to the characteristic scale of chiral symmetry breaking A ^ 
Anfji^l GeV we can still consider the chiral symmetry as approximate. The bottom 
of the Mexican hat then gets distorted on a scale small compared to the deepness 
of the potential, and the Goldstone bosons acquire small masses. In this case it is 
more appropriate to speak of pseudo-Goldstone bosons. One might still hope to 
describe the system at low energies by an effective theory which is built on the 
underlying chiral symmetry, although this symmetry is strictly speaking broken. 
The mass matrix M, now nonvanishing, is required to transform just as the chiral 
field U, i.e., 

M^gLMgl. (3.69) 
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We require the chiral Lagrangian to be invariant under G. The kinetic term and the 
mass term of the resulting effective theory are 



-cTx[M^lJ ^MU\ 



(3.70) 



where the trace is taken over flavor space. The two constants f-^^ and c have to be 
fitted to experimental values, similarly to the constants of the Walecka model. In 
principle, higher order terms in U are allowed but shall be neglected here. Note that 
the Goldstone fields themselves appear in the exponent of the field C/, i.e., they are 
already present to all orders. 

In the context of compact stars, we do not only want to describe isolated mesons. 
We also need to include baryons and their interactions. The baryon octet fields are 
given by the matrix 



B 



\ 



A 

7e 



E+ 
V2 \/6 



2^ 



(3.71) 



which includes the proton p, the neutron «, and the hyperons from Table 13. II A 
simple way to understand the structure of this matrix is as follows. Consider the 
baryons as composed of a diquark and a quark. The diquarks form an antitriplet, i.e., 
one can think of the columns of the matrix as labelled by (t7, c/, s) which corresponds 
to the quark content {ds,us,ud). Then the rows are simply labelled by the flavors 
in the fundamental representation {u,d,s), and one easily checks that this yields the 
quark content of the baryons as given in Table 13.1] 
The free baryon Lagrangian is 



^B=Ti[B{ifd^-mB)B], 



(3.72) 



where — 1.2GeV is the SU{3)l x SU{3)r symmetric baryon mass. To write 
down the interaction between baryons and the mesons it is convenient to decompose 
the chiral field into left- and right-handed fields, 



where, without loss of generality, we may choose 



(3.73) 



(3.74) 



such that 



(3.75) 
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We now add the meson-baryon interaction terms 13] |5| (see also chapter 7 of Ref. 
||6| for more details), 

= iTi[Br^[J^,B]]+DTT[BY^Y5{4,B}]+FTi[BY^Y5[4,B]] 
+aiTr[B'{^M^+^'^M^^'^)B] +fl2Tr[B^B(<^M<^ + <^^M'''^^)] 
+a3TT[B''B]TT[M^^+M\^^)^] , (3.76) 

with the additional constants D, F, ai, 02, 03, and the vector and axial-vector cur- 
rents 

J^^\{^^d^^+^d^^^), (3.77a) 
J^^^i^^d^^-^d^^^). (3.77b) 

The Lagrangian is an expansion in M/A and d /A with the scale of chiral symmetry 
breaking A . Higher order terms in these parameters are omitted. In summary, we 
have the Lagrangian 

= ^f/+^B + ^/. (3.78) 

Later we shall also add electron and muon contributions, but they are simple and we 
ignore them for now to keep the notation brief. 



3.3.3 Kaon-nucleon matter 

Since we expect (charged) kaon condensation in a compact star rather than any other 
meson condensation (possibly there is pion condensation) let us for simplicity drop 
all mesons other than the kaons. We can then write 

f/ =e'2 = cos e + 'Sing, (3.79) 

with 

7 / 04-!05\ 

Q=Y.^o^a= <p6-i<h\, (3.80) 

"=4 \04 + '^5 06 + '</>7 ) 

where we have absorbed /jt into the fields 0a = Oa/ Ik such that the 0a's are dimen- 
sionless. We can now compute a simple expression for the matrix U . To this end we 
first verify by explicit matrix multiplication 

Q' = fQ, (3.81) 



where 



(3.82) 



46 



3 Basic models and properties of dense nuclear matter 



From Eq. ( 13.811 ) we obtain (for instance via complete induction) 



for all « > 1 , which can be used to compute 
cosg — 1 

= 1 



2(«-l)^2 









4! 


— fi - 


- COS(/)) 



(3.83) 



4! 



(3.84) 



and 



sinG = G( 



Q- 



5! 

3! 5! 



— sinffl . 



2n2 



5! 



(3.85) 



As a further simplification let us now drop the neutral kaon fields, = = 0, 
because we expect charged kaon condensation. Then, from Eqs. ( 13.84b and (I3.85l l 
we obtain 

' """-^ ^ • sin^ > 



U = 



V 



cos (j) 


. 04 + '05 



i- 
1 



sin0 






cos (p 



(3.86) 



Now we interpret the fields ^4 5 as vacuum expectation values, 04 5 (04.5), and 
neglect the fluctuations around this background. The general procedure to describe 
Bose-Einstein condensation, including fluctuations, is explained in appendix lA. 1 I for 
the 0'* model. The condensates are assumed to be constant in space and to have the 
time dependence 0(r,x) (j)e^'^'^', i.e., our ansatz is 



{K-) = (04)+/(05> 
{K+) = (04) -/(05) 



: 0e-'>"^' , 



(3.87a) 
(3.87b) 



The real, constant (i.e., space-time independent) value of has to be determined 
later from minimizing the free energy; jiK plays the role of a kaon chemical po- 
tential, as we shall see more explicitly below. More precisely, jik is the chemical 
potential for K while —jJ-K is the chemical potential for K^. We thus arrive at 



U = 



cos ie''^'^' sin \ 
1 

/g-'MK'sin^ COS0 / 



(3.88) 
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We are now prepared to evaluate from Eq. ( 13.701 ). We shall neglect the masses 
of the up and down quarks such that M ~ diag(0,0,mi). We also define the kaon 
mass 

m\ = —^. (3.89) 

J-K 

Rather than c, we shall later use the kaon mass w/f ~ 494 MeV as a parameter of the 
model. This yields 

= -Vi(^) (3.90) 

with the tree-level potential 

V{<p) = -^sin2,/)+m|/2(l -cos(/)), (3.91) 

where we have subtracted the constant vacuum contribution y((/) = 0). This potential 
contains the kaon condensate to all orders. We shall work with this expression below, 
but it is instructive to expand it up to fourth order in 0, 



^ '^^-^ {M f + M_^(/,0)4 . (3.92) 



24/2 

This is the familiar expression from a <^'^ model for the free energy of a Bose 
condensate with chemical potential /i/f, mass niK, and effective coupling (4/i|^ — 
m\) I (6/^), see for instance Eq. ( IA.18I ) in the appendix. As expected, condensation 
occurs for yi\ > nvj^ because in this case the quadratic term is negative and the quar- 
tic term positive, i.e., we have recovered the Mexican hat potential (where we have 
already picked one direction since is real). 

For the baryonic Lagrangian we only keep the lightest baryons, the proton and 
the neutron. From Eq. ( I3.72l i we thus obtain 

= Piif'dyi - ruB + /Mp)/5 + n{ifdf, - mg + flin)n , (3.93) 

where we have added the proton and neutron chemical potentials jXp and /i„. For the 
interaction terms we need 



/ cos(^/2) ;e'>"^' sin(^/2) \ 

1 

V /e-Wsin(^/2) cos(^/2) / 



(3.94) 



which obviously fulfills = U. By inserting this into Eqs. ( 13.77b we see that the 
spatial components of the currents vanish, Jv = Ja = (since there is no spatial 
dependence in the condensate), and the temporal components are 
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/-I 00 

4 = iliKsm^i<l>/2) 00 
V 1 

/ -e'>'^' \ 
J^^ iliKCOs{(j)/2)sm{(j)/2) i . (3.95a) 

y^-W Q / 

Hence the various nonzero terms needed for ^ in Eq. ( |3.76t become 

iTT[Byo[4,B]] ^ HK{2p^p + n^n)sm^{<j)/2), (3.96a) 

aiTi[B'^{^M^+^'^M^^^)B] = -2aim,p^ psm^{^/2) , (3.96b) 

a2Tr[B^B{^M^ + ^''M^ ^^)] 2fl2»J.v(/p + nMcos^(0/2) , (3.96c) 

a3Tr[B^B]Ti[M^^+M\^''f] = 2a3ms{p^ p + nn)[l - 2sm^{^/2)] i3.96d) 

It is left as an exercise to verify these resuhs. Inserting this into Eq. ( 13.761 ). and 
putting together the contributions from the chiral field, the baryons, and the interac- 
tions between them, the total Lagrangian can be written as 

if = -y (0) +p[iY^d^ -mB + f\^p + pi;)]p 

+n[iy^d^ - MB + y + lK,)> ■ (3.97) 

Similar to the Walecka model in Sec. 13.11 the effect of the kaon condensate on the 
nucleons can be absorbed into an effective chemical potential. In a slightly different 
notation than in Sec. 13. II (where pi was absorbed into ^*), we have kept the actual 
thermodynamic chemical potentials separate, and we have 

jU* = 2(fl2 + fl3)m, + [2jU/f-2(fli +fl2 + 2fl3)OT.,]sin^(^/2), (3.98a) 
PL* = 2(fl2 + fl3)'Wi + [Ai/f-2(fl2 + 2fl3)m.v]sin2(^/2). (3.98b) 

We can now, analogously to Sec. 13.11 evaluate the partition function at T = to 
obtain the thermodynamic potential density Q. = —T/V InZ which can be written 
as 

Q=V{(j))+eB-{M*,+l^„)n„-{M*p + l^p)np, (3.99) 
with the nucleon number densities n„ and rip, and the nucleon energy density 

£b = 2 f -^Jk^Ws@{kF,i-k), (3.100) 

i=p.,n-' y^^l 

where kp^t are the respective Fermi momenta. 

Before adding the lepton contribution we need to find the relations between the 
various chemical potentials through the conditions of chemical equilibrium. The 
leptonic processes including nucleons are 
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■p + l+Vi, p + e^n + V(. (3.101) 



Here, £ = e^ji can either be an electron or a muon. We also have the purely leptonic 
processes, 

e ^ jd + V^ + Ve, n ^e + Ve + V^, 

and the processes involving kaons. 



(3.102) 



n^p + K', e^K' + Ve. (3.103) 
These processes lead to the independent conditions 

lie ^ llK ^ ll^x , lln= pip + lie ■ (3.104) 

The system is thus characterized by two independent chemical potentials, say jif, and 
jj.,, . We implement the constraint /i„ — pLp + pie by rewriting the terms containing the 
nucleon chemical potentials in the potential ( 13.991 ) as /i„n„ + pLpiip = llnng — l-ienp- 
Since we want to work at fixed nB = n„+np, we perform a Legendre transformation 
of £2 with respect to the variables jj.,, and = — hb. This amounts to adding the 
term /X„«b to £2 which yields the relevant free energy for the baryons and the kaon 
condensate, 

QB,K = V{^) + eB- [{li;-lie)Xp + {l -Xp)li;,] IIB. (3.105) 

Here we have introduced the proton fraction 

Xn=—, (3.106) 

which has to be determined dynamically from minimizing the free energy. We can 
now add the lepton contributions to arrive at 

Qb,k,I: = £iB.K + Ee- l^e'le + - m]^){.£li - jU^-Hm) ' (3.107) 

where we have used lie = li^, where are the lepton energy densities (£ ~ e,li), 
and where 

ne = f^, n^='-^-^—^— (3.108) 

are the corresponding lepton number densities. The step function in the muon con- 
tribution accounts for the fact that muons only appear if /i^ is larger than their mass 
TO^ = 106 MeV. On the relevant energy scale, electrons are massless to a very good 
approximation and thus are present for any nonzero /i^. 

For a given baryon number «b, the variables of CIb.k/ ^re the proton fraction Xp, 
the kaon condensate 0, and the chemical potential for (negative) electric charge /i^. 
They are determined by minimizing the free energy with respect to Xp and and by 
requiring charge neutrality. 
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B.K.i 



B.K.i 



dQ, 



'.K.f: 



dXp d(j) 



= 0. 



(3.109) 



It is straightforward to compute the various derivatives, and after a few lines of 
algebra the result can be written as 



1 



«sCOs2(0/2) dxp 



— 2aimstcLn^{(j) /2) , 



(3.110a) 



= cos - —I 
All 



riB 



-(1- 



I — (fliXp + a2 + 2a^)ms 



= fllle?<in- (p - riB [xpCos2(0/2) - sin2(0/2)] 



, (3.110b) 



(3.110c) 



The second equation has been obtained after dividing both sides by sin^. This 
means that (/) = is always a solution and Eq. ( 13.1 lObI ) is only valid for nonvan- 
ishing condensates. In the third equation we recover the various contributions to the 
electric charge density: the first term on the right-hand side is the pure contribution 
from the kaon condensate. It gives a positive contribution to the negative charge 
density for jie > 0. The second term on the right-hand side arises from the nucleons 
and their interactions with the kaon condensate. Only for 0=0 does it give the pure 
proton contribution —Hp = —XpiiB- Finally, the other two terms are the expected 
contributions from the leptons. 

The onset of kaon condensation can be determined by setting = in all three 
equations. This yields three equations which can be solved for x^, /i^ , and n^, where 
rfg is the critical density beyond which there is a condensate and Xp, /ig the values 
of the proton fraction and the charge chemical potential at this density. Since Eq. 
( 13.1 lObl i is only valid for 7^ 0, this has to be understood as approaching «g from 
above. 

We leave the numerical evaluation of the critical density and the general evalu- 
ation for all yiB as an exercise, see problem [3. 3 1 An important modification, which 
we have neglected for simplicity, has to be taken into account for this evaluation. 
Namely, the energy density £b has to be modified due to interactions among nu- 
cleons. It is beyond the scope of these lectures to derive this modification, see Ref. 
||5l and references therein for more details. Here we simply quote this modification 
which is needed in order to get physically sensible results. One needs to use an 
expansion of £b around symmetric nuclear matter Xp ~ 1 /2 of the form 



£b £0 + «b(1 - 2xpf-S{u) , 



no 



(3.111) 



Here Eq is the energy density of symmetric nuclear matter, whose form is not rel- 
evant because we only need the derivative of Eb with respect to Xp. The nuclear 
saturation density is denoted by «o, and 
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Fig. 3.5 Density fractions of neutrons («), protons (p), electrons (e ), muons (/i ), and the kaon 
condensate (K^) froin Eqs. B.UOt with modified as given in Eqs. 13. UU . )3.112t . The param- 
eters are (see Ref 0) a^m, = -67MeV, aom, = 134 MeV, a^m, = -222MeV, = 106MeV, 
= 93 MeV, niK = 494 MeV, mg = 1200 MeV. We see that for this parameter choice the onset of 
kaon condensation is at about three times nuclear saturation density, itg ~ 3.2;io- 



5(«) = (22/^-1) 



,3(37rW2)2/3 ^ 



u^l^^F{u) +SoF{ii) 



(3.112) 



is the nuclear symmetry energy (see Ref. (|2l for a discussion of the nuclear symme- 
try energy in the context of the maximum mass of neutron stars). For the numerical 
evaluation shown in Fig. 13.51 the nuclear symmetry energy at the saturation point 
So = 30MeV has been used, as well as the function F{u) = u. See caption of the 
figure for the choice of the other parameters. 



3.4 From hadronic to quark phases: possibiUty of a mixed phase 

We have already mentioned the possibility of a hybrid star, i.e., a star with a quark 
matter core surrounded by nuclear matter. How does the interface between these 
two phases look? Is it a sharp interface or is there a shell in a hybrid star where the 
hadronic and quark phases coexist in a mixed phase? If the former is true, there will 
be a jump in the density profile of the star, while the latter allows for a continuous 
change in density. 

Mixed phases are a very general phenomenon. In the context of compact stars, 
not only the mixed hadronic/quark matter phase is of relevance. Also in the inner 
crust of a hybrid or neutron star one may find mixed phases. There one expects a 
neutron superfluid coexisting with a lattice of ions, i.e., a mixed phase of neutron 
matter and nuclei. In these lectures, we shall not discuss the properties of the crust 
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Fig. 3.6 Illustration of the possibility of a inixed phase. The pressure P of two phases A and B 
is given by the respective curves as a function of a chemical potential /i. Note that dP/d^l has to 
increase with increasing /i (increasing /i cannot lead to a decrease of the corresponding charge; this 
would lead to an instability). Suppose is the electric charge chemical potential and we require 
charge neutrality. Then the squares mark the points at which a given phase is charge neutral. The 
circle in the left panel marks a point where the two phases have equal pressure and opposite charge. 
Since this point has higher pressure than either of the squares, this is the ground state (neglecting 
surface tension and Coulomb energy). In this state, phase A and B coexist and occupy different 
volume fractions, detennined by the different slopes of the curves. In the right panel, there is no 
point where both phases have equal pressure and opposite charges. Therefore, the square on the 
curve B is the ground state. 



of a compact star in detail. See Sec. l6.2l for a brief discussion and Ref. fS] for an 
extensive review. Other examples of mixed phases in different systems are liquid- 
gas mixtures or simply a solid, which is a mixture of an electron gas and nuclear 
matter (sitting in the lattice of ions). 

In Fig. l3.6l the possibility of a mixed phase is illustrated. We see that the condition 
of charge neutrality plays an important role here. It is important that in a compact 
star charge neutrality is required globally, not locally. In other words, certain regions 
in the star may very well have a nonzero electric charge as long as other regions have 
opposite charge to ensure an overall vanishing charge. 

It is plausible that such a mixed phase will have a crystalline structure. For in- 
stance, one phase may form spheres sitting at the points of a lattice which is im- 
mersed in the other phase. Other possibilities are rods or slabs f9l, such that the 
mixed phase looks like spaghetti or lasagna, wherefore astrophysicists have termed 
such phases nuclear pasta. In any case, if a mixed phase is possible because of a 
general argument such as given in Fig. 13.61 this does not mean that it is indeed re- 
alized. One has to take into account Coulomb forces (which seek to break charged 
regions into smaller regions) and surface forces (which seek to minimize the surface 
and thus work in the opposite direction). We shall not discuss these forces quantita- 
tively but rather give some general arguments about mixed phases. 

We start from the simple picture that at small quark density (or quark chemical 
potential pC) the hadronic phase is preferred and that there is a first-order phase 
transition to the quark matter phase at some critical chemical potential. The question 
is whether there is a mixed phase between these two pure phases. The pressures of 
the two phases and depend on the quark chemical potential and 
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the charge chemical potential He (we work at zero temperature). Phase coexistence 
is possible when the pressures of the two phases are equal, 

Pl,{fl,l^e)=Pg{fi,l^e). (3.113) 

Now suppose the neutrality condition were local (which it isn't in our context). Then 
the charge must vanish in each phase separately, 

Qi,{li,lie)^Q,{fl,lie)=0. (3.114) 

These two conditions yield fie for each phase separately as a function of ji, ^^{ll) 
and Consequently, the condition of equal pressure, 

P,,(M,Mi'(M))=^,(Ai,M.^(M)), (3.115) 

yields a unique /i. Only at this /i do the phases coexist. This amounts to a sharp 
interface at a given value for the pressure, where on both sides of the interface the 
pure hadronic and the pure quark phases exist with different densities, i.e., there is 
a density jump in the profile of the star. 

Now we impose the weaker (and realistic) condition of global charge neutrality. 
This means that in any mixed phase only the total charge has to vanish. We denote 
the quark volume fraction by 

yh + Vq 

where Vq and V/, are the volumes occupied by the quark and hadron phases, respec- 
tively. Then, neutrality reads 

[\-Xc)Qh{H,He)+XqQci{ll,He)^Q- (3.117) 

This yields a function pie{Xq,pt^) which is then inserted into the condition of equal 
pressure, 

Ph{lJ-,IJ-e{Xq,l^))=Pq{lJ-:MXq:l^))- (3.118) 

The result is a chemical potential as a function of Xq, p^iXq)- Thus there is a finite 
interval on the jU-axis where a mixed phase is possible. We see that the looser con- 
dition of global charge neutrality allows for a shell with a mixed phase in a hybrid 
star. These formal arguments become more transparent in a geometric picture, see 
Fig.lll] 

We shall not go into the details of an explicit calculation of the quark/hadron 
mixed phase because, even neglecting surface tension and Coulomb energy, this 
calculation eventually has to be performed numerically. Instead we show the result 
of such a calculation LIO] in Fig. |3.8l (cf. also Fig. |3.4| where we have akeady seen a 
mixed phase). One recovers the (projection of the) topology of Fig. |377l in Fig. 13.81 
The figure shows the mixed phase being the preferred phase in a certain interval 
without taking into account Coulomb energy and surface energy. In the complete 
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Fig. 3.7 Schematic picture of a hadron-quark mixed phase in a finite interval of /i. Left panel: 
the pressures of the two phases define two surfaces parametrized by /i and ^u^. The intersection 
of the two surfaces forms a line where coexistence of the two phases is possible. Right panel: the 
neutrality condition for each of the phases defines a curve in the ^l-fl^ plane, and thus a curve on 
the respective surfaces (for illustrative purposes let the charge be nonzero - denoted by c - since 
for zero charge there would have to be a valley of the pressure). A mixed phase may exist from A 
(where Xq = 0) to B (where Xq = 1). provided that, for a given /i, the pressure on this line is larger 
than the respective pressure on the neutrality curves of each phase. In this segment none of the 
phases is neutral separately, but they may combine to a globally neutral phase. Note that the extra 
direction /i^ is crucial to have a finite segment along the axis where a mixed phase is possible. 
If the mixed phase is realized, the arrows indicate the ground state for increasing values of (the 
pressure also has to increase along this line). 



calculation one finds that a relatively small surface energy is needed to destroy the 
mixed phase. It thus appears unlikely that a mixed phase of quarks and hadrons 
exists in a hybrid star. 



Problems 

3.1. Binding energy and saturation density in the Walecka model 

Solve Eq. (13.23b at zero temperature numerically for different values of the baryon 
density. Use the solution to compute the binding energy per nucleon and check 
that the values ( 13.361 ) are obtained upon using the values of the coupling constants 
g^/(47r) = 14.717, g^/(47r) = 9.537. In other words, reproduce the results from 
Fig. 13.21 If you are a bit more ambitious you can also do it the other way around; set 
up and solve the two equations that are needed to determine the coupling constants 
from the conditions (13.361 1. 
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H (MeV) 

Fig. 3.8 Figure from Ref. 1101 showing the transition from nuclear matter (NM) to a mixed phase 
(mix) to a quark matter phase (CFL) (color-flavor locking (CFL) is explained in Sec. |4.2) . In the 
mixed phase, fie is lowered in order to make the nuclear phase positively charged and the CFL 
phase negatively charged. Taking into account Coulomb energy and surface energy shows that 
the IX interval for the mixed phase shrinks with increasing surface tension (J until it completely 
disappears for c > 40 MeV / fm" . The exact value of C7 is not known but it is likely to be larger than 
that limit value such that a mixed phase appears unlikely. The limit value does not depend much 
on whether the mixed phase has spheres, rods, or slabs. 



3.2. Walecka model with scalar interactions 

Reproduce the result of Fi g . 1 3 . 3 1 numeric ally. 

3.3. Onset of kaon condensation 

Solve equations ( 13.1101 ) - with the modifications given in Eqs. ( 13. 11 11 1 and ( I3.112l i 
- numerically to determine the density fractions of nucleons, kaons, and leptons at 
r = as a function of baryon density. In particular, compute the critical baryon 
density for the onset of kaon condensation. See caption of Fig. l3.5l for the values of 
the parameters and compare your result to the plot in this figure. 
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Chapter 4 

Superconductivity and superfluidity in a 
compact star 



In our discussion of interacting nuclear matter we have so far ignored a very im- 
portant physical effect. We have not included the possibility of superfluidity and/or 
superconductivity, although we have briefly mentioned the effect of superconduc- 
tivity on the equation of state of quark matter, see Sec. 12.31 In the following, we 
shall discuss these effects in more detail. But first let us recapitulate what super- 
conductivity is. Once we have introduced the basic concept we shall see that it may 
appear in several variants in a compact star And we will see that it is crucial for the 
understanding of transport properties of dense matter And the transport properties 
of dense matter, in turn, are related to the phenomenology of the star. 

Consider a system of fermions at zero temperature with chemical potential ji and 
free energy 

Q=E-plN. (4.1) 

Now first suppose the fermions are noninteracting. Then, adding a fermion with en- 
ergy /X, i.e., at the Fermi surface, leaves the free energy Q. obviously unchanged; the 
energy E is increased by jU, but the second term subtracts the same amount since 
we add N = I fermion. Now let us switch on an arbitrarily small attractive interac- 
tion between the fermions. Then, by adding two fermions at the Fermi surface, we 
can actually lower the free energy because the attractive interaction will lead to an 
energy gain from the binding energy. Therefore, the Fermi surface we have started 
with is unstable. A new ground state is formed in which pairs of fermions are cre- 
ated at the Fermi surface. Since two fermions formally can be viewed as a boson, 
these fermion pairs will form a Bose condensate^ This formation of a condensate 



In fact, the fermions are correlated in momentum space, not in real space. Consequently, in the 
weak-coupling limit, the fermion pairs are not spatially separated bosons. The typical size of a 
pair is rather lai'ger than the mean distance between fermions. Therefore, one apparently has to be 
careful to describe the pairs as bosons. However, recent experiments with cold fermionic atoms 
show that there is no phase transition between the weak-coupling limit (where the pairs are wide 
spread) and the strong-coupling limit (where the pairs are actual difermions, i.e., bosons). This is 
the so called BCS-BEC crossover. This observation suggests in particular that it is not too bad to 
think of the fermion pairs as bosons even in the weak-coupling limit. 



57 



58 



4 Superconductivity and superfluidity in a compact star 






2A 







Fig. 4.1 Left panel: particle and hole excitations (solid lines) in a system of noninteracting ultrarel- 
ativistic fermions with chemical potential ^l. The dashed lines are the antiparticle and antiparticle 
hole excitations. Right panel: quasiparticle excitations after switching on small interactions which, 
via Cooper's Theorem, give rise to an energy gap A according to Eq. j4.2t . here chosen to be 
A = 0.5/i. What were pure particle and pure hole excitations in the left panel have now become 
momentum-dependent mixtures of particles and holes. 



of fermion pairs due to an arbitrarily small interaction is called Cooper's Theorem 
and the fermion pairs are called Cooper pairs. 

This mechanism is completely general, i.e., it holds for arbitrary fermions with 
a Fermi surface as long as their interaction is attractive. It holds for electrons in a 
usual superconductor, i.e., a metal or alloy, for ^He atoms in superfluid helium, for 
fermionic atoms in an optical trap etc. In our context, it can be applied to protons, 
neutrons, and quarks. Anticipating that the Cooper mechanism leads to superfluidity 
for neutral fermions and to superconductivity for charged fermions, we thus expect 
(;) neutron superfluidity, (/;) proton superconductivity, and (//;) quark superconduc- 
tivity to be in principle possible in a compact star. Quarks are of course a bit more 
complicated since they not only carry electric charge but also color charge. There- 
fore, we need to make more precise what we mean by quark superconductivity, see 
Sec. 1121 

Let us first stay on a very general level and discuss the basic consequences of 
Cooper pairing. A Cooper pair is held together by a sort of "binding energy" (al- 
though it is not a bound state), i.e., one needs a finite amount of energy to break a 
pair. Consequently, the single-particle dispersion relation acquires an energy gap A, 

e, = ^{E,-^f- + A^, (4.2) 

with Zi^ = \/W+n? as in the previous chapters. One might think that does not 
reproduce the usual dispersion Ej^ — pi for a vanishing gap, rather £<; — > \Ej^ — 
This is no contradiction after taking into account the fermion hole excitations, 
such that in the ungapped system e^. = ±(£'<: — /i) to which the A = Q limit of 
= ± ^J {Eli — nY -\-A'^ is indeed equivalent. The excitation described by Eq. (I4.2l i 
is also called quasiparticle since it contains the interaction of the original particles 
in an effective way. To excite a quasifermion in a superconductor, a finite amount of 
energy is needed, while a fermion at the Fermi surface of a noninteracting system 
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can be excited by an infinitesimally small energy, see Fig. 14.11 The energy gap in 
the dispersion relation is responsible for most of the phenomenological properties 
of a superconductor. For instance, it gives rise to the frictionless charge transport in 
an electronic superconductor, since (sufficiently low energy) scattering of electrons 
off phonons cannot excite a single-electron state. Or, in the context of superfluid- 
ity, the energy gap explains the frictionless flow in the same way. For quantitative 
predictions it is thus crucial to compute the magnitude of A . We shall perform this 
calculation within perturbative QCD for quarks in Sec. 14.31 

The energy gap is in general a temperature-dependent quantity. It typically de- 
creases with temperature and becomes zero at and above a certain critical temper- 
ature r^. This critical temperature indicates the phase transition from the supercon- 
ducting to the non-superconducting phase, as we shall demonstrate with the discon- 
tinuity of the specific heat in the following section. Since the onset of supercon- 
ductivity or superfluidity is a phase transition, there must be a symmetry which is 
spontaneously broken below the critical temperature. In particular for quark matter, 
the symmetry breaking pattern is very useful to characterize the superconductor, see 
Sec.lU 



4.1 Specific heat for isotropic and anisotropic superconductors 

As a example of the effect of A let us compute the specific heat of a superconductor! 
The specific heat is easy to compute and shows characteristic features of a super- 
conductor. We start from the free energy of a superconductor made of fermions with 
two degenerate (spin-) degrees of freedom, 

where the quasiparticle energy is given by Eq. ( 14.2b . We shall, for simplicity, 
consider massless fermions, = k. The entropy (density) is given by the derivative 
with respect to the temperature (with respect to the explicit temperature dependence 
only, there is also an implicit temperature dependence in A) 

with the Fermi distribution 

^^ = ^- ^'-'^ 
To derive Eq. (I4.4l i one uses the identities 



^ More precisely, here we compute the fermionic contribution to the specific heat. There may be 
fight Goldstone modes which dominate the specific heat at small temperatures. In this section we 
ignore such modes for the purpose of illustrating the effect of the fermionic energy gap. 
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y =ln(l-/,)-ln/,, ln(l+e-^*/^) =-ln(l-/,). (4.6) 
From the entropy we then compute the specific heat (at constant volume) 

cv = T^=2 -—^ek^. (4.7) 
For the temperature dependence of the gap we assume the following simple form, 



A{T)=©{T,:-T)Ao^l-^, (4.8) 

such that the zero-temperature gap is Aq, the gap approaches zero at T = Tc and 
vanishes for all temperatures larger than T^. Then, for T < Tc we have 

dA _ AlT _ dSk _ T Al dfk_l e^^l^ ( ej A^ 



dT T^A dT E^T} dT ek{e^klT + if\r- T} 
and consequently 



(4.9) 



cv^2 -\^ + ^\. (4.10) 



r2 T} 



We are only interested in temperatures much smaller than the chemical potential, 
T jj.. Then, the main contribution comes from the Fermi surface, and we can ap- 
proximate dkk^ j^^dk. We introduce the new variable x = [k — ii)/T, and define 



<P=^. (4.11) 



This yields 



where we have approximated the lower boundary by —jx/T ~ —00 and have used 
that the integrand is even in x (which gives rise to the new integration boundaries 
[0,°o] and a factor 2). We have not yet performed the 9 integral since we shall al- 
low for anisotropic gaps. From this general expression we easily get the limit of a 
vanishing gap, (p = Aq ^ 0, i.e., the result for the non-superconducting state, 

dx-, 7- = ^- (4.13) 



n- Jo 1 + coshx 
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Before evaluating the specific heat in the superconducting phase at small tempera- 
tures, let us discuss the behavior of cy at the critical temperature. This is best done 
by looking at Eq. ( 14.101 ). Approaching T(. from above, cy is given by setting Aq and 
A (T) (appearing in e^^) to zero in that equation. In the superconducting phase, ap- 
proaching Tc from below, we only set A in ei^ to zero. Consequently, at 7^ there is a 
jump in the specific heat which is given by 



Acv = 2 



dx- 



1 



2, ,2 



K^Tc Jo l+coshjc n^Tc 



(4.14) 



where we have assumed the gap to be isotropic. This jump is a typical signature for 
a second-order phase transition, since the specific heat is the second derivative of 
the thermodynamic potential. 

Next we evaluate Eq. (|4.12t for temperatures much smaller than the gap, i.e., in 
the Umit cp ^ °°. First we consider an isotropic gap. We can approximate 



Consequently, 



.,-2 



(4.15) 



cv 



dxx^e ^ 



r3/2 



(4.16) 



where we used 



dxx^e ^ = / — 



dxe ^ = ^ 



1/2 



(4.17) 



The main result is that the specific heat is exponentially suppressed by the factor 
e^f = e^^/^ for temperatures much smaller than the gap. The suppression of the 
specific heat in a superconductor provides a good example to get some intuition for 
the properties of superconductors. To this end, note that the specific heat is a mea- 
sure of how many degrees of freedom are available to store heat. A large number of 
degrees of freedom means a lot of "storage room" and thus a large specific heat. A 
small specific heat, such as for a superconductor at sufficiently small temperature, 
thus means there are very few states available. This is a direct consequence of the 
energy gap which obviously leads to a region in the energy spectrum with no al- 
lowed states. Only by increasing the temperature does the exponential suppression 
disappear because temperature provides the energy to populate states above the gap 
which in turn are then available to store thermal energy. 
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Next let us assume an anisotropic gap of tiie form 

A^Asine. (4.18) 

In a compact star, anisotropic gaps may be realized in neutron superfluidity and 
possibly for quark superconductivity. The reason is very different in the two kinds of 
matter: at large density, the i-wave interactions between neutrons become repulsive 
and thus only interactions in the /7-wave channel can lead to superfluidity (this is 
in contrast to protons which do form i-wave superconductors). In the case of quark 
matter, anisotropic gaps may occur due to a mismatch in Fermi momenta of the 
quarks that form Cooper pairs; anisotropics then arise either because the mismatch 
allows only for pairing in certain directions in momentum space or because pairing 
occurs in the spin-one channel which does not suffer from the mismatch. In either 
case, there are several possibilities for the specific form of the angular dependence 
of the gap and it is not entirely clear which one is reaUzed in the relevant density 
regime. For more details see for instance Refs. |[Tl|2l for nuclear matter and Ref. Q 
for quark matter. 

With 9 being the angle between the momentum and the z-axis, the form (14.181 1 
implies point-like nodes of the gap function at the north and south pole of the Fermi 
sphere. In other words, although there is a nonzero order parameter for superfluidity, 
there are directions in momentum space where quasiparticles can be excited with in- 
finitesimally small energy. For sufficiently small temperatures, these directions give 
the dominant contribution to the specific heat. Therefore, in the low-temperature 
approximation, we only integrate over angles in the vicinity of the nodes. We re- 
strict the angular integration by requiring the quasiparticle energy (with respect to 
the Fermi surface) to be at most of the order of the scale set by the temperature, 

/iosin0<7rr, (4.19) 

which, for small angles 9 and small temperatures implies 9 <n/(p. Therefore, the 
specific heat becomes (note the factor 2 since we obtain the same result for north 
and south pole) 



cv 



^ ' dx / d99(x~ + (p-9-) 



1 + coshx Jo 

5k^ pL^T 1 



(4.20) 



We see that instead of an exponential suppression we now get a power-law suppres- 
sion o= {T / AY of the specific heat compared to the non-superconducting result. In 
this sense, the specific heat measures how effectively the quasiparticle excitations 
are suppressed by the gap. Our result shows that the dimensionality of the zero- 
energy excitations in momentum space translates into the temperature dependence 
of the specific heat: in the normal phase, there is a two-dimensional Fermi surface 
that contributes at T = 0, while for an isotropic gap, this Fermi surface is, simply 
speaking, gone. The anisotropic gap ( 14.181 ) is an intermediate case, its suppression 
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lies between the normal and the completely gapped phase. One may thus expect 
that between the zero-dimensional point nodes and the fully gapped spectrum there 
is another intermediate case, namely one-dimensional line nodes, see problem [4.1| 

The low-temperature results for the specific heat are relevant for the physics of 
compact stars because the superconducting gap of either nucleonic superconduc- 
tivity/superfluidity or quark superconductivity may well be much larger than the 
temperature of the star. In particular, the specific heat is important in the context of 
the cooling of the star, for example through neutrino emissivity Ey. With Cy being 
the energy loss per unit time and volume through neutrino emission (for example 
through the processes (12.271 ) in nuclear matter or the processes ( I2.40l i in quark mat- 
ter), the relation between Cy, cy, and the change in temperature is 



(The minus sign is needed since a positive Ey is an energy loss, i.e., the temperature 
will decrease, dT /dt < 0.) Integrating this relation from a time to (with temperature 
Tito) = To) yields 



This shows that the ratio of the specific heat and the neutrino emissivity enters the 
cooling behavior of the star. Typically, for a given phase, the neutrino emissivity will 
exhibit a similar behavior as the specific heat. For instance, in a superconductor, the 
emissivity as well as the specific heat are exponentially suppressed in which case 
the subleading behavior becomes important. In a real compact star, however, there 
is most likely not just a single phase and the phase that dominates the behavior of 
the emissivity is not necessarily the one that dominates the specific heat. 

The neutrino emissivity is much more difficult to compute than the specific heat, 
and we devote a whole chapter to its discussion and to a detailed calculation for the 
case of quark matter, see chapter|5] 



4.2 Color-flavor locked (CFL) quark matter 

In our discussion of superconductivity and superfluidity in compact stars we first 
focus on a density regime where we can perform rigorous calculations from first 
principles. This is the regime of asymptotically large densities, where we deal with 
weakly coupled, deconfined quark matter^ The quarks are weakly coupled due to 
asymptotic freedom, which says that the coupling of QCD becomes weak for large 
exchanged momenta. For our purpose, the QCD coupling can be considered as a 
function of the quark chemical potential and becomes arbitrarily small for large 




(4.21) 




(4.22) 



^ We shall not go into details of neutron superfluidity and proton superconductivity. For a detailed 
review of these matters, see Ref. j^. A shorter discussion can be found for instance in Sec. 3.2 of 
Ref [5). 
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chemical potentials. In other words, quarks at infinite chemical potential are free. 
Because of this important property of QCD we may use perturbative methods at 
high densities. The high-density region of the QCD phase diagram shown in Fig. 
II. H is therefore maybe the best understood regime of QCD. The other regimes in 
that phase diagram are more complicated: we have seen that for nuclear matter one 
usually relies on phenomenological models; the high-temperature, small-density re- 
gion, where the QCD coupling also becomes small, has subtle nonperturbative ef- 
fects because of infrared degrees of freedom; first-principle QCD calculations via 
computer simulations (lattice QCD) are much more complicated than perturbative 
physics at high densities and are so far restricted to vanishing chemical potential. 

This possibility of understanding a region of the phase diagram rigorously from 
first principles is a good theoretical motivation to study ultra-dense quark matter 
However, for our astrophysical purposes we need to point out that these studies are 
valid at densities much larger than expected in compact stars. In a compact star, the 
quark chemical potential is at most of the order of < 500 Me V. The perturbative 
calculation of the energy gap 4, to be discussed in Sec. 14.31 can be estimated to be 
reasonable at chemical potentials of the order of /i > lO^MeV (!) Given this dif- 
ference of many orders of magnitude, extrapolation of perturbative results down to 
compact star densities may seem bold. However, the (rough) quantitative agreement 
of these extrapolations with different approaches, using phenomenological models, 
gives us some confidence that the ultra-high density calculation may be of relevance 
for astrophysical calculations. Furthermore, we shall also apply general arguments, 
based on symmetries, which we can expect to hold even at moderate density where 
the coupling is strong. In summary, the following discussion, strictly speaking only 
valid for extreme densities, is of theoretical interest and may also give us insight 
into compact star physics. 

At this point we may remember that we have already discussed the approach to 
compact star densities from the opposite side. In the Walecka model of Sec. l3.1l we 
have constructed the model such that we have reproduced properties of nuclear mat- 
ter at densities accessible in the laboratory. These densities are lower than the ones 
in compact stars. We had to extrapolate up to higher densities to obtain predictions 
of astrophysical relevance. Therefore, we learn that matter inside compact stars is 
quite hard to tackle; we have to approach it from different sides, and currently we do 
not have rigorous control over our approaches. This reflects the discussion begun in 
the introduction: it shows that the question "What is the matter composition inside 
a compact star?" is, due to our lack of understanding of dense, strongly-interacting 
matter, not only an application of QCD but also relevant to understand QCD. 

From this somewhat philosophical discussion now back to superconductivity in 
quark matter. Cooper's Theorem tells us that an attractive interaction, however small 
it may be, leads to the formation of a quark Cooper pair condensate. At asymptoti- 
cally high densities, this attractive interaction is provided by single-gluon exchange. 
We can formulate quark pairing in terms of representations of the color gauge group 
SU{3)c, 

5t/(3),: [3],®[3],= [3]^©[6]^ (4.23) 
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On the left-hand side we have two quarks in the fundamental representation, i.e., two 
complex three-vectors since the number of colors is three, Nc — 3. They interact in an 
antisymmetric (A) anti-triplet channel and a symmetric (S) sextet channel which are 
attractive and repulsive, respectively. The attractive channel thus provides an anti- 
triplet of diquarks which has (anti-)color charge. The attractiveness of this channel 
can be understood for instance from the existence of baryons. Namely, in a simple 
picture a baryon contains a diquark in the [3]^ representation. If it is made of, say, a 
red and a green quark it has color anti-blue. The baryon is then color-neutraUzed by 
combining this anti-blue diquark with a blue quark. 

An obvious property of a quark Cooper pair is that it is color-charged. Therefore, 
it breaks the color symmetry SU{3)c spontaneously. In analogy to electronic super- 
conductors, which break the electromagnetic L'(l)emi quark Cooper pairing is thus 
termed color superconductivity. For an extensive review of color superconductiv- 
ity see Ref. ||6l. The order parameter of color superconductivity is the expectation 
value of the quark-quark two-point function (yfyf). The color structure of this ob- 
ject has to be antisymmetric because the antisymmetric representation [3]^ is the 
attractive channel. The flavor structure is governed by the chiral symmetry group 
discussed in Sec. 13.3.11 For now, we may consider these sym- 
metries to be exact, since at the high densities we are working we may neglect 
all three quark masses compared to the chemical potential. Each of these global 
5t/(3)'s leads to the same representations as the color group, 

SUi3)f : [3]/® [3]/ = [3]^e [6]}, (4.24) 

with f ^ L,R. Since the overall wave function of the Cooper pair has to be anti- 
symmetric and since pairing in the antisymmetric spin-zero channel is preferred, we 
need to pair in the flavor [3]^ channel. In other words, the color-flavor structure of 
the Cooper pair is 

(W)e [3]>[3]^. (4.25) 

More specifically, with A, a,j5 < 3 being color indices and B,i,j < 3 being flavor 
indices, 

(V/,«C75rf>-e"'^^e,7fi</'l- (4-26) 

Here, we have added the Dirac structure with the charge-conjugation matrix C = 
i'f-'f' , leading to even-parity, spin-singlet pairing. The 3x3 matrix ^ now de- 
termines the specific color-flavor structure within the given antisymmetric repre- 
sentations. This shows that there are in principle many different possible color- 
superconducting phases. They are distinguished by different pairing patterns, i.e., by 
which quark pairs with which other quark. (At asymptotically large densities, where 
the flavor symmetries are exact, many pairing patterns are equivalent by symmetry 
and only a few physically distinct phases exist.) In particular, one may construct 
phases in which some of the quarks are paired while some others are not. 



As already mentioned in the introduction, we neglect the heavy quark flavors although in this sec- 
tion we consider asymptotically large densities. Since we are ultimately interested in extrapolating 
our results down to compact star densities, we only take u, d, and s quarks into account. 
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At high densities, the favored phase is the color-flavor locked (CFL) phase Q. 
We can characterize it by the following properties, 



(i) The CFL order parameter is given by 



(4.27) 



(ii) In the CFL phase, all quarks are paired with pairing pattern rd — gu, bu — rs, 
bd — gs, ru — gd — bs (where rd is a red down quark, gu a green up quark etc.), 
and there are 8 quasiparticles with gap A and 1 quasiparticle with gap 2A . 

(iii) The CFL phase has the following symmetry breaking pattern. 



SU{3)c X SU{3)r X SU{3)l x U{\)b ^ SU{3),+l+r x Z2 . 



(4.28) 



These three properties are in fact equivalent. Before discussing their physical 
implications, many of which can be read off from properties (ii) and (iii), let us 
show how the physical statement (ii) follows from the more abstract statement (i). 
To get a clear picture of the matrix structure of the order parameter, let us denote 



the bases of the color and flavor antitriplets [3]^ and [3]'^ by (y* 
(fe)(7 = —i£ijB- Then, we can write Eq. (14.27b as 
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(4.29) 



This 9x9 matrix is obviously symmetric, as required (the color-flavor structure is 
symmetric, giving overall antisymmetry through the antisymmetric Dirac structure). 
Its rows and columns are labelled with the nine quarks, ru, rd, rs, gu, gd, gs, bu, bd, 
bs. A nonzero entry indicates that the corresponding quarks pair We see that the 
matrix has a block structure with three 2x2 blocks and one 3x3 block. This leads 
to the pairing pattern given in point (ii). Note that this is a basis dependent statement. 
In particular, since the color symmetry is a gauge symmetry, (y/Cysy/) is a gauge 
variant object. The physically relevant statement, however, is the second part of 
point (ii) about the quasiparticle excitations. This statement is gauge invariant. The 
gap structure is given by the eigenvalues of the square of the above 9x9 matrix. 



(4.30) 



where Xr are the eigenvalues of 
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We shall prove the form of the quasiparticle excitations ( 14.30b in Sec. 14.31 Here we 
simply compute the eigenvalues A, . They are given by the solutions of 



det(A-L) =0. 



This can be rewritten as 

— exp [Trln(A — L)\ = exp 



Tr InA - y — 



(4.32) 



(4.33) 



We now have to compute L". First note that 

(J • I);;.'^ = -E^f^'^e.jA = -5,"8f + dfdf ^ = 5"H,j + . (4.34) 
This result can be used to compute 



=5L-A. 



(4.35) 



Consequently, all powers of L only have the matrix structures L and 1. Thus we 
make the ansatz 

L"^a„L + bn. (4.36) 
Multiplying both sides of this equation by L and using Eq. ( 14.35b yields 

fl„+i = 5a„ + bn , b„+\ = -4fl„ . (4.37) 

These recursion relations can be solved with the ansatz a„ = p". This yields the 
equation =5/7 — 4 which is solved by p\ =4 and p2 — I. Consequently, the 
general solution is the linear combination 

a„ = ap'l + l3p"2 = 4"a + l5. (4.38) 

From above we know ai = 1 and fl2 = 5 which yields a = — /3 = 1 /3. Hence 



Inserting this into Eq. (14.33b yields 

'L-1 



= exp < Tr 



L — 4 

ln(A-4) — ln(A-l] 



3 ' ' 3 
Now we use Trl = 9 and, from Eq. ( 14. 34b . TrL = 12. Thus we have 
= exp[ln(A-4)+81n(A-l)] = (A-4)(A-1)^ 



(4.39) 



(4.40) 



(4.41) 



68 



4 Superconductivity and superfluidity in a compact star 



Consequently, the eigenvalues of L are 1 (8-fold) and 4 (1-fold). Physically speak- 
ing, together with Eq. ( 14.301 ) this means that in the CFL phase 8 quasiparticle exci- 
tations have a gap A and 1 quasiparticle excitation has a gap 2A . This is the second 
part of point (ii). Of course, this discussion says nothing about the magnitude of A, 
which has to be computed from the QCD gap equation, see subsequent section. We 
leave it as an exercise to show that (iii) follows from (i), see problem l4~2l 

Points (ii) and (iii) reveal many important physical properties of the CFL state. 
Since these points are solely based on symmetry considerations, they are indepen- 
dent of the details of the interaction. Therefore, they can be expected to hold also at 
lower densities where perturbative QCD is not applicable. First, one may ask why 
CFL is the ground state and not any other order parameter given by a different matrix 
(^^. The simple answer is that the CFL order parameter is the only one in which all 
quarks participate in pairing, as we have seen. All other possible order parameters 
leave several excitations ungapped. Therefore, the CFL phase leads to the largest 
condensation energy and thus is the ground state at high densities (at lower densi- 
ties the situation is much more complicated). A more formal argument is that the 
CFL phase is the color superconductor with the largest residual symmetry group. It 
is thus a particularly symmetric state which also indicates that it is preferred over 
other color superconductors, although this is not a rigorous argument. 

From (iii) we read off the following properties of CFL, 

• CFL breaks chiral symmetry. We see that the CFL symmetry breaking pattern 
( I4.28l l is, regarding chiral symmetry, the same as in Eq. ( 13.65b . However, the 
mechanisms are different. The latter is caused by a chiral condensate of the form 
(WrVl), while the CFL condensate has the form (v/fli//^) (and the same with R 
L). At first sight, the CFL condensate thus preserves the full chiral symmetry, 
i.e., apparently one can still do separate L and R rotations without changing the 
ground state. However, the symmetry breaking occurs through the "locking" with 
color, i.e., in order to leave the order parameter invariant, a color rotation has to be 
undone by equal rotations in the left- and right-handed sectors. Although caused 
by different mechanisms, the two scenarios lead to similar physics. As for the 
usual chiral symmetry breaking, the CFL phase also has an octet of Goldstone 
modes. Since all fermions acquire an energy gap, these Goldstone modes become 
very important for the phenomenology of the CFL phase. Moreover, at lower 
densities, where the strange quark mass cannot be neglected, kaon condensation 
is expected in the CFL phase, not unlike its nuclear matter relative discussed in 
Sec. 13.31 The kaon-condensed CFL phase is usually called CFL-ZT" and will be 
discussed in the next subsection. Sec. 14.2.11 

• The color gauge group is completely broken. While spontaneous breaking of a 
global group leads to Goldstone bosons, spontaneous breaking of a gauge group 
leads to masses for the gauge bosons. Here, all gluons acquire a Meissner mass, 
just as the photon acquires a Meissner mass in an electronic superconductor A 
nonzero Meissner mass for a gauge boson is the field-theoretical way of saying 
that there is a Meissner effect, i.e., that the magnetic field can penetrate the super- 
conductor only up to a certain penetration depth. The inverse of this penetration 
depth corresponds to the Meissner mass. In the CFL phase, one linear combi- 



4.2 Color-flavor locked (CFL) quark matter 



69 



nation of a gluon and the photon remains massless. In other words, there is an 
unbroken U{1)q C SU{3)c+l+R7 generated by Q which is a Hnear combination 
of the original charge generator Q and the eighth gluon generator Tg (if you have 
done problem [4.2l vou can easily show this and determine the exact form of the 
Unear combination). This phenomenon is also called rotated electromagnetism. 
Since the admixture of the gluon to the new gauge boson is small, one may say 
that the CFL phase is a color superconductor but no electromagnetic supercon- 
ductor This is of relevance for compact stars since it implies that the CFL phase 
does not expel magnetic fields. 
• The CFL phase is a superfluid since it breaks the baryon number conservation 
group U{\)b- This is important since this is an exact symmetry, even at lower 
densities where finite quark masses become important. Therefore, there is always 
one exactly massless Goldstone mode in the CFL phase. 



4.2.1 Kaon condensation in CFL quark matter 

We have pointed out that chiral symmetry is not only broken in the hadronic phase, 
but also in CFL. This is by itself an interesting fact since it means that in QCD 
chiral symmetry is spontaneously broken at very low and very high densities. How 
about the region in between? This is unknown, but the possibility remains that chiral 
symmetry is, at small temperatures, broken for all densities. Since the symmetry 
breaking patterns of nuclear matter and CFL are identical (note that in a neutron 
superfluid also the t/(l )b is broken), this implies that possibly there is no real phase 
transition at moderate densities and small densities in the QCD phase diagram. In 
Fig. IL II this corresponds to the possibility that the "non-CFL" region is absent, at 
least at r = 0. 

Now let us use the chiral symmetry breaking of CFL for a concrete calculation. 
Since in the CFL phase all (quasi)fermions acquire energy gaps of at least A - 
whose magnitude we compute from first principles in Sec. 14.31 - the physics of the 
CFL phase at temperatures smaller than A is determined by the pseudo-Goldstone 
modes associated to chiral symmetry breaking (and the exact Goldstone mode from 
breaking of U{\)b which we do not discuss here). As we discuss below, A can be 
expected to be of the order of 10 MeV at densities present in compact stars. This 
is large enough to make fermionic excitations in a possible CFL phase in a star 
essentially irrelevant. Therefore, for astrophysical applications, the discussion of 
the physical properties of the Goldstone modes is crucial. 

In the context of kaon condensation in nuclear matter. Sec. 13.31 we have used an 
effective theory for the chiral field U and its interactions with nucleons. Also for 
the mesons in CFL we can write down such a theory. In this case, the chiral field is 
given by 

Z = <^l<^R, (4.42) 
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where (pi and (pR are the 3x3 matrix order parameters in the left- and right-handed 
sector. In our above discussion we have not distinguished between 0i and since 
in "pure" CFL we have = ^r = 1. For unitary 3x3 matrices 0£ and ^r, E is 
unitary, Z GU{3).lt thus contains 9 degrees of freedom, one of which one usually 
ignores since it corresponds to the ?]' which is heavy due to the explicitly broken 
U{1)a- Eight degrees of freedom remain, E G SU{3), and we can identify them as 
pions, kaons etc. just like in hadronic matter, see Eq. ( 13.67b . Despite the similarities, 
there is an important difference to hadronic matter: as one can see from the defini- 
tion of the chiral field ( 14.421 ). a meson in CFL is composed of two fermions and two 
fermion holes (each (p in Eq. ( 14.42b represents a diquark). For example, a neutral 
kaon should be viewed as an excitation K'^ ~ iisdu. Note that this "CFL kaon" has 
the same quantum numbers as the "usual kaon", composed of a particle and an an- 
tiparticle, ^ sd. Hence, if you want to construct a CFL kaon from a usual kaon 
you need to replace s ^ dil and d — s- us. This identification reflects the a«r/-triplet 
representation in Eq. ( 14.24b . As a consequence, the meson masses in CFL are or- 
dered inversely compared to the usual mesons. To see this, first note that the quark 
flavors {u,d,s), ordered with increasing mass, niu < < rris, have the anti-triplet 
counterpart {ds,us,ud). Here the masses (squared) have become ordered in the op- 
posite way, mdnis > mi,ms > mi,mii. Therefore, in nuclear matter (and ignoring finite 
density effects), OTj^o < m^o because nij^o °= m,, + m^i and nij^o °= + m^, whereas in 
CFL nij^ < nij^o because m^Q <^ niumd + m,ims and m^Q °^ mjms + niuins- We shall 
verify the form of the kaon mass in CFL below within the effective theory. 
The effective Lagrangian for mesons in CFL is given by 

^ ^ ^Tr[Voi:Voi:''' - vldiEdiE^] + ^detMTr[M-i(i: , (4.43) 

with 

yQE = dQZ + i[A,Z], (4-44) 

where M = diag{mii,m^i,ms) is the quark mass matrix. The matrix A enters the the- 
ory as the temporal component of a gauge field; it plays the role of an effective 
chemical potential for the field X. We shall see below how this translates into effec- 
tive chemical potentials for the neutral and charged kaons. 

The original works where this Lagrangian has been proposed are Refs. iHIH. 
There you can find detailed explanations about the structure of the effective La- 
grangian and its differences to the effective meson Lagrangian for hadronic matter 
(13.70b . Comparing with Eq. ( 13.701 ) we see that in CFL we do not have a term linear 
in the quark masses, rather only quadratic, ' detM °^ (and higher even powers 
which we have neglected). We also have different coefficients in front of the tem- 
poral and spatial part of the kinetic term, originating from the breaking of Lorentz 
invariance in a medium, V;i = 1/ V3. As for the hadronic phase, there are two con- 
stants f„ and a. This reminds us of the nature of effective theories like the ones given 
by Eqs. ( 13.701 ) and ( 14.43b : they are expected to give at least a qualitatively correct 
description even beyond the regime where the theory can be tested experimentally 
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or from first-principle calculations. The reason is that they are almost entirely deter- 
mined by symmetries. Only the coefficients have to be taken from the experiment or 
an underlying microscopic theory. The former is done in the effective theory of the 
hadronic phase. The latter, namely fixing the constants /;r and a from perturbative 
QCD, is done for the effective theory of CFL. In particular, one can expect that, if 
CFL is the ground state of dense quark matter at densities relevant for compact stars, 
the effective theory is a powerful tool to compute the phenomenology of a potential 
quark core of the star. 

Although terms of higher order in the fields and the mass matrix have already 
been neglected in Eq. ( 14.43b . the Lagrangian still looks complicate. The meson fields 
da appear in the exponent of I,, 



(4.45) 



with the Gell-Mann matrices A„, and thus they appear to all orders even in the given 
truncated theory. Let us first rewrite the Lagrangian by abbreviating Q = OaKi/ fn 
such that 

Z = e'2 = cos e + / sin Q . (4.46) 
Then, the various terms of the Lagrangian become 

i:i[d()EdQE^ = Tr[(5ocose)2 + (^osine)^] , (4.47a) 
i:r[dil.dil.^] = Tr[(Vcose)' + (Vsine)2] , (4.47b) 
Tr[[A,i:][A,i:]^] = 2Tr[A'-(Acose)2-(Asin2)2], (4.47c) 
ii:x[-dQl.[A,I.]' + [A,r\dQi:'^] = 2/Tr[(5ocose)[A,cose] 

+ (aosine)[A,sine]], (4.47d) 

and thus 

^ = ^Tr[A2-(Acose)2_(^sine)2+2fl(detM)M"'cose] 

+ ^Tr[(5ocose)2 + (c(osine)2 - [(Vcose)^ + (Vsine)^]] 

+i y Tr[(5o cos Q) [A, cos Q\ + (5o sin Q) [A , sin Q\\ . (4.48) 

Let us first interpret Q as a constant background, i.e., as the meson condensate, and 
neglect the fluctuations. This will allow us to compute the values of the various 
condensates at zero temperature. In general, all mesons may condense and the pa- 
rameters of the theory determine which of the condensates becomes nonzero. We 
recall from the above discussion that we expect the kaons, not the pions, to be the 
lightest mesons in CFL. Therefore, let us simplify Q by setting all fields except the 
kaon fields to zero. 
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7 / (/)4-/</>5\ 

e= £</>«!«= (l)6-i<h], (4.49) 

«=4 y ^4 + 06 + / 

with the dimensionless condensates (pa = Qa/ fn- With this ansatz we shall be able to 
construct a zero-temperature phase diagram that contains regions of no condensates, 
charged kaon condensates, neutral kaon condensates, and possibly coexistence of 
both. This is exactly the same ansatz as we have made in Sec. l3.3.3] for kaon con- 
densation in nuclear matter, see Eq. (13.80b . We can thus follow the steps below Eq. 
(llSOl l to obtain 

cose= 1-^(1 -COS0), (4.50) 



and 



where 



Q 

sine=|^sin0, (4.51) 



^^ = f^ + f^ + ^l + ^^. (4.52) 



Since we assume our condensates to be constant in time and space, only the first 
line of the Lagrangian ( 14.481 ) survives. The tree-level zero-temperature free energy 
is the negative of this Lagrangian and becomes 



U = ^Tr 
2 



2^^J^ (fl(detM)M-'G2 ^A^Q" 



1 COS0 ( „( A^*. Tk/f\Tk/^—\ \7./^2\ 

\ 

04 ^""^ > ' 02 



(1 — cos0)^ , , ^9, 9 sin"0,,^,9 



(4.53) 



where we have subtracted the "vacuum" contribution 

f/cFL = U{E = 1) = -/2fldetMTr[M-i] , (4.54) 

such that the state without kaon condensates, i.e., the pure CFL state has free energy 
f/ = 0. With the definitions of the matrices A and Q in Eqs. ( I4.44l i and ( |4.49l l, the 
notations 

0|+ = 02 + 02 , 02o = 02 + 02 , (4.55) 

and abbreviating A = diag(fli,fl2,a3), the various traces are 

Tr[a(detM)M->e2_^2g2j ^ - ^|+)0|+ + (m^o - ju2^)02o 

-2fl3(fli0|+ +fl202o), (4.56a) 
Tr[(Ae2)2] ^ (flj02^ +fl2</>lo)^ + fl30^ (4.56b) 
Tr[(A2)2] = 2fl3(fli0|+ +02010), (4.56c) 
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where we have defined the kaon chemical potentials and masses 

''2 2 2 

M..-^, (4.57a) 

= aniiiiins + m„) , = flm„(;7z,v +mj) . (4.57b) 

It will become clear below that these quantities really act as masses and chemical 
potentials for the kaons. For simplicity we have omitted the electric charge chemical 
potential in the Lagrangian which would have appeared in as an additional 
contribution. Inserting Eqs. (14.56b into Eq. ( 14.53b . we can write the free energy as 



J 71 



(1 — COS0) 



2(l-cos</,)^(^Mi|| + M2|i') . (4.58) 



Here and in the following we use, for notational convenience, the subscript 1 for 
and 2 for K^. To understand the expression for the free energy we consider the limit 
case of small condensates, 0, <C fn, i-e., (j>i = Qi/ fx ^ 1 for / = ^K^. Then we 
can expand U (0i , 02) up to fourth order in the condensates to obtain 



2 ' 2 
4^f + f02^ + f0r0|- (4.59) 



with 



{i-K-,K^), a.^-(^^. (4.60) 
6/^ 2 4/- 

We have thus reduced the effective theory to a two-component 0'* theory, cf. Eq. 
(IA.18I ) in appendix lA.il with effective coupling constants for the self-coupling of 
the kaons and an effective coupling constant a for the interaction between charged 
and neutral kaons. 

We may come back to the full free energy ( 14.58b to find the ground state of the 
system for arbitrary chemical potentials /ij, /i2. To this end, one has to minimize the 
free energy through the equations 

dU _ dU _^ 
(901 (902 

By construction, the free energy of the CFL state without kaon condensation, 0i = 
02 = 0, is given hy U ~Q. If one of the condensates vanishes, say 02 = 0, one of the 
equations ( 14. 61b is automatically fulfilled, and the other one becomes 
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0=1^ 
/I d^i 



= sin ^1 (;7Zj — /if cos ] 



(4.62) 



This has a nontrivial solution for m\ < /ip 



1 for OT[ > /ij 



cos^i = 



for ml < jx^ 



(4.63) 



and the free energy density becomes 

for mj > fif 

Ui<h=0) = l fl[,n\-lilf 



(4.64) 



for mf < /ij 



By symmetry, we find the same solution for 02 if we set ^\ = 0. Equating the free 
energies of the two phases 0i = 0, 02 7^ and 0i 7^ 0, 02 = one finds the condition 
for coexistence of two condensates. 



„2\2 



„2\2 



(4.65) 



This condition can also be obtained by assuming two nonvanishing condensates in 
Eqs. (I4.6II 1. As a result we obtain the phase diagram shown in Fig. 14.21 where we 
restrict ourselves io pL\,ii2> without loss of generality. 

What are the values of the kaon chemical potentials in the real world? In other 
words, where in the phase diagram of Fig. l4.2l does a compact star sit? Let us first see 
whether in a star we can expect the kaon chemical potentials to be larger than their 
mass, i.e., whether kaon condensation is possible. As discussed above, for quanti- 
tative predictions of the effective theory we rely on the results for the constants /jt 
and a at asymptotically large densities and their extrapolation down to densities in 
a compact star. This extrapolation yields 
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21-81n2 



18 



27r2 



(lOOMeV)^ 



:0.03, 



(4.66) 



where we used a quark chemical potential /i ~ 500 MeV and a fermionic energy 
gap A ~ 30MeV. Then, from Eqs. ( 14.57b we conclude that both kaon masses are of 
the order of ~ nif^o ~ (a '"light — 5MeV, where we used a quark mass 
for u and d quarks m\[o^i — 5MeV and a strange quark mass ~ 150MeV. The 
kaon chemical potential then is ~ /i^o ~ m^/(2/i) ~ 20 MeV. This suggests 
that the interior of the star sits outside the rectangle given by the dashed lines in 
Fig. |4.2[ i.e., if there is a color-flavor locked core in a compact star it is likely to be 
kaon-condensed CFL matter rather than "pure" CFL matter (we shall confirm this 
conclusion below for nonzero temperatures). Does this matter contain a charged or 
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Fig. 4.2 Zero-temperature phase diagram for kaon condensation in the /i^+-/i^o -plane. No con- 
densation occurs if the chemical potential is smaller than the meson mass. Coexistence of the two 
condensates is only possible along the (solid) line that separates the CFL-AT" from the CFL-K^ 
phase. This line is given by Eq. <4.65t and marks a first order phase transition. For large chemi- 
cal potentials, it approaches the line fi^^ = ^f^. The (dashed) lines separating either of the two 
meson-condensed phases from the pure CFL phase are second order phase transition lines. In the 
condensed phases, the condensate and free energy are given by Eqs. )4.63) and <4.64) , respectively. 




a neutral kaon condensate? Firstly, the slightly heavier d quark compared to the u 
quark makes the slightly heavier than the K'^. This asymmetry is taken into ac- 
count in Fig. 14. 21 Moreover, the electric charge of a potential condensate would 
require the presence of electrons to neutralize the system, which further disfavors 
the charged kaon condensate. We thus expect the CFL- A"" phase to be the most likely 
meson-condensed phase in CFL. 

As a second application of the effective theory for mesons in CFL let us compute 
an estimate of the critical temperature of (neutral) kaon condensation. This is im- 
portant to answer the question: if there is CFL matter in a compact star and if there 
is kaon condensation at zero temperature, at which temperature (i.e., at which point 
in the life of the star) does condensation set in? 

The full temperature-dependent theory defined by the effective Lagrangian is 
very complicated. We therefore expand the Lagrangian ( 14.481 1 up to fourth order in 
the matrix- valued field Q to obtain 



^ = :^Tr 



+^Tr[(*S)2-,.J(Ve)2 + 2i(*e)|A,01] , (4.67) 

where we have neglected terms of fourth order in Q which contain derivatives such 
as Tr[(2(3o2)^] etc, and where we have dropped the contribution constant in Q, 
which serves to normalize the free energy of the pure CFL state to zero, see remark 
below Eq. ( 14.53b . Next, one has to separate the condensate from the fluctuations, as 
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demonstrated in appendix lA.il The resulting Lagrangian has the same structure as 
given in the appendix for the 0'^ model, see Eq. (IA.17I ): a tree-level potential; terms 
of second order in the fluctuations which define the tree-level propagator; terms 
cubic in the fluctuations which correspond to interactions due to the presence of 
the condensate; and finally terms quartic in the fluctuations. Here we do not discuss 
the explicit structure of these terms in general, for details you may consult Ref 
ifTOll . We rather restrict ourselves again to the kaon degrees of freedom. As a further 
simplification, we set the charged kaon condensate to zero, 64 = O5 =0 (but keep 
the charged kaon fluctuations). This is motivated by the above discussion about the 
more favorable neutral kaon condensate. For the neutral kaon condensate we choose, 
without loss of generality, a direction in the degeneracy space of the condensate by 
setting ©7=0, and we denote 9 = = 06- The tree-level potential from Eq. (|4.59l l 
then simply becomes 



m 



2 ,,2 



u{9)= ' ^^' e' + !^e\ (4.68) 

The kaon sector of the inverse tree-level propagator is block diagonal, 
where 

f -K^ + m^^-li^ + aO^ -lijiikQ \ 

^01 = 9,9 9 ' (4.70a) 

with the abbreviation = k^ — v\k^. The verification of this form of the kaon 
tree-level propagator is left as an exercise, see problem 14.31 Analogously to the 
calculation in the appendix, we obtain the kaon dispersion relations. They are given 
by the poles of the propagator Do, which are the zeros of the determinant of the 
inverse propagator ' . The dispersion for the charged kaon is 



ef{k) = ^Jvlk^+m\ + a^^Tll\■ (4.71) 

We see that the condensate gives a contribution to the mass of the . For the 
neutral kaon we obtain 



{k) = \ El + Ai| T ^^lilEl + 5M^ , (4.72) 



where 



(4.73) 



4.2 Color-flavor locked (CFL) quark matter 



77 



Since kaon condensation breaks a global symmetry of the system, namely the U{1) 
associated to conservation of strangeness, we expect a Goldstone mode. (Notice the 
two-fold condensation process: due to the condensation of quark Cooper pairs, chi- 
ral symmetry is broken and pseudo-Goldstone bosons appear in the system; on top 
of that, these pseudo-Goldstone modes - here the neutral kaons - condense them- 
selves, breaking the global symmetry further and giving rise to another Goldstone 
mode.) This mode is expected to be gaplessH To check this expectation, we first 
compute the condensate from the tree-level potential (I4.68t . The nontrivial mini- 
mum of this potential is 

02 = M_Zi. (4.74) 

P2 

This impUes A^^^I^q + SM"^ = (S/il — '^2)" ^"^^ ^l=o + = ^^2 ^ "^2 which we 
can insert into the kaon dispersion ( 14.72b . The result is 

e2+(/t = 0)=0, (4.75) 

confirming the existence of a gapless mode. 

Following the calculation in the appendix, we can immediately write down the 
thermodynamic potential at finite temperature. 

In order to extract an estimate for the critical temperature, we expand the potential 
for large T, 



2 



2 Uvl ' J" ' A 

2k^ 4 , m\ + ml-2{^\ + ^) . 

45 vr 12v3 ' ^^-"^ 

The r'* term is easy to obtain and has also been discussed in appendix lA.il For the 
term we have neglected 8M in the neutral kaon dispersions ( 14.72b . Then they 
assume the same form as the ones for the charged kaons ( 14.71b and we can use 
the expansion for the pressure of a noninteracting Bose gas, see for instance the 
appendix of Ref . IfTTI . 

We have arrived at a potential with terms constant, quadratic, and quartic in d. 
Since we assume the existence of a condensate at T = 0, we have /Z2 > «J2, which 
we have argued to be realistic for densities in compact stars. Therefore, the quartic 
term is always positive, while the quadratic term starts from a negative value at 



^ Due to the weak interactions this mode acquires a small energy gap in the keV range which we 
neglect here. 
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r = and becomes positive for sufficiently large temperatures. Consequently, the 
nontrivial solution for the condensate ceases to exist when the coefficient in front 
of the quadratic term vanishes. This yields the condition for the critical temperature 
which we thus estimate to be 

With the definitions ( 14.601 1 we can express as a function of kaon chemical po- 
tentials and masses]! Before we interpret the result we point out a problem of the 
current approach. We have seen that at zero temperature, with 0{T =0) given by 
Eq. ( |4.74t . we have e/(A: = 0) = 0. At finite temperature we expect the condensate 
to melt, i.e., 9{T) < 9{T =0) for all T. In this case, however, the excitation of the 
Goldstone mode (which should remain gapless for all T < Tc due to the Goldstone 
theorem) becomes imaginary if written in the form ( 14.72b . This is clearly unphys- 
ical and due to the approximation we have made. The solution to this problem is 
to set up a more elaborate approximation scheme which evaluates the thermal kaon 
masses self-consistently. This is beyond the scope of these lectures, see Ref. ifTOl 
for such a treatment. 

It turns out that our estimate of the critical temperature coincides with the self- 
consistent calculation. We can therefore use Eq. (14.78b for a physical conclusion. 
With the definition of the effective coupling constants a and j52 in Eq. (14. 60b and 
the approximate numbers for the kaon chemical potentials and masses discussed 
below Eq. ( 14.66b we obtain Tc ~ 60MeV. This is of the order of or even larger than 
the critical temperature T^^^ for CFL itself. We do not aim to compute the critical 
temperature of CFL in these lectures. We simply give the (mean-field) result, 

j.CFL_2i/3.o.57zi, (4.79) 

where A is the zero-temperature gap. This relation differs by a prefactor of order one 
from the relation obtained from the usual Bardeen-Cooper-Schrieffer (BCS) theory, 
Tc ~ 0.574; see remark below Eq. ( 14.991 ) for the origin of this prefactor. For our 
present purpose it is sufficient to notice that the critical temperature in a supercon- 
ductor is typically of the same order as the zero-temperature gap. Since A is also of 
the order of tens of Me V, we may apparently conclude that the kaon condensate does 
not melt before the CFL phase itself melts. However, we need to remember that our 
effective theory is only valid for temperatures smaller than the gap A. Therefore, 
the estimated critical temperature for kaon condensation is close to or beyond the 
limit of validity of our effective description. Nevertheless, as a tentative conclusion 
we can say that as soon as quark matter is cold enough to be in the CFL state, we 
also expect it to be cold enough for kaon condensation, provided that the parame- 
ters are such that kaon condensation is present at zero temperature. In other words. 



* Notice that for a + lfii < the critical temperature fomially becomes imaginary, i.e., the con- 
densate apparently "refuses" to melt. This situation cannot occur for realistic parameters in our 
case but is an interesting theoretical possibility. See appendix C in Ref. 1101 and references therein 
for more information. 
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upon decreasing the temperature, one encounters the transition from unpaired quark 
matter to CFL-K^, not from unpaired quark matter to CFL and then to CFL-A"". The 
critical temperature we have found is larger than all temperatures we are interested 
in for compact star applications. Therefore, we have learned that the temperature 
inside a compact star is, for all times in the life of the star, sufficiently low for the 
CFL-K° phase. 



4.3 Color-superconducting gap from QCD 

Let us now go through a true QCD calculation from first principles. Our goal is to 
compute the gap A with perturbative methods. As explained above, this calculation 
can be expected to be strictly valid only at densities much larger than present in 
compact stars. 

In the theoretical treatment of superconductivity one introduces charge-conjugate 
fermions, which can be thought of as hole degrees of freedom. A hole is left in the 
Fermi sea if you remove a fermion. One might thus say that introducing fermion 
holes leads to an overcounting of the degrees of freedom because if the theory knows 
about all fermions it also knows about where a fermion is missing. And indeed, we 
have formally doubled the degrees of freedom. However, since in a superconductor 
quasiparticles are mixtures of fermions and fermion holes, this is a necessary ex- 
tension of the theory. The fermion spinors become spinors in the so-called Namhu- 
Gorkov space and the fermion propagator becomes a 2 x 2 matrix in this space. The 
Cooper pair condensate is taken into account in the off-diagonal elements of this 
propagator, i.e., it couples fermions with holes. The inverse tree-level propagator in 
Nambu-Gorkov space is 

with the inverse tree-level fermion and charge-conjugate fermion propagators 

[G±]-i = y^'K^ ± Ai7o = L I^o ± (m - efc)]7oA,±^ , (4.81) 



where 



A±* = ^(l±e7o7-k) (4.82) 

are projectors onto positive and negative energy states. Since our QCD calculation 
applies to asymptotically large densities, we can safely neglect all quark masses. 
See appendix IA.2l for a derivation of the tree-level fermion propagator and its rep- 
resentation in terms of energy projectors. From Eq. ( 14.81b we immediately get the 
tree-level propagators 
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The full inverse propagator S is obtained from a Dyson-Schwinger equation 

5-'=5o'+i:, (4.84) 

with the self-energy 



<P 

<P+ 



(4.85) 



In principle, E also has nonvanishing diagonal elements which we neglect here. The 
off-diagonal elements contain the gap function A (K), 

0+{K) ^A{K).^Y5 , = -A{K)^''y5 , (4.86) 

where ^ specifies the color-flavor structure of the color-superconducting phase; in 
the CFL phase ^ = J • I, see Eq. ( 14.291 ). From the Dyson-Schwinger equation ( I4.84l i 
we obtain the inverse propagator, which we formally invert to obtain the propagator, 

S=(T+^-), (4.87) 



F+ G 



with 



G± = ([Gj]-^-0^Gj0±) \ (4.88a) 
^ -G'^'P^G^ . (4.88b) 

The off-diagonal elements are termed anomalous propagators. They are typi- 
cal for all superconductors, see for example Ref. 1121 . From their structure ( I4.88bb 
we see that they describe the propagation of a charge-conjugate fermion that is 
converted into a fermion through the condensate (or vice versa). One can thus 
think of the condensate as a reservoir of fermions and holes, and the quasiparti- 
cles are not just single fermions but superpositions of states with fermion number 
-5-3-1135 

Inserting Eqs. ( 14.8 11 1, (14.83b . and ( 14.86b into Eq. ( I4.88ab . we compute the diagonal 
elements of the propagator (for simplicity we assume ./#^ ~ which is true in the 
CFL phase, but may not be true in other phases). 



. e=± 



ko±{n- ek) - 



-1 

A^L 



koTip- ek) 



A+^o } , (4.89) 



with L = as defined for the CFL phase in Eq. (14.31b . Now we write L in its 
spectral representation, 

L= £ K^r, (4.90) 



'•=1,2 



with Xr being the eigenvalues of L, Ai = 1, A2 = 4, and the projectors onto the 
corresponding eigenstates. 
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L-4 



L-1 



(4.91) 



3 ' 3 

Obviously, these projectors are complete, + 3^2 = 1; they are also orthogonal, 
3^2 ~ 0, as one can see with the help of Eq. ( |4.35t . We obtain 



^0 ± ~ e^) 



XrA' 



koTil^- ek) 



fco ± (m - ek) 



koTil^- ek) 



e,r "O \^k,rl 



with 



£1 



k.r 



{ek-^)^ + XrA^ 



(4.92) 



(4.93) 



The poles of the propagator are = i^^r' ^kr dispersion relations of 

the quasiparticles (e = +) and quasiantiparticles (e — — ). We have thus confirmed 
Eq. (14.30b . in particular we now understand why the eigenvalues of L appear in 
the excitation energies. Note that the structure of the dispersion relations is thus 
determined entirely by the color-flavor (and Dirac) structure of the order parameter, 
and thus ultimately by the symmetry breaking pattern. Only the calculation of the 
magnitude of A goes beyond simple symmetry considerations and depends on the 
form of the interaction between the fermions. 

Using the result K92h for G± and Eq. (I4.88bl ). one easily obtains the anomalous 
propagators. 



e.r "O 



(4.94) 



The gap equation is a self-consistent equation for the off-diagonal elements of the 
self-energy E. We shall not discuss the detailed derivation of the gap equation (see 
Sec. IV.A in Ref. ||6l for this derivation). The gap equation reads 



<I>HK)=g^-'£rTjF+iQ)fnDf,iK-Q) 



(4.95) 



where g is the QCD coupling constant, which will be our expansion parameter, 
where D'^y is the gluon propagator, and where r„ = A„/2 (a = 1, ... ,8) with the 
Gell-Mann matrices A^. In Figs. 14.31 and |4~4| we show the self-energy and the gap 
equation diagrammatically. 

The first step is to transform the matrix equation ( |4.95t into an equation for the 
scalar gap function A(K). To this end, we multiply both sides of the gap equation 
with 75.-#A^ from the right and take the trace on both sides. Furthermore, we ne- 
glect the antiparticle contribution e = — (and denote Ek ,- = £^^) and use the fact that 
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E = 





V 



Fig. 4.3 Diagrammatic representation of the one-loop self-energy in Nambu-Gorkov space. Curly 
lines are gluon propagators, double lines correspond to G+ (left-pointing arrow) and G (right- 
pointing arrow), single lines to Gq (left-pointing arrow) and Gq (right-pointing arrow), and the 
circles are the gap matrices <J>+ (cross-hatched) and <P^ (hatched). The vertices have the form 
87^ Ta with the QCD coupling g. 




Fig. 4.4 Diagrammatic representation of the gap equation which arises as follows. On the one 
hand, the one-loop self-energy is given by cutting a fermion line in the corresponding two-loop 
diagram of the effective action. In Nambu-Gorkov space, this yields the matrix of four diagrams 
shown in Fig. 14.31 On the other hand, the self-energy is given by Eq. <4.85t . Equating these two 
matrices leads to the gap equation in the off-diagonal elements. The algebraic form of the gap 
equation is given in Eq. )4.95) . It is a self-consistent equation for <P^ (equivalently, one may solve 
the equation for 'P^), and thus for the gap function A {K). 



the gluon propagator can be taken to be diagonal in color space, D'^^ — d"''Df^^ 
This yields 

^(^) - f^^LE4%Tr[7''75A-7^75A+]Tr[r/.^,^,7L.^]D,v(P) 



3 ¥■ 



2 AjQ) ^ 1 AjQ) 



3 ql - 



9.^ J 



<TT[fYsA^-Y'Y5A+]D^y{P), 



(4.96) 



where we abbreviated P = K Q, and where we have used the results for the color- 
flavor traces 



TT[Tj.J(.'3^iTa.J(] = 2Tr[r/.^^27;,.^] = ■ 



(4.97) 



It is left as an exercise to verify these traces. With the gluon propagator in Coulomb 
gauge, 
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DooiP) ^ D(:{P) , Doi{P)^0, Dij^{dij^piPj)D,{P), (4.98) 
where and D, are the longitudinal and transverse components, we have 

J y Q 

X [(l+q-k)A'(P)-2(l-p-qp-k)A(i^)] . (4.99) 

Again, it is left as an exercise to verify this result by performing the trace in Dirac 
space. The two terms on the right-hand side arising from i and e^ 2 are due to 
the two-gap structure of CFL. Let us for simplicity ignore this structure in the fol- 
lowing, i.e., we replace e,^ 2 by e^.i (for more details about the QCD gap equation 
for CFL, see Ref. 1141 ). This simplification does not change the main result which 
is the dependence of the gap on the QCD coupling g. The two-gap structure has 
a nontrivial effect for instance on the relation between the critical temperature and 
the zero-temperature gap, see Eq. ( |4.79t . In fact the 2'/^ in that equation is actually 
{Xf^^X^^^y^^ where the exponents 2/3 and 1/3 are the prefactors in front of the two 
fractions in Eq. ( 14.99b . 

For the sake of brevity, let us now skip a few steps in the calculation. One inserts 
the specific form of the longitudinal and transverse gluon propagators (in the so- 
called hard-dense loop approximation), performs the Matsubara sum and the angular 
integral. Details of all these steps can be found for instance in Ref. ||T3l . and one 
obtains 

Ate. A r' dcj ^ tanh ^ f In 1^ + In 1^ + i In -J^ ] . (4. 100) 

Here, the three terms in parentheses arise from static electric gluons, non-static mag- 
netic gluons, and (Landau-damped) soft magnetic gluons, respectively. The last of 
these terms is responsible for the leading behavior of the gap which will turn out to 
be different from the usual BCS behavior in electronic superconductors. The reason 
is the existence of a long-range interaction mediated by the magnetic gluons in QCD 
for which there is no analogue in the interaction of electrons in a metal. We have 
deiinedml=Nfg^li^/{67Z^) (Nf being the number of flavors), andM^ = {3n/4)m^^, 
and we have restricted the momentum integral to a small vicinity around the Fermi 
surface, q £ [jx — d,n + 5] (5 <^ jx), where we expect the gap function to be 
peaked. The three logarithms can be combined to obtain 



2 AjQ) ^ 1 A{Q) 



with 
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and where we have taken the zero-temperature limit tanh ^ ^ I. The logarithm can 
be approximated by 



1 , 

2^" 



bji 



@{k-q)\n — + 0{q-k)\n — 



Moreover, we define the new integration variable 

2b IX 



y = gln 



q-ll + eq 



and abbreviate 



x = g\n 



IbpL 



k- jl + Ek 



Ibix 
gin—, 



xo = gln — , 



(4.103) 



(4.104) 



(4.105) 



where A is the zero-temperature value of the gap at the Fermi surface, A = Aq=^. 
We have 



dy: 



-d(q-n), 



bile 



-y/g 



{q-ll + £qY 



(4.106) 



With the latter relation we approximate ln(^;U /e^) ~ y/g, \n{bii / ~ x/g to obtain 
A{x) =x r dyA(y)+ f dyyA{y). (4.107) 



We can rewrite this integral equation as a second-order differential equation. 



dA 



d^A 



dx 

This equation is solved by 



— = / dyA{y) ^ =-Aix) 



A{x) = Zicos(jc* — x) 



(4.108) 



(4.109) 



such that the value of the gap at the Fermi surface (which corresponds to jc = jc*) is 
A, and such that the first derivative of the gap at the Fermi surface vanishes, since the 
gap peaks at the Fermi surface. To compute the value of the gap at the Fermi surface, 
we insert the solution ( 14. 109b back into the gap equation ( 14.107b and consider the 
point X ~x*, 



A = A dyycos{x* —y) = A [cos(jc* —y) — ysin(x* —y)]y^ 
= A[l — cos(jc* —xq) +xosin(x* — xq)] . 



y=x' 



(4.110) 
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Since xq is of order g, we approximate cos(jt:* —xq) = cos jc* cosxq + sinx* sinxo ~ 
cosjc* +xosmx*, sin(x* — xq) = sinx* cosxq — cosx* sinxg ~ sinx* — xqcosx*, and 
thus 

4 ~/i(l-cosx*). (4.111) 

Hence, cosx* ~ and thus 

A=2biiexp(^-^^ . (4.112) 

This important result, first derived in Ref . ifTSll , shows that the color-superconducting 
gap is parametrically enhanced compared to the BCS gap in conventional supercon- 
ductors. In BCS theory there is a contact interaction instead of gluon exchange, and 
the resulting gap equation has the form 

Aocg^ [^d{q-^)-. (4.113) 

Jo Eq 

Here the gap does not depend on momentum and one obtains A °^ exp(— const/^^), 
i.e., the coupling appears quadratic in the denominator of the exponential. This is 
in contrast to the color-superconducting gap ( 14.1 12b where the coupling appears 
linear in the denominator of the exponential. As mentioned above, this is due to the 
long-range interaction from magnetic gluons. For more details and a more general 
solution of the QCD gap equation see Sec. IV in Ref. ||6l and references therein. 

The solution of the QCD gap equation is a weak-coupling result and thus only 
vaUd at very large chemical potentials where the QCD coupling is sufficiently small. 
It is nevertheless interesting to extrapolate this result to larger couplings. Of course 
one should keep in mind that this extrapolation has no theoretical justification. We 
show the gap as a function of the coupling in Fig. 14.51 We see the exponentially 
small gap at small coupling and observe a maximum of the gap at a coupling of 
about g ~ 4.2. For compact stars we make the following rough estimate. Accord- 
ing to the two-loop j8 -function (which should not be taken too seriously at these 
low densities), the coupling at /i = 400 MeV is g ~ 3.5. From Fig. 14.51 we then 
read off A ~ 80 MeV. However in our derivation of the result we have ignored a 
subleading effect which yields an additional prefactor ~ 0.2. Therefore, we can es- 
timate the color-superconducting gap for compact star densities to be of the order of 
A - lOMeV. 

This result suggests that the critical temperature of color superconductivity is also 
of the order of Tc ^ lOMeV, cf. Eq. ( 14.791 1. Remember that compact stars have tem- 
peratures well below that value (only in the very early stages of the life of the star, 
temperatures around lOMeV are reached). This suggests that color superconductors 
are viable candidates for the matter inside the star. More precisely, if there is decon- 
fined quark matter inside the star, it is very likely that it is in a color-superconducting 
state. 

We conclude this chapter about color superconductivity by noticing that, besides 
the strong-coupling nature, other interesting questions arise at lower densities. We 
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Fig. 4.5 Color-superconducting gap A over quark chemical potential /i as a function of the QCD 
coupling g. The curve shows the result from Eq. )4.ll2t with Nf = 2, predicting a weak-coupling 
behavior A/^loc exp(— const/g). The values of A/ fi for large coupling is a simple (and in principle 
unreliable) extrapolation of the weak-coupling result. 



have seen in Sec. 12.2.21 that in unpaired quark matter the Fermi momenta of up, 
down, and strange quarks split apart, see Fig. 12. II This is due to the nonzero strange 
quark mass and the conditions of neutraHty and weak equilibrium. In our discussion 
of superconductivity we have always assumed that the fermions that form Cooper 
pairs have identical Fermi momenta. This is true in the region of asymptotically 
large densities where the strange quark mass can be neglected. It is not true, how- 
ever, at lower densities. The different Fermi momenta rather impose a "stress" on 
the pairing Q It is a quantitative question whether the pairing gap is large enough to 
overcome this stress. Roughly speaking, if the gap is larger than the mismatch in 
Fermi momenta, the usual pairing is still possible. It is therefore conceivable that 
the CFL phase persists down to densities where the transition to hadronic matter 
takes place. If the gap is too small, however, or the mismatch too large. Cooper pair- 
ing in the conventional way is not possible anymore. There are several versions of 
unconventional pairing which may take over and constitute one or several phases 
between the CFL phase at high densities and hadronic matter. Some of them break 
rotational invariance and may lead to nodes of the gap in certain directions in mo- 
mentum space as discussed in the context of the specific heat in Sec. 14.11 Others 
even break translational invariance and exhibit crystalline structures. All of the un- 
conventional phases have in common that there is less, and less symmetric, pairing 
than in the CFL phase. There is less pairing because the CFL phase is the only color 
superconductor where all quarks are gapped in all directions in momentum space. 
There is less symmetric pairing because the CFL phase is the color superconductor 
with the largest residual symmetry group. In the phase diagram of Fig. II. H all color 

' Cooper pairing with mismatched Fermi momenta is an interesting general phenomenon and not 
only relevant for quark matter, but also in condensed matter physics and atomic physics. See for 
instance Ref. 1161 where mismatched pairing of fermionic atoms is investigated experimentally in 
an optical trap. 
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superconductors other than CFL are collectively denoted by non-CFL. From what 
we just said it is clear that this region of the phase diagram may either be completely 
absent or, if present, may itself contain several phase transition lines separating dif- 
ferent color superconductors. More details about stressed pairing in quark matter 
and unconventional color superconductors can be found in Ref. Q. 

In summary, we emphasize that not only the strong-coupling nature but also the 
less symmetric situation (due to the finite strange quark mass) complicates our un- 
derstanding of quark matter in compact stars. This supports the theme of these lec- 
tures that we need to compute properties of candidate phases and check them for 
their compatibility with astrophysical observations. In the following section we shall 
turn to one of these properties, namely the neutrino emissivity. 



Problems 

4.1. Specific heat for anisotropic superfluid 

Compute the low-temperature behavior of the specific heat for a gap function with 
fine nodes, i.e., instead of Eq. ( |4.18t , take A — > Zi|cos0| and apply analogous ap- 
proximations as for the case of point nodes. 

4.2. Symmetries of CFL 

Show that from the structure of the CFL order parameter given in Eq. ( I4.27l i it fol- 
lows that the CFL symmetry breaking pattern is given by Eq. (14.281 1. Hints: it is suf- 
ficient to treat the chiral group SU (3)l x SU{3)r as one single flavor group SU (3) f. 
A color-flavor transformation {U,V) G SU{3)c x SU{3)f with U = exp(i^,^T„), 
V = exp{i(l)iTa) acts on the order parameter as {U,V){J ■ I) = [UJ^U'^) {VIaV^). 
One then has to show that only SU (3 )c + / transformations leave the order parameter 
invariant. 

4.3. Kaon propagator 

Derive the inverse tree-level propagator for neutral and charged kaons given in Eqs. 
(14.69b and ( 14.70b from the Lagrangian (14.67b . 
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Chapter 5 

Neutrino emissivity and cooling of the star 



We have seen in Sec. |2] that measuring mass and radius of a compact star is not 
sufficient to deduce the matter composition inside the star; it is neither conclusive 
for a distinction between nuclear matter and quark matter nor between unpaired 
quark matter and color-superconducting quark matter. We now turn to an observable 
which is more sensitive to the microscopic properties of dense matter, namely the 
temperature of the star. More precisely, its cooling curve, i.e., the temperature as a 
function of the age of the star. Approximately one minute after the star is born, the 
temperature has cooled below 1 MeV and the star becomes transparent for neutrinos. 
Consequently, neutrinos (and antineutrinos) which are produced in the star can leave 
the system and carry away energy. Neutrino emission is thus the dominant cooling 
mechanism of a compact star in about the first million years of its life. After that, 
photon emission takes over. We shall not be concerned with this late regime here. 

A very detailed review about neutrino emissivity in nuclear matter is Ref. fH. 
If you are interested in a shorter review, also discussing quark matter, I recommend 
Ref. ID. Before turning to the microscopic calculation of the neutrino emissivity gy, 
let us discuss its importance for the cooling curves. First of all, as already discussed 
briefly in Sec. 14. II it is not only the emissivity which is important for the cooling. 
Once you know how much energy per time and volume is carried away, you need 
to know how this affects the temperature of the star. Hence you also need to know 
the specific heat. The specific heat cy is a thermodynamic quantity and thus much 
easier to compute than the neutrino emissivity. We have done so in Sec. 14. H and have 
seen that superconductivity has a huge effect on cy, namely, due to the energy gap, 
cv is exponentially suppressed at sufficiently small temperatures. We shall see that 
superconductivity has a similar effect on the neutrino emissivity. Besides Cy and cy, 
also the heat conductivity is important for the cooling behavior Most forms of dense 
matter are very good heat conductors, such that the star becomes isothermal. As a 
consequence, in a realistic star which may have layers of different phases of dense 
matter, cooling tends to be dominated by the phase with the highest emissivity and 
the phase with the highest specific heat. 
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5.1 Urea processes in nuelear matter 

In Fig. l5.ll we show some data and schematic comparison with calculations for the 
cooling curves. We see that there are different classes of processes which lead to 
significantly different cooling scenarios. The most efficient process is the so-called 
direct Urea proeess which leads to a very fast coolingQln nuclear matter, the direct 
Urea processes are 

p + e + Ve, p + e-^n + Ve- (5.1) 

We have discussed these processes in the context of j3 -equilibrium, where they serve 
to establish the relation jj.p + jXe = jj.„, assuming that neutrinos and antineutrinos es- 
cape from the star, /iy = 0. Here we are interested in the question how both processes 
contribute to the neutrino emissivity. Since it does not matter for the energy balance 
whether neutrinos or antineutrinos are emitted, both processes contribute - in chem- 
ical equilibrium - equally to the emissivity. For the neutron, proton, and electron, 
the dominant contribution in momentum space to the processes comes from the mo- 
menta close to the Fermi momentum. The neutrino momentum is of the order of the 
temperature T which can be neglected compared to the Fermi momenta. Therefore, 
momentum conservation for both processes in Eq. ( 15.11 ) reads 

kf ,„ = kf + kf ,e . (5.2) 

In other words, the Fermi momenta k^ „, k^p, and kpe must form a triangle. For 
this triangle to exist, the triangle inequality has to be fulfilled, 

kfM < kp.p + kp^e- (5.3) 

We know that in a neutral system we have kp^p — kp^e, and thus the triangle inequality 
becomes 

kF.n<2kF,p. (5.4) 

Consequently, with °^ k^p- (i = n,p), 

n„ <%np ^ — > ^ , (5.5) 

i.e., the proton fraction has to be larger than 11%. We have seen in Sec. 12. 11 that 
this is not the case for noninteracting nuclear matter. Interactions can change this, 
especially for very large densities. At lower densities, this means that the direct Urea 
process is strongly suppressed in nuclear matter 

This brings us to a second class of processes which are less efficient than the 
direct Urea process, but may be the most efficient ones to emit neutrinos when 
the direct Urea process is suppressed. Momentum conservation can be fulfilled by 



This process is as efficient in sucking energy out of the star as the Casino de Urea in Rio de 
Janeiro is in sucking money out of the pockets of the gamblers. Hence the name. 



5.1 Urea processes in nuclear matter 
36.0 
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Fig. 5.1 Effective surface temperature and luminosity L, vs. age of compact stars, taken from 
Ref. |3). Observed values are compared with different cooling scenarios, each represented by a 
band that reflects the large uncertainties in the microscopic calculations. 



adding a spectator neutron or proton. This is the so-called modified Urea proeess, 

N + n^N + p + e + Ve, N + p + e ^ N + n + Ve , N = n,p. (5.6) 

As can be seen from Fig. 15. II this process typically results in a much slower cooling. 
The cooling is thus very sensitive to the proton fraction of nuclear matter, especially 
around the threshold of 11%. In other words, this sensitivity provides a good check 
on the equation of state. Phenomenological models with equations of state which 
predict the proton fraction to be above this threshold can be excluded since the star 
would cool too fast. There are several other neutrino emissivity processes in nuclear 
matter which we shall not discuss here. Some of these processes happen only with 
superconducting protons and superfluid neutrons, and are due to constant formation 
of Cooper pairs. 
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5.2 Direct Urea proeess in quark matter 

The direct Urea processes in quark matter are 

d^u + e + Vf,, u + e^d + Ve, (5.7a) 
s^u + e + Ve, u + e^s + Ve- (5.7b) 

These processes obviously require the availabiHty of single quarks. If quarks are 
paired in Cooper pairs one first has to break a pair. This costs energy. Therefore, 
in a phase where all quarks are paired (gapped), such as the CFL phase, we can 
expect the direct Urea process to be strongly suppressed. As for the specific heat, 
we expect an exponential suppression at temperatures small compared to the gap 
(at larger temperatures, but still below the color-superconducting phase transition, 
thermal energy is available to break the pairs). Recall that the gap is of the order of 
10 MeV, and the temperature of the star is well below that. Therefore, the exponen- 
tial suppression exp(— Zi/r) forbids any sizable effect of the Urea process. Other 
processes coming from Goldstone modes dominate the neutrino emissivity in the 
CFL phase 13 . However, their contribution is much lower than that of the unsup- 
pressed direct Urea process. Therefore, if the core of a hybrid star is made of CFL 
quark matter, with outer layers of nuclear matter where any kind of Urea process is 
possible, the cooling properties are utterly dominated by these outer layers. 

We have briefly discussed that at lower densities the CFL phase may not be the 
ground state anymore. Any other color-superconducting phase will have ungapped 
modes0The simplest example is the so-called 2SC phase where all blue and strange 
quarks are ungapped while the others are gapped. There are more complicated can- 
didate phases with ungapped modes only in certain directions in momentum space. 
In any case, the neutrino emissivity of these phases will be dominated by these un- 
gapped modes, and thus will be comparable to the emissivity of unpaired quark 
matter. We are thus interested in the neutrino emissivity of unpaired quark matter. 
To be a bit more ambitious, let us discuss the emissivity in the 2SC phase. From this 
calculation we will obtain the result for the unpaired phase "for free" because of the 
unpaired modes in the 2SC phase. Furthermore, we learn something about comput- 
ing reaction rates in a superconductor which show some interesting features. And 
also we will see in an actual calculation why the emissivity of the gapped modes is 
exponentially suppressed. In other words, the goal of this section will be to under- 
stand 

• the role of the Cooper pair condensate and the energy gap on the Urea process 
(we shall estimate this qualitatively) 

• the result of the emissivity of unpaired (ultrarelativistic) quark matter (we shall 
compute this quantitatively). 

All we shall need from the 2SC phase is the propagator. From Eq. ( |4.92t we know 
that the general form of the propagator can be written as 



^ A possible exception is the color-spin locked phase which has Cooper pairs with total angular 
momentum one and which we do not discuss here. 
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G^-fA^l^^r'^p^, (5.8) 



where we used Eq. (14.831 1 and where we dropped the antiparticle contribution. Note 
that this form of the propagator assumes that all flavor chemical potentials are the 
same. For the neutrino emissivity we need to drop this assumption. The order param- 
eter in the 2SC phase is characterized by 0^ = 5a35^^ where is the color-flavor 
matrix from Eq. ( 14.261 ). For simplicity, we drop the strange quarks and consider only 
a two-flavor system of up and down quarks^ Then, the color-flavor structure of the 
gap matrix is 

^ = 72/3, (5.9) 

with the second Pauli matrix To in flavor space and /3 in color space, as defined 
above Eq. ( |4.29t . The color-flavor structure of the 2SC phase is much easier to deal 
with than the one of the CFL phase because color and flavor matrices factorize. 
Since = 1, we have ~ J^, whose eigenvalues are Ai = 1 (4-fold) and A2 = 
(2-fold). This is the formal way of saying that in the 2SC phase quarks of one color, 
say blue, remain ungapped. The projectors onto the corresponding eigenspaces in 
color-flavor space are 

^1=73^, ^2 = 1-4- (5-10) 

They are trivial in flavor space and project onto red and green quarks (which are 
gapped) and blue quarks (which are ungapped), respectively. 

In a neutral two-flavor system, up and down chemical potentials are different, 
namely /i„ + jif, = fi^i, where jXe turns out to be nonzero due to the neutrality con- 
straint. The generalization of the propagator (I5.8l l to this case can be written in terms 
of the flavor components (see problem fS.ll l 

with 



e,,,^y'(A-^)2 + A,42, 5^^^^, fl^^^, (5.12) 

and Xr, S^r as above (^^ now being only matrices in color space since the fla- 
vor components are written separately). This structure of the propagator and the 
resulting quasiparticle dispersion relations are interesting on their own, since they 
describe Cooper pairing with a mismatch in Fermi momenta, as discussed at the 



The weak interaction between u and s quarks is suppressed compared to tlie one between u and d 
quarks due to the Cabibbo angle. However, the finite strange quark mass may partially compensate 
this effect because it leads to a larger phase space for the Urea process. Here in these lectures we 
do not want to deal with these complications and thus simply consider a system of massless up and 
down quarks, and thus only the processes j5.7at . 
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Fig. 5.2 Neutrino self-energy Zy and W-boson polarization tensor 11 needed for the neutrino emis- 
sivity from the quark Urea process. 



end of Sec. 14.31 However, in the present context of neutrino emissivity, we are only 
interested in the qualitative features of the gapped modes. Thus we shall ignore this 
complicated structure of the propagator and temporarily set ji,, = /i^/. Only when we 
compute the emissivity from the unpaired modes we shall reinstate the difference in 
up and down chemical potentials. 

Next we need to set up the equation that determines the neutrino emissivity. One 
possible formalism is the finite temperature real-time formalism. We shall not ex- 
plain this formahsm but refer the reader for more details to the textbooks IS) and 161 . 
For our purpose it is enough to know that the real-time formalism can be used for 
nonequilibrium calculations. Therefore it is well suited for transport properties and 
neutrino emissivity. Since these properties are always close-to-equilibrium proper- 
ties, one often simply uses an equilibrium formalism, such as the imaginary-time 
formalism, and adds whatever is needed as a small out-of-equilibrium feature by 
hand. In the real-time formalism we can start from the kinetic equation 

i^TT[YoG<{Py)] = -Ti[G>{Py)L<{Py)-Z>{Py)GUPv)] , (5.13) 

where Gy and Gy are the so-called "greater" and "lesser" neutrino propagators, and 
Py is the neutrino four-momentum. The greater and lesser propagators are obtained 
from the retarded propagator in the same way as given in Eqs. ( 15.141 ) for the case of 
the W-boson polarization tensor. The trace in Eq. ( 15.131 1 is taken over Dirac space. 
The two terms on the right-hand side correspond to the two directions of both pro- 
cesses (I5.7ab . i.e., there is a neutrino gain term from d ^ u + e + Ve, u + e ^ d + Ve, 
and a neutrino loss term from u + e + Ve ^ d, d + Ve ^ u + e. Since neutrinos, once 
created, simply leave the system, only the gain terms, namely the directions given in 
Eq. ( |5.7at . contribute. The neutrino self-energies Ey '^ are given by the diagram in 
Fig. 15.21 The present formalism amounts to cutting this diagram. The figure shows 
that a cut through the internal m, d, and e lines produces two diagrams which rep- 
resent the Urea process. One part of the neutrino self-energies are the W-boson 
polarization tensors FI^'^, as shown diagrammatically in Fig. 15.21 They are defined 
through the imaginary part of the retarded polarization tensor ImHff, 



5.2 Direct Urea proeess in quark matter 



95 



n>iQ) = -2/[l+/B(<?o)]ImnK(e), (5.14a) 
n<{Q) = -2ifB{qo)lmnR{Q) , (5.14b) 

with the Bose distribution function fsix) = 1/ (e^/^ ^1)- We shall discuss the cal- 
culation of ImiTfl in detail below. The kinetic equation (15.13b becomes 

^/v f,Pv = \ -pr-T-^ i-Ad P«,Pv 

at 8 J (znj^pvPe 

■Xfripe - l^e)fB{pv + l^e- Pe)lmn^'' {Q) , (5.15) 

where fp (x) = 1 / (e'^/-^ + 1 ) is the Fermi distribution function, where, due to four- 
momentum conservation, 

Q^iPe-Pv-l^e,Pe-Pv), (5.16) 

and where 



L^'^(p„pv) = Tr [(70/7, - 7. p,) 7-(l - ^){yoPv - 7 • Pv) / (1 - 7^) 



(5.17) 



(In this section, Lorentz indices are denoted A, (J, . . . in order to avoid confusion 
with the subscript v which indicates neutrino quantities.) If you are interested in the 
details of the derivation of Eq. ( 15.15b or more details about the real-time formahsm, 
see Ref. 17] and references therein. In this reference the neutrino emissivity is com- 
puted in the same formalism; however, for anisotropic phases, which leads to more 
complicated calculations than we shall present here. The following is equally under- 
standable if you simply start with Eq. (15.15b whose features are physically plausible 
as we explain now. 

The left-hand side of Eq. ( 15.15b is the change of the neutrino occupation number 
in time. It is related to the emissivity by 

ev = -2^|^Pv/v(f,Pv), (5.18) 

where the factor 2 accounts for the contribution from antineutrinos. The neutrino 
emissivity is thus the change in energy per unit time and volume. Our task is to 
compute the right-hand side of Eq. ( 15. 15b and integrate over the neutrino momentum 
according to Eq. ( 15. 18b to obtain Cy. To understand the right-hand side of Eq. ( 15.15b 
first note that the vertex for the processes d ^ u + and e ^ v + W^ is given 
by 

r' = - ^ J. ^ y^(i-y^), (5.i9) 

2v2sin Ow 

with the Weinberg angle 9w- (For the process d i-^ u + there is an addi- 
tional factor Vuci from the Cabibbo-Kobayashi-Maskawa (CKM) matrix; however, 
V„ci ~ 1.) The W-boson propagators can be approximated by the inverse Vt^-boson 
mass squared since all momenta we are interested in are much smaller than this 
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mass ~ 80GeV. Thus, pulling out the constant factors of the vertices in the W- 
boson polarization tensor, we obtain the overall factor with the Fermi coupling 
constant 

Gf = — ^-A = 1.16637- 10"" MeV-2. (5.20) 

ml sin^ 0w 

The additional factors in the trace of Eq. ( 15.17b come from the electron and neutrino 
propagators. And finally, the distribution functions in Eq. ( 15.151 ) belong to the elec- 
tron and the VF -boson. Eventually, the Bose distribution of the W will drop out since 
the W-boson polarization tensor will turn out to be fg^, see below. This makes 
sense because the W does not appear in the initial or final state of the process we are 
interested in. 



5.2.1 W -boson polarization tensor 



Next we need to compute ImH^"^ for which we first compute 



^ K 



(5.21) 



where the trace is taken over Dirac, color, flavor, and Nambu-Gorkov space. We 
have defined P = K + Q,K and P will play the role of the u and d quark momentum, 
respectively. The weak vertices in Nambu-Gorkov space are 











-/(I + 7^)t^ 



(5.22) 



where t± = (ti ± 1^2) /2 are matrices in flavor space, constructed from the Pauli 
matrices Ti, T2. They take care of the fact that a u and a d quark interact at the 
vertices. Recall that, for notational convenience, we have pulled out the constants 
of the weak vertices already and absorbed them in the overall factor G^. With the 
quark propagator S from Eq. (14.87b . the trace over Nambu-Gorkov space yields 



K 

-Tr 



n^"{Q) = ^i:{Tr[/(l-75)T_G+(/^)7"(l-7^)T+G+(P) 

y^{\ + y^)T+G-{K)f{l + f)T-G- (P) 
/(I - y^)T-F-{K)f{l + f)T^F+{P) 
Yil + '^)T+F+{K)fil-'^)T+F-{P)]] 



-Tr 

-Tr 



(5.23) 



We see that there is a contribution from the anomalous propagators . The corre- 
sponding diagram in Fig. l5.3l is an instructive example for processes in a supercon- 
ductor which are only possible due to the Cooper pair condensate, see explanation 
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<ud> 




<ud> 

Fig. 5.3 Anomalous contribution to ttie IV-boson polarization tensor 11. The loop consists of two 
anomalous fermion propagators, according to Eq. j4.88bt . The lower propagator consists of a full 
fermion propagator (double line), the condensate (hatched circle), and a charge-conjugate free 
propagator (single line in opposite direction), and analogously for the upper one. Electric charge 
conservation at the weak vertices determines the flavor content of each line. As a consequence, one 
reads off that the condensate acts as a reservoir that can convert a u quark into a d quark hole and 
vice versa. 



in the caption of the figure. The anomalous contribution is thus only present for the 
gapped modes. We shall ignore it here for simplicity (it is smaller than the contri- 
bution from the normal propagators, but not negligibly small 18J). This leaves us 
with the first two traces in Eq. ( |5.23t which are the contribution of the normal prop- 
agators (they of course also contain the superconducting gap). It turns out that both 
traces are identical which we use without explicit proof. We thus continue simply 
with twice the first term, 

n^-(e) 2^^^Tr [/(I - r')T-ro^A"r"(i - 7')T+ro^.vAp 

K r.s 

X Tt T 9 T J y-'-'^^) 



7^2 p2 r,2 p2 

1 ^k.r "O 



where we have inserted the propagator ( 15.81 ). (Recall that we have set jj.,, = /ij tem- 
porarily to avoid complications; this is sufficient to discuss the effects of supercon- 
ductivity qualitatively, but eventually we shall reinstate the difference in ji,,, and 
to compute the result for unpaired quark matter In principle, for the 2SC phase we 
would have to use the propagators given in Eq. ( 15. lib ). The color-flavor traces are 

Tr[T_^iT+^i] = 2, (5.25a) 

Tr[T_^iT+^2] = 0, (5.25b) 

Tr[T_^2T+^i] = 0, (5.25c) 

Tt[t^^2t:+^2] = I- (5.25d) 

Recalling that projects onto the gapped red and green quarks and ^2 onto the 
ungapped blue quarks, this is easy to interpret: the weak interaction cannot change 
colors. Therefore, the quark loop in the polarization tensor - see right diagram in 
Fig. 15.21 - contains an up quark and a down quark of the same color. They are either 
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both gapped (then they are red or green, hence the result 2 in Eq. ( |5.25a| i). or they 
are both ungapped (then they are blue). There is no term involving one gapped and 
one ungapped quark. We thus get two contributions, 



^ K 



kQ-iji-k) po-{p-p) ko-iji-k) po-{n-p) 



pI - £2 



lA _ c2 



7 



"p,2 J 



where we abbreviated the Dirac trace 



,(5.26) 



(5.27) 



We notice that the second contribution in Eq. (15. 26b is obtained from the first upon 
setting 4=0. Thus, for notational convenience, let us simply continue with one 
color degree of freedom, say the first term without the factor 2, and denote e^^ = Ejtj . 
In the end it is then straightforward to get the full result. 

Next one has to perform the sum over the fermionic Matsubara frequencies. This 
technique is discussed in detail for a simple example in appendix lA. 2. II Here we 
need the more complicated result from problem [A.2l 



Po-£p 



^ h + ei (M - k)] [ep + e2(M - P)] /f (-61 £^)/F(e2£p) ^g) 



e[.e2 



qQ-eiek + e2ep 



fB{-eiek + e2ep) 



(Remember P = Q + K.) We comment on the physical meaning of the sum over 
the signs ei,e2 = ± below. To obtain the retarded polarization tensor, we need to 
replace qo^ qo^ ' i?. Then, the imaginary part is obtained by using the identity 



lim — ^ — = ^- =F in5{x) , 
where ^ denotes the principal value. This yields 



(5.29) 



imn^'^(e) 



e\ei 



{my 



«1 De2 



/F(-ei£fc)/F(e2£p) 



fB{-ei£k + e2£p, 
where we have defined the Bogoliubov coefficients 



Bi = \\\ 



pL-k 
£k 



(5.30) 



(5.31) 
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These coefficients appear in the theory of any kind of superconductor or superfluid, 
see for instance Ref. ||9l- Inserting the result ( |5.30l l back into Eq. ( 15. 15b yields 

xfF{pe-lie)fF{-eiek)fF{e2ep)5{qa-e^ek + e2ep). (5.32) 

As expected, the Bose distribution from Eq. ( 15.15b cancels with the denominator 
from Eq. ( 15.30b since on the one hand = Pe — Pv — fJ-e according to Eq. ( 15.16b and 
on the other hand qo — eiEii — e2£p according to the 5-function. 



5.2.2 Effect of superconductivity on Urea process 

Eq. (15.321) describes the change in the neutrino occupation number due to the pro- 
cess u + e ^ d + Vf,. (The other relevant process d ^ u + e + yields the same 
result and is taken into account by the factor 2 in Eq. ( 15.18b .) For this process one 
expects Fermi distributions of the form /«/,,( 1 —/</), the factors fe and /„ standing 
for the incoming fermions, and the factor I — fd standing for the outgoing fermion 
(for the neutrino, fv ~ 0). So what is the meaning of the sum over ej, e2? With 
/(-x) = 1 -/(x) it seems that all combinations /<./„/rf, /e/„(l -f^), fe{l-fu)fd, 
and /g(l — /„)(1 — fd) appear. In other words, also processes where both the up 
and down quark are created or annihilated apparently give a contribution. More pre- 
cisely, the quasiparticles, which are mixtures of up and down particles and holes, are 
allowed to appear on either side of the reaction. This is an interesting property of a 
superconductor or superfluid where particle number conservation is spontaneously 
broken and particles can be created from or deposited into the condensate. 

To see explicitly that in the unpaired phase only one of the four subprocesses sur- 
vives, let us define the new Bogoliubov coefficients and the new dispersion relations 

B^=^l+e^^, h^sgn{k-n)£k. (5.33) 

Then we use that for any function F we have 

^ rdkBiF{eek) = rdkBlF{-eh) . (5.34) 

g Jo eJQ 

This reformulation is useful to understand the mixing of particles and holes, which is 
manifest in the Bogoliubov coefficients. Had we taken the limit 4—7-0 with the orig- 
inal formulation in 5^, e^., we would have obtained the excitation energy = |A: — /i | 
which describes a hole for k < ji and a particle for k > jj.. The more conventional 
excitation e^ — k — pL which describes a particle for all k is only obtained as a limit 
using fi^, Ek (both formulations are of course physically equivalent). Now, since in 
the unpaired phase = 1 , B^T = 0, we see that only the subprocess with e\ =62 = 1 
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survives in the unpaired phase. The other three subprocesses are only possible in the 
superconducting phase. 

The general result in the superconducting phase has to be computed numerically. 
Here we proceed with a discussion of the behavior at temperatures much smaller 
than the gap, T <C 4 . The neutrino emissivity is obtained by integrating Eq. ( I5.32l i 
over the neutrino momentum according to Eq. ( 15.181 ). For the purpose of a simple 
estimate we may consider the expression 



(e'^2v^W + i^ ^ (5.35) 



where we have abbreviated 



(p=-, / = / dvv^ / dx / dy, (5.36) 

1 Jv^x.y Jo Jo Jo 

and introduced the new dimensionless variables 

^ ^ P-I^d ^ = — (5 37) 

The integration over the electron momentum has been rewritten as an integration 
over the t/-quark momentum. We shall discuss the calculation more explicitly for 
the case of unpaired quark matter below. Especially the angular integral, i.e., the 
phase space for the process, needs to be considered in detail. For now we are only 
interested in the suppression due to the gap. In the integrand of Eq. ( 15.351 ) one 
recovers the distribution functions for the electron, the M-quark, and the d-quaik. 

We may now perform the sum over ei and e2 and approximate gV^^^ ^ 1 and 

evT+ip^ ^ since (p ~^oa for small temperatures. Then the four terms, in the order 

(ei,e2) = (+,+),(-,-),(-,+),(+,-), become 

/ 1 1 



1 1 \ 



(5.38) 



The terms where ei, e2 assume different signs, i.e., the third and fourth term, yield 
contributions of the order of e^^f. They are thus even stronger suppressed than the 
first two terms which are identical and yield contributions proportional to e^f, 
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2(pe-'f I 2 -21-27 (pe^'P .(5.39) 



_ .2 



In the last step we have performed the remaining integral numerically which yields 
a numerical factor, unimportant for our present purpose. The main result is the ex- 
pected exponential suppression of the neutrino emissivity for case of gapped u and 
d quarks, Ev °= e^"^/^ for T ^ A. The full numerical solution, also taking into ac- 
count the temperature dependence of the gap A, shows that this approximation is 
valid up to temperatures of about T < Tc/3 where Tc is the critical temperature of 
superconductivity. 



5.2.3 Result for unpaired quark matter 

With the help of the new Bogoliubov coefficients (15.331 1 we can easily take the limit 
of unpaired quarks. For an explicit calculation of the emissivity for this case we 
need the following ingredients. First we need to perform the remaining traces in 
Dirac space and do the contraction over Lorentz indices. This is done in problem 
I5.2l with the result 

LxoiPe , Pv ) ^^"^ (k, P) = 64{pe Pe k) [pv Pv p) ■ (5.40) 

Next we observe that the result for the right-hand side of Eq. (15.32b would be zero 
without further corrections: we have to take into account so-called Fermi liquid 
corrections which are induced by the strong interaction. We have mentioned these 
corrections briefly in Sec. 12.31 see Eq. ( 12.78b . To lowest order in the strong coupling 
constant - which is related to the coupling g from Sec. l4.3l bv a.v = g^/ {4n) - we 
have 

We illustrate in Fig. l5.4l how these corrections open up the phase space for the direct 
Urea process. As a consequence, there is a fixed angle 0,,^/ between the u and d 
quarks, and the 5-function in Eq. (15.32b can be approximated by 

2 

5{pe-pv+k-p)c:^-^5{cOSe„d-COSeo), 00500 = I -K-^. (5.42) 

(We have reinstated the different chemical potentials /i„, /i^/.) We denote the angle 
between the neutrino and the d quark by Oy^ and approximate the factor 



{pe - Pe ■ k) {pv - Pv • p) 2^ePv - COS Oyd) ■ 



(5.43) 
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Fig. 5.4 Illustration of how Fermi liquid effects from the strong interaction open up the phase 
space for the direct Urea process in unpaired quark matter. Right-hand side (solid Fermi momenta): 
without Fermi liquid corrections, the Fermi momenta of the ultrarelativistic quarks are given by 
Pf,u = ^iu. PF.d = i^d- Start with the momentum of the up-quark, P[,. The circle with center at its 
tip indicates possible endpoints of the electron momentum p,,. Since pf,e = /i^ and + = /i,; 
(/3 -equilibrium), one cannot form a triangle with p„, p,, and the down-quark momentum p^;, unless 
one chooses the three vectors to be collinear. In this case, the triangle collapses to a line and the 
phase space for the Urea process vanishes. Note that the neutrino momentum pv ~ is negligibly 
small on the scale of the figure. Left-hand side (dashed Fermi momenta): the strong interaction 
changes the quark Fermi momenta to pf ,, ~ — k), pf j ~ /i,/(l — k) with K = 20^/(371). In 
other words, both Fermi momenta are reduced, but the down-quark Fenni momentum is reduced 
by a larger absolute amount. Since the electron Fermi momentum is not changed, a finite region 
in phase space opens up. The resulting triangle has a fixed angle between up- and down-quark 
momenta given by the values of the chemical potentials and K. 



This factor vanishes for the case of colHnear scattering. The effect renders it 
nonzero, hence this factor and in consequence the total neutrino emissivity is pro- 
portional to tty. Putting all this together and changing the integration variable from 
Pe to the d quark momentum p yields 




X 5 (cos 0„,/ - cos 0o) /f {Pe - Mf')/^ - M») [1 - ./f (p - l^d)] ■ (5.44) 
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Since we have taken only one color degree of freedom from Eq. (15.261 1, we have 
reinstated a factor A^^ = 3. Next we introduce the dimensionless variables x, y, v 
from Eq. ( 15. 37b . and with the definition (|5.18l l of the total neutrino emissivity we 
obtain 

T /■ c/X2„„ f dQn f dQir 



[InfJ {iKfJ {2nf 
X (1 — COS 0v(/)^(cos 0,,^/ — cos 00 ) 



ffF{v+x-y)fF{y)[l-Mx)]. (5.45) 

Here we have approximated the lower boundaries by —l-L„/ii/T — °°- With the inte- 
gral 

fjx j°°jyfp{y+x-y)fF(y)[l -fF{x)] = (5-46) 

and the (trivial) angular integral 

f dQ„^ f dQr, f dQi, , ,1 

JilifJil^Jil^^'^'''' ^-^^('^"^ - '^"^ - 16^ ' ^'-''^ 

we obtain the final result 

457 

£v ~ —a,Glne^„lidT^ ■ (5.48) 
63U 

This result has first been computed by Iwamoto in 1980 ifTOl . 



5.3 Cooling with quark direct Urea proeess 

From the result for the neutrino emissivity we can now get a simple cooling curve 
for unpaired quark matter Of course we shall ignore a lot of details of realistic stars. 
The result will simply show how a chunk of unpaired two-flavor quark matter cools 
via the direct Urea process. Nevertheless, the result is very illustrative and shows 
that the direct Urea process is indeed an efficient cooling mechanism. We use Eq. 
(14.221 1, which relates the temperature as a function of time to the emissivity and the 
specific heat. For the emissivity we use the result ( |5.48b . For the specific heat, recall 
the result (|4.13l l which is valid for two fermionic degrees of freedom, taking into 
account spin; we thus have to multiply this result by the number of colors and add 
up the contributions of u and d quarks, 



cy = (M,? + AiJ)r. 



(5.49) 




where we have defined 



Then, performing the integration in Eq. ( |4.22| l yields 

_315 M^HmJ 1 
9Ua,GlliePLuPLdT^'^' 

To get an estimate for this characteristic time scale, we assume pLj = 500 Me V, 
jU„ = 400 Me V, lie = 100 Me V, a, = 1, an initial temperature of To = 100 keV at an 
initial time fo = lOOyr, and use the value of the Fermi coupling ( 15.201 ) to obtain 

T~ 10"^yr~5min. (5.52) 

This is a very short time compared to the astrophysical time scales we are interested 
in. The function T{t) is plotted in Fig. 15.51 We see the rapid drop in temperature on 
a time scale of minutes down to a few keV. We thus recover the shape of the direct 
Urea cooling from Fig. 15. II For late times t S> fo, we have T{t) °^ f^'/"*. 



Problems 



5.1. 2SC propagator 

Show that for the case of different flavor chemical potentials the fermion propagator 
of the 2SC phase is given by Eqs. ( 15. lib and ( I5.12I I. 

5.2. Trace over Dirac space 

Show that 
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L^^aiPe , Pv ) ^^"^ (k, p) = 64{p, - p, • k) (pv - Pv • p) , (5.53) 
with Lifj{pe,pv) and 5^''^'^(k,p) defined in Eqs. ( I5.17l i and ( 15.271 ). respectively. 
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Chapter 6 
Discussion 



Let us summarize what we have learned about compact stars and dense matter, hav- 
ing in mind the two questions we have formulated in the preface. In addition, let us 
also list a few things which would in principle have fitted into these lectures topic- 
wise. We haven't discussed them in the main part because either I found them not 
suitable for a concise, and yet pedagogical, introduction or because they are simply 
beyond the scope of these lectures, such as some of the theoretical approaches to 
dense matter listed at the end of Sec. 16.21 And, well, some selection has to be made, 
so for some of the following points there is no good reason why they appear here 
and not in the main part. The volume of the main part is chosen such that it should 
conveniently fit into a one-semester course, maybe dropping one or two of the more 
specialized subsections. In Sec. 16.21 1 will give some selected references where inter- 
ested readers can find more information about the questions we haven't addressed 
in the main part. 



6.1 What we have discussed 

• Astrophysical observables and their relation to microscopic physics. The first 
thing you should have learned in these lectures is in which sense compact stars 
are laboratories for the understanding of dense matter The experiments we can 
do in this laboratory are less controlled as for example tabletop experiments in 
condensed matter physics. This means we cannot always measure the quantities 
we would like to know, or at least not to an accuracy we would need for our 
purposes. And it means that it is often impossible to switch on or off certain un- 
wanted effects at will, which would be desirable to extract an exact value for a 
given quantity. For instance we would ideally like to have a precise look into the 
interior of a compact star, but these kind of observations will always be indirect at 
best since the information we get from the interior is filtered through the surface 
and the atmosphere of the star. However, in spite of these restrictions (which, to 
some extent, have been and will be overcome through improvements in observa- 
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tional technology), we have seen that our observational data of compact stars can 
be closely linked to the properties of dense matter. Examples we have discussed 
in detail are the mass-radius relation which is related to the equation of state and 
the cooling curve which is related to the neutrino emissivity and the specific heat. 
• Theoretical approaches to dense matter We have emphasized at several points 
in these lectures that the density regime which is of interest for the physics of 
compact stars is very difficult to tackle. After all, this difficulty led us to consider 
compact stars not only as an application of QCD but also as an important means 
to understand QCD. The main reason for the theoretical difficulty is the strong- 
coupling nature of QCD. We have discussed attempts to approach the relevant 
density regime from two sides, coming from lower and higher density. 
First, we have discussed nuclear matter, for which we have solid knowledge at 
low densities. This knowledge is strongly built upon experimental data. In princi- 
ple, even a single nucleon is theoretically a very complicated object if considered 
from first principles. First-principle calculations, at least for sufficiently simple 
properties of nucleons, are possible in computer simulations, but effective theo- 
ries remain an important tool to describe nuclear matter, and they work well (by 
construction) at sufficiently low densities. One of the basic examples we have 
discussed is the Walecka model. However, finding the correct description of nu- 
clear matter at high density is a challenge, and astrophysical data can be used to 
rule out or confirm certain models. 

Second, we have discussed QCD from first principles in the context of decon- 
fined quark matter This approach is rigorous at asymptotically high density and 
therefore is interesting on its own right. We have discussed that it predicts the 
CFL state. Whether CFL persists down to densities relevant for compact stars is 
unknown. We have discussed that, to get a rough idea about the low-density re- 
gion, one may simply extrapolate the rigorous results. But this of course stretches 
the results beyond their range of validity. We have also introduced a more pow- 
erful approach to deduce intermediate-density properties from the high-density 
calculations. This approach relies on the symmetries of the CFL state. Building 
on these symmetries, one can construct an effective theory which can give us at 
least qualitative insight into the properties of CFL at lower densities, although 
this approach cannot tell us whether CFL is indeed the ground state of matter at 
densities present in the core of a compact star 



6.2 What we could have, but haven't, discussed 

• r-modes - bulk/shear viscosity. We have said little about the rotation of a compact 
star except for stating that it can rotate very fast, up to almost a thousand times 
per second. For the purpose of our lecture, however, the rotation frequency is 
a very interesting observable because it is sensitive to the microscopic physics. 
One of the reasons is as follows. 
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Certain non-radial oscillatory modes of a rotating star, in particular the so-called 
r-modesQ are generically unstable with respect to gravitational radiation. The 
reason can be understood in a rather simple argument. Consider the situation 
where the star rotates counterclockwise, seen from the polar view, and where an 
observer in the co-rotating frame sees non-radial oscillations which propagate 
clockwise. These modes lower the total angular momentum of the star, i.e., if 
the star's angular momentum is positive, the oscillations have negative angular 
momentum. Now assume that these propagating modes are seen from a distant 
observer to move counterclockwise, i.e., they are "dragged" by the star's rota- 
tion or, in other words, their angular velocity in the co-rotating frame is smaller 
in magnitude than the angular velocity of the star, seen from a distant observer 
The pulsations now couple to gravitational radiation. The emitted radiation has 
positive angular momentum since a distant observer sees the pulsations move 
counterclockwise. Consequently, the total angular momentum of the star must be 
lowered. This, however, means that the angular momentum of the oscillations, 
which is already negative, is increased in magnitude (becomes more negative). 
Therefore, the emission of gravitational radiation tends to increase the amplitude 
of the pulsation which in turn leads to a stronger gravitational radiation etc. This 
is the r-mode instability. Note that the rotation of the star is crucial for this argu- 
ment. In a non-rotating star, the effect of gravitational radiation is dissipative, i.e., 
the non-radial oscillations would be damped. For a nice pedagogical introduction 
into this general relativistic effect see Ref. |1 1. 

The energy loss from gravitational radiation due to the r-mode instability makes 
the star spin down drastically and quickly. Consequently, the observation of suf- 
ficiently high rotation frequencies implies that some mechanism must be at work 
to avoid the instability. The above argument for the instability is generic for all 
rotating perfect fluid stars. If there is dissipation, i.e., if the matter inside the star 
has a nonzero viscosity, the instability can be damped. Put differently, in order 
to rotate fast the star has to be viscous. This statement seems paradoxical at first 
sight but makes sense with the above explanation. In Fig. 16. II we show an ex- 
ample for critical frequencies of hybrid and quark stars, derived from viscosity 
calculations. 

In hydrodynamics, there are two kinds of viscosity, shear and bulk viscosityH 
Bulk viscosity describes dissipation for the case of volume expansion or com- 
pression while shear viscosity is relevant for shear forces. Both kinds of viscosi- 
ties are relevant for the damping of the r-mode instability, typically they act in 
different temperature regimes, bulk viscosity at rather large, shear viscosity at 
rather small temperatures. What is the microscopic physics behind the viscosity? 
Let us explain this for the case of the bulk viscosity. 

Imagine a chunk of nuclear or quark matter in thermal and chemical equilib- 
rium in a volume Vq. Now we compress and expand this volume periodically, 
y(f) = Vq + SVocoscot. In the astrophysical setting, these will be local vol- 

' Oscillatory modes of compact stars are classified according to their restoring force. In the case of 
r-modes, this is the Coriolis force. 

^ In the case of a superfluid, there are in fact several bulk viscosities. 
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T(K) 

Fig. 6.1 Critical rotation frequency (normalized to the Kepler frequency Qk, the upper limit for 
the rotation frequency beyond which the star would start shedding mass from its equator) as a 
function of temperature for hybrid and quark stars. If a star is put somewhere above the respective 
curves, the r-mode instability will set in and the star will spin down quickly. "APR" stands for 
a certain nuclear equation of state, "Bag" denotes unpaired quark matter in the bag model and 
the box labelled LMXB indicates the location of observed low-mass X-ray binaries. Within the 
given calculation they are located in a stable region for both hybrid and quark stars. For more 
explanations and details see Ref. (2) where this figure is taken from. 



ume oscillations where O) is typically of the order of the rotation frequency of 
the star Through the volume change the matter gets out of thermal and, pos- 
sibly, chemical equilibrium. The latter may happen if the matter is composed 
of different components whose chemical potentials react differently on a den- 
sity change. An example is unpaired quark matter with massless up and down 
quarks and massive strange quarks. The system now seeks to reequilibrate. For 
instance, if the compression has increased the down quark chemical potential 
compared to the strange quark chemical potential (in chemical equilibrium they 
are equal), the system reacts by producing strange quarks, for instance via the 
process u + d ^ u + s. If it does so on the same time scale as the external os- 
cillation, there can be sizable dissipation (think of compressing a spring which 
changes its spring constant during the process; you will not get back the work 
you have put in). Consequently, the calculation of the bulk viscosity requires the 
calculation of the rate of processes such as u + d ^ u + s which indeed turns 
out to be the dominant process for the bulk viscosity in unpaired quark matter. 
Other processes which contribute are leptonic processes, such as the direct Urea 
process we have discussed in the context of neutrino emissivity in Sec. 15.21 It is 
important to note that again the weak processes are the relevant ones. In princi- 
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pie, also strong processes contribute to the bulk viscosity since they reequilibrate 
the system thermally. However, they do so on time scales much smaller than 
the external oscillation. Therefore, the system reequilibrates basically instanta- 
neously during the compression process and no energy is dissipated. These argu- 
ments also show that the bulk viscosity is a function of the (external) frequency. 
Maximum bulk viscosity is obtained when the rate of the respective microscopic 
process (which is a function of temperature) is closest to this frequency. Hence, 
it may well be that for a certain temperature regime a superconducting state has 
larger bulk viscosity than a non-superconducting state. This may sound coun- 
terintuitive but note that the (partial) suppression of the rate of the microscopic 
process by exp(— zi /T) may actually help the viscosity if it brings the rate closer 
to the external frequency. See for instance Sec. VII in Ref. Q for a brief review 
about viscosity in quark matter phases, and Refs. ||4]|5] for examples of detailed 
microscopic calculations of bulk viscosity in quark matter. 
• Magnetic fields. We have mentioned in the introduction that compact stars can 
have huge magnetic fields, the highest magnetic fields measured for the surface 
of a star (then called magnetar) are about lO'^ G. The first question one might 
ask is what the origin of these magnetic fields is. The conventional explanation is 
that they are inherited from the star's progenitor, a giant star that has exploded in 
a supernova. While the magnetic flux is conserved in this process, the magnetic 
field is greatly enhanced because the magnetic field lines are confined in a much 
smaller region after the explosion. 

Other questions regarding the magnetic field concern their interplay with dense 
matter. We have learned that nuclear matter can contain superconducting protons. 
Protons form a type-II superconductor where the magnetic field is confined into 
flux tubes. Since at the same time the rotating neutron superfluid forms vortices, 
a complicated picture emerges, where arrays of flux tubes and vortices intertwine 
each other. Their dynamics is complicated and relevant for instance for the ob- 
served precession times of the star, see for instance Ref. IS). This issue is also 
related to pulsar glitches, see below. 

In the main part we have only touched the interplay of color superconductors 
with magnetic fields. We have stated without calculation that the CFL phase is 
not an electromagnetic superconductor, i.e., a magnetic field can penetrate CFL 
matter. More precisely. Cooper pairs in CFL are neutral with respect to a certain 
mixture of the photon and one of the gluons. Because of the smallness of the 
electromagnetic coupling compared to the strong coupling, the gluon admixture 
is small and the new gauge boson is called "rotated photon". There are color su- 
perconductors which do expel magnetic fields, for instance the color-spin-locked 
(CSL) phase. In this case. Cooper pairs are formed of quarks with the same flavor, 
and a Cooper pair carries total spin one (instead of zero in the CFL phase). The 
CSL phase is an electromagnetic superconductor. It is of type I, i.e., expels mag- 
netic fields completely. For a short review about spin-one color superconductors 
in compact stars and their effect on magnetic fields see Ref. Q. 
Magnetic fields also play a role in the cooling of the star since they have an effect 
on the heat transport, resulting in an anisotropic surface temperature, see Ref. ||8l 
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and references therein. An extensive review about magnetic fields in compact 
stars is Ref. 

• Crust of the star. The crust of the star is a very important ingredient for the un- 
derstanding of observations. In the conventional picture of a neutron star there is 
an outer crust with an ion lattice, and an inner crust with a neutron (super)fluid 
immersed in this lattice. This crust typically has a thickness of about 1 km. A lot 
about the crust can be found in Refs. ifTolfTTl . In our discussion of neutron stars 
vs. hybrid stars vs. quark stars it is important that the crust provides a crucial 
distinction between an ordinary neutron star (or a hybrid star) and a quark star. 
How does the crust of a quark star look? Several scenarios have been suggested. 
First suppose that the surface of a quark star exhibits an abrupt transition from 
strange quark matter to the vacuum. This is possible under the assumption of the 
strange quark matter hypothesis we discussed in Sec. 12.2.11 because, if the hy- 
pothesis is true, strange quark matter is stable at zero pressure. "Abrupt" means 
that the density drops to zero on a length scale of about 1 fm, given by the typi- 
cal length scale of the strong interaction. Now recall that (unpaired) three-flavor 
quark matter contains electrons. They interact with quark matter through the elec- 
tromagnetic interaction, therefore their surface will be smeared (several hundred 
fm) compared to the sharp surface of the quark matter. As a consequence, an 
outward-pointing electric field develops (i.e., at the surface positively charged 
test particles are accelerated away from the center of the star). This electric field 
can support a thin layer of positively charged ions, separated from the quark mat- 
ter by a layer of electrons. Hence a "normal" crust for a quark star is conceivable, 
consisting of an ion lattice. In contrast to the crust of a neutron star, such a crust 
of a quark star would be very thin, at most of the order of 100 m. See Ref. lfT2l 
for more details about this picture of the surface of a quark star (and for other 
properties of quark stars). This picture may be challenged by the possibility of a 
mixed phase at the surface of the star. Here, mixed phase refers to a crystalline 
structure of strangelets immersed in a sea of electrons. In this case, there would 
be no electric field and thus no possibility for a "normal" crust. The quark matter 
would rather have its own crystalline crust. Estimates in Ref. ifTSl show that it is 
unlikely that such a mixed phase is formed once surface tension is taken into ac- 
count. In any case, a rigid crust, if at all present, will be much thinner in a quark 
star than in a neutron star or a hybrid star 

This difference is relevant in the context of "magnetar seismology". Quasi- 
periodic oscillations observed in the aftermath of X-ray bursts from magnetars 
can be related to typical oscillation frequencies of the crust. In other words, 
"star quakes" have significantly different properties depending on whether one 
assumes the star to be a neutron star or a quark star In fact, the ordinary crust ex- 
plains the data quite well while the crust of a quark star seems to be incompatible 
with the observed phenomenology [141. 

• Pulsar glitches. Pulsar glitches are an interesting phenomenon related to the ro- 
tation frequency, the crust (more precisely the crystalline structure of the crust), 
and superfluidity. For spinning-down pulsars one observes sudden spin-ups, i.e., 
in the overall trend of a decreasing rotation frequency, the frequency increases 
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in irregular intervals significantly on a very short time scale. This is convention- 
ally explained through superfluid vortices in the neutron superfluid that pin at the 
lattice sites of the inner crust ifTSi . 

To understand this statement and the consequences for glitches, we recall the 
following property of superfluids. A superfluid, be it superfluid helium, super- 
fluid neutron matter, or any other superfluid, is irrotational in the sense that the 
superfluid velocity has vanishing curl. Therefore, if the superfluid is rotated it 
develops regions where the order parameter vanishes, i.e., where it becomes a 
normal fluid0 The angular momentum is then "stored" in these regions which 
are called vortices. An array of vortices, which are "strings" in the direction of 
the angular momentum, is formed with the total angular momentum of the su- 
perfluid being proportional to the density of vortices (because each vortex carries 
one quantum of circulation). Consequently, if the rotation frequency decreases, 
the array of vortices becomes sparser, i.e., the vortices move apart. 
The next ingredient in the glitch mechanism is the pinning of the vortices at the 
lattice of nuclei in the inner crust. Generally speaking, there is an effective in- 
teraction between the vortices and the nuclei, resulting in a certain path of the 
vortex string through the lattice which minimizes the free energy of the system. 
You may think of this preferred configuration as follows. Superfluidity, i.e., neu- 
tron Cooper pairing, lowers the free energy of the system. Therefore, the system 
may want to put the vortices, where there is no Cooper pairing, through the lattice 
sites because they are not superfluid anyway. Otherwise, i.e., by putting them be- 
tween the lattice sites, one loses pairing energy. The actual details of the pinning 
mechanism are complicated and, depending on the density, the preferred path of 
the vortices may in fact be between the sites, in contrast to the above intuitive 
argument. However, this does not matter for our argument for the mechanism of 
glitches: 

In a rotating neutron star, the neutron vortices pin at the lattice of the inner crust. 
Now the star spins down. On the one hand, the vortices "want" to move apart. 
On the other hand, there is an effective pinning force which keeps them at there 
sites. Hence, for a while they will not move which implies that the superfluid 
(the vortex array) is spinning faster than the rest of the star At some point, when 
the tension is sufficiently large, the vortices will un-pin, move apart and thus 
release their angular momentum which spins up in particular the surface of the 
star whose rotation is observed. Then, they re-pin and the process starts again. 
An alternative scenario, where nuclear matter is replaced by quark matter has 
been suggested tl6|. In our discussion of the CFL phase we have seen that quark 
matter can be a superfluid. This means that one of the conditions for the mecha- 
nism of vortex pinning is fulfilled. The second condition, a sufficiently rigid lat- 
tice, may be provided by one of the unconventional color-superconducting phases 



It is instractive to view this phenomenon in analogy to a type-II superconductor. There, a inag- 
netic field (if sufficiently large but not too large) penetrates the superconductor through flux tubes. 
It partially destroys superconductivity, i.e., in the center of the flux tubes the order parameter is 
zero. Hence the analogy is superfluid - superconductor; angular momentum - magnetic field; vor- 
tices - flux tubes. 
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which are possible in the case of mismatched Fermi momenta, see discussion at 
the end of in Sec. 14. 31 Some of these phases indeed exhibit a crystalUne structure. 
Such a quark crystal is of very different nature than the ion lattice because it is 
the energy gap from superconductivity which varies periodically in space, giving 
rise to crystals characterized by surfaces where the gap vanishes. It remains to be 
seen in the future which of these scenarios passes all observational constraints 
and can explain the pulsar glitches or if there is a yet unknown mechanism for 
these curious spin-ups. 
• Other theoretical approaches to dense matter What are the alternatives to under- 
stand QCD at large, but not asymptotically large, densities? Lattice QCD, i.e., 
solving QCD by brute force on a computer, is by now a powerful tool for strong- 
coupling phenomena at zero chemical potential. However, at nonzero chemical 
potential, one encounters the so-called sign problem which renders lattice calcu- 
lations unfeasible. Progress has been made to extend lattice calculations to small 
chemical potentials, more precisely to small values of fi/T. But calculations at 
large ji and small T, as needed for compact stars, are currently not within reach. 
See Ref. ifTTl for a non-technical recent overview article about lattice QCD, in 
particular its contributions to the QCD phase diagram and about the sign prob- 
lem; you may also try Ref. ifTsl . 

Because of the problems of lattice calculations at finite chemical potential one 
has to rely on model calculations or on extrapolations similar to the ones dis- 
cussed in these lectures. One model for quark matter we have not discussed is 
the Nambu-Jona-Lasinio (NJL) model. This model does not contain gluons and 
describes the interaction between quarks by an effective pointlike interaction. 
It has been used to compute the QCD phase diagram at intermediate densities. 
Since the result depends strongly on the parameters of the model, it should be 
taken as an indicator for how the phase diagram might look, not as an accurate 
prediction. Due to its simplicity it is widely used and can indeed give some inter- 
esting results which serve as a guideline for the understanding of QCD. For an 
extensive review about the NJL model in dense quark matter see Ref. |fT9l ; for 
an application of the NJL model in the context of compact stars, see for instance 
Ref. Il20l . 

Finally, we point out that arguments for large numbers of colors A^^ may be ap- 
plied to gain some insight to QCD where A^^ = 3. In particular, it has been ar- 
gued that at Nc ^ °° an interesting novel phase, termed quarkyonic matter, pop- 
ulates the T-jj. phase diagram 1211. The (yet unsolved) problem is to find out 
whether this phase, or some modification of it, survives for Nc = 3. More gener- 
ally speaking, the large Nc approach is another approach where calculations can 
be performed in a regime where everything is under rigorous control. From these 
rigorous results one then tries to get closer to the regime one is interested in. In 
this sense, this approach is not unlike the perturbative approach. In view of the 
possible, but not at all obvious, relevance of large-A'^ physics to Nc = 3 physics, 
one can also apply the duality of certain string theories to field theories similar to 
QCD, based on the so-called AdS/CFT correspondence. For pedagogical reviews 
see Refs. 1221 123| . This somewhat speculative but popular approach to QCD has 
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recently been pursued especially for large-r, small-^u physics, but is, in certain 
variants, also suited for the physics at finite chemical potential. 
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Appendix A 

Basics of quantum field theory at finite 
temperature and chemical potential 



Many of the discussions in the main part of these lectures rely on field-theoretical 
methods, in particular on quantum field theory at finite temperature and chemical 
potential. One purpose of the following basic discussion is therefore to explain how 
a chemical potential is introduced in quantum field theory. We shall also discuss how 
finite temperature enters the formalism, although for most quantities we discuss in 
these lecture notes we consider the zero-temperature limit, which is a good approx- 
imation for our purposes. For instance in the discussion of the Walecka model. Sec. 
13.11 we give the finite-temperature expressions, based on appendix I A. 21 before we 
set r = in the physical discussion. In other parts, we do keep T 7^ in our results, 
for instance when we are interested in the cooling behavior of dense matter, see 
chapter|5] 

We shall start with the Lagrangian for a complex bosonic field and derive the 
partition function in the path integral formalism, taking into account Bose-Einstein 
condensation. This part is particularly useful for our treatment of kaon condensation 
in CFL quark matter, see Sec. 14.2.11 We shall in particular see how bosonic Matsub- 
ara frequencies are introduced and how the summation over these is performed with 
the help of contour integration in the complex frequency plane. In the second part 
of this appendix we shall then discuss the analogous derivation for fermions. 



A.l Bosonic field 

We start from the Lagrangian 

^O^d^(p*d^(p-m^\(p\--X\(p\\ (A.I) 

with a complex scalar field (p with mass m and coupling constant A . We shall first 
show how a chemical potential jj, is introduced. This will lead to a new Lagrangian 
Jff, wherefore we have denoted the Lagrangian without chemical potential by J^'q. 
The chemical potential jj. must be associated with a conserved charge. We thus need 
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A Basics of quantum field theory at finite temperature and chemical potential 



to identify the conserved current. From Noether's theorem we know that the con- 
served current is related to the symmetry of the Lagrangian. We see that J^q is 
invariant under U{1) rotations of the field. 



(p^e-'"(p. (A.2) 



This yields the Noether current 



d{d^(p) 5a d{d^(p*) 5a -r -r -r j, 

with dj^jf^ = 0, and the conserved charge (density) is 

/ = /((p*aV-(pA*)- (A.4) 

In the following we want to see how the chemical potential associated to enters 
the Lagrangian. The partition function for a scalar field is 



Z = Tr( 



2i% I Si(p exp 

periodic 



X 



(A.5) 



This equation should remind you that the partition function can be written in the 
operator formalism in terms of the Hamiltonian H and the charge operator .^V, or, 
as we shall use here, in terms of a functional integral over (p and the conjugate 
momentum n, with the Hamiltonian and the charge density o/K = j'^. We have 
abbreviated the space-time integration by 







dxjcPx, (A.6) 



where the integration over "imaginary time" t = it goes from to the inverse tem- 
perature j3 = \/T. In the following, the four- vector in position space is denoted 
by 

X=(f,x) = (-/T,x). (A.7) 



The term "periodic" for the (p integral in Eq. ( IA.5b means that all fields (p over which 
we integrate have to be periodic in the imaginary time direction, ^(0,x) = 9(j3,x). 
This is essentially a consequence of the trace operation in the first line of Eq. ( IA.5I ): 
the partition function is formally reminiscent of a sum over transition amplitudes 
which have the same initial and final states at "times" and /3 . 
Let us, for convenience, introduce the two real fields (pi, (p2, 

(p = ^{(pi + i(p2) . (A.8) 



Then, the Lagrangian becomes 
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2n2 



The conjugate momenta are 



d{do(pi] 



i= 1,2. 



(A.9) 



(A. 10) 



Consequently, with f = <P2'^\ — (Pi'^2, which follows from Eqs. (IA.4I) . (lA.SI l. and 
jA.lOl) . we have 



-IJ.{(p2ni-(pi7t2). 



(A. 11) 



The integration over the conjugate momenta Tti , 7t2 can be separated from the inte- 
gration over the fields (pi , (p2 after introducing the shifted momenta 



ni=7ti-do(pi- IJ.(P2, n2 = 7t2-do(p2 + IJ.(pi. 



This yields 



where the new Lagrangian ^ now includes the chemical potential, 
1 



(A. 12) 



(A.13) 



^ 2 



(A. 14) 



Thus we see that the chemical potential produces, besides the expected term ji f, 
the additional term + 9t )/2. This is due to the momentum-dependence of f. 

In terms of the complex field (jO, the Lagrangian reads 



if =|(ao-/M)<Pl'-|V(pp-m2|(p|2-A|(pr 



(A.15) 



which shows that the chemical potential looks like the temporal component of a 
gauge field. We can now insert Eq. ( IA.13I ) into the partition function ( lA.SI l. The 
integration over conjugate momenta and over fields factorize, and the momentum 
integral yields an irrelevant constant A^, such that we can write 



Z = N I 9(pi9(p2exp f ^ . 

-/periodic Jx 



(A. 16) 
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In order to take into account Bose-Einstein condensation, we divide the field into 
a constant background field and fluctuations around this background, (p, 0,- + 9,. 
A nonzero condensate 0i + i(j)2 picks a direction in the U{1) degeneracy space and 
thus breaks the symmetry spontaneously. We can choose 02 = and thus may denote 
([) = Then, the Lagrangian (IA.14b becomes 

with the tree-level potential 

^/(0^) = ^^^0' + ^0^ (A.18) 

and terms of second, third, and fourth order in the fluctuations, 

■^^^^ = [-d^(pid^(pi-d^(p2d^((>2-2n{(p2do(pi-(pido(p2) 

+ {m^-H^){(pj + (p^) + 3X(p^(pf + X(p^(p^] , (A.19a) 
= -A(?)(pi((pf + (p|), (A. 19b) 

^^'' = -|(<pf + (A.19C) 

We have omitted the Unear terms since they do not contribute to the functional inte- 
gral. Note that the cubic interactions are induced by the condensate. 

In this appendix we are only interested in the tree-level contributions U{(j)'^) and 
in order to explain the basic calculation of the partition function for the sim- 
plest case. We therefore shall ignore the cubic and quartic contributions and 
. We introduce the Fourier transforms of the fluctuation fields via 

(p{X) = -^^Y^e-^'-^ipiK) = ^^.'('»"^+k.x)^(^)^ (^20) 
V J V ^ V IV 

with the four-momentum 

K={ko,k) = {-i(o„,k), (A.21) 

and with the Minkowski scalar product K X ^ koXQ — k • x = — (tco,, + k • x). (Al- 
though for convenience we have defined the time components with a factor / and thus 
can use Minkowski notation, the scalar product is essentially Euclidean.) The nor- 
malization is chosen such that the Fourier-transformed fields (p{K) are dimension- 
less. The 0-component of the four-momentum is given by the Matsubara frequency 
(On- To fulfill the periodicity requirement 9(0,x) = (p(j3,x) we need e"""'' = 1, i.e., 
a)„j3 has to be an integer multiple of 2k, or 

(On = 2nnT , neZ. (A.22) 



With the Fourier transform ( IA.20I ). and 
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s'^^'^YSKfi, (A.23) 
we have 

with the free inverse propagator in momentum space 

^0 ~ \ 2i^ko -K^ + m2 + A</)2 - ^2 j ■ (A.25) 

With Eqs. (IA.16b . (IA.24b and using that (p{K) = (p*{-K) (because q){X) is real) we 
can write the tree-level thermodynamic potential as 

Q. T 
— = InZ 

y V 



T 

;7(^^)-— In / ^(pi^(j[>2exp 



= U{f^) + — \ndet^^, (A.26) 

where the determinant is taken over 2x2 space and momentum space. Here we 
have used the general formula 

^D^^-^x-ix ^ (2;i:)^/2(detA)-i/2 , (A.27) 

for a Hermitian, positive definite matrix A, which is a generalization of the one- 
dimensional Gaussian integral 

r dxe---"'^ = . (A.28) 
7-oc V a 

To further evaluate the thermodynamic potential, we first define the tree-level 
masses 

m\=m^ + 3X^^, (A.29a) 

ml = m^+X^^ . (A.29b) 

Then, we obtain 
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lndet^<L^ = in]-[ -L [[-K^ + m\- ii^){-K^ + ml m') - 4^'^^] 



. j4 



L 

K 



J2 



J2 



(A.30) 



where we defined the quasiparticle energies 



with 



2 2 

m\ — 



(A.31) 



(A.32) 



Even at tree-level, the quasiparticle energies (lA.Bll i look complicated, but become 
simple in the noninteracting limit. 



A =0: 



and for vanishing chemical potential. 



At = 0: 



'2/1 ■ 



(A.33) 



(A.34) 



Further properties of these quasiparticle energies are discussed in the context of 
kaon condensation in CFL, see Sec. 14. 2.11 Next, we perform the sum over Matsubara 
frequencies in Eq. ( IA.30b . We use the result 



J2 



|+21n(l-e-^^/^ 



-const, 



(A.35) 



for a real number e^^, and where "const" is a temperature-independent constant. Be- 
fore we prove this result via contour integration in the complex plane, we use it to 
compute the final result for the tree-level thermodynamic potential. We insert Eq. 
( IA.35b into Eq. ( IA.30l l, the result into Eq. ( IA.26I ), and take the thermodynamic limit 
to obtain 



2T 



In 1 - e-^k 



4/T 



In 1 - e-^k 



£,7/7- 



(A.36) 

From this expression we can for instance compute the pressure P — — i2 /V. The first 
term in the integrand yields an infinite contribution which however is temperature- 
independent. We may thus use a renormaUzation such that the vacuum pressure van- 
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ishes. Then, for sufficiently large temperatures, where in particular (j) =0, particles 
and antiparticles yield the same contribution and we obtain 



A.1.1 Summation over bosonic Matsubara frequencies 

Here we prove Eq. (IA.35b via contour integration in the complex frequency plane. 
Especially for more complicated Matsubara sums this is a very useful technique as 
can be seen by applying the following method to the Matsubara sums in problems 

|A3]and|A2] 

First, in order to get rid of the logarithm, we write 

We now perform the sum in the integrand which, denoting = Tx, we write as a 
contour integral. 

The second identity follows from the residue theorem, 

^//^/W=i:Res/Wlz^z„' (A.40) 

where Zn are the poles of f{z) in the area enclosed by the contour C. If we can write 
the function / as f{z) ^ (p{z)/ V^(z), with analytic functions (p{z), "/(z), the residues 
are 

Res/(z)|„ =4H- (A-41) 

The contour C in Eq. (IA.39I I is chosen such that it encloses all poles of coth[a)/ {2T)] 
and none of 1/(0)^ - e^'). The poles of coth[a)/(2r)] are given by e'^/^r _ ^-m/iT ^ 
0, i.e., they are on the imaginary axis, (O = icOn with the Matsubara frequencies (On- 
In the above notation with the functions (p and y/, 

1 „(a/{2T) ,„-(a/{2T) 

cpico) = If — — , v/(«)=.'"/(^^)-.-'"/(^^^ 

(p(i(X)„) 1 

^ ^T7^ = 2 ' (A.42) 

\l/'{i(0„) co^ + £k 
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from which Eq. ( |A.39l l follows immediately. Next, we may deform the contour C 
(which consists of infinitely many circles surrounding the poles) and obtain 

V- 1 1 , 11,0) 

T} =■ = / dco—r T--coth — 

^o}},+el 2niJ-i^+-q co^-e^l 2T 

1 /■-'■"-'? 1 1 CO 
/ dco—^ 2oCOth — 

= / dco—:, ^coth — , (A.43) 

where we have changed the integration variable (O ^ — coin the second integral. We 
now use the residue theorem a second time: we can close the contour in the positive 
half-plane at infinity and pick up the pole at (O ~ Ek, 

(note the minus sign from clockwise contour integration). Here, 

(A.45) 

is the Bose distribution function. We have thus found 

1 ^ 1 _ 1 /'I 1 \ 
r^(2«7r)2+x2~7^(^2 + ^j • ^^-"^^^ 

Now we insert the result into the original expression ( |A.38l l and integrate over to 
obtain (with const denoting T-independent constants) 



^ Ji ^ X \2^ e^-l 

which is the result we wanted to prove 



const 



Y+2\n{\-e-^''i' ) + const , (A.47) 



A.2 Fermionic field 



To describe a system of non-interacting fermions with mass m we start with the 
Lagrangian 
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where \j/ ~ V^^}^- As for the bosons we are interested in adding a chemical poten- 
tial to this Lagrangian. To this end, we determine the conserved current as above, 
i.e., we first identify the global symmetry of the Lagrangian which is given by the 
transformation \j/ e^'" y/. The conserved current is 

which yields the conserved charge (density) 

/ = rV- (A.50) 



The conjugate momentum is 



-i\lf^ . (A.51) 



This means that in the case of fermions we need to treat i/a and i/^ as independent 
variables. The partition function for fermions is 



Tre 



^V^'^^V'sxp 

antiperiodic 



X 



(A.52) 



This has to be compared to the analogous expression for bosons, Eq. dA.SI l. Recall 
that the periodicity of the bosonic fields is a consequence of taking the trace in the 
operator formalism. In other words, the partition function in the path integral formal- 
ism can be derived from a transition amplitude with identical initial and final states. 
In the case of fermions, the fields in the path integral are Grassmann variables, as 
a consequence of the anticommutation relations of creation and annihilation oper- 
ators. In this case, the trace involves a transition amplitude where initial and final 
states differ by a sign. Therefore, in the fermionic partition function the integration 
is over antiperiodic fields V''(0,x) — ~\if{l5,x) and (0,x) = -v/^(j3,x). 
With the Hamiltonian 



^ = TT^or-^o =¥('T- V + m)v/, (A.53) 

(here and in the following we mean by the scalar product y ■ V the product where 
the Dirac matrices appear with a lower index yd we thus obtain 

(A.54) 

In this case we cannot separate the n ^ Xj/'' integration from the iff integration. Re- 
member that, in the bosonic case, this led to a new Lagrangian which contained 
the chemical potential not just in the term y"/!. Here, the Lagrangian with chemical 
potential simply is 



^antiperiodic 
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^ = \jf{ifd^ + fn-tn)\i/. (A.55) 

Note that again the chemical potential enters just like the temporal component of 
a gauge field that couples to the fermions. Analogously to the bosonic case, we 
introduce the (dimensionless) Fourier-transformed fields 

Y{X) = -Ly.^-''^''v{K) , W{X) = -^Y.^'^^'WiK) , (A.56) 
yV K yv K 

(note the different dimensionality of fields compared to bosons; here the field \lf{X) 
in position space has mass dimension 3/2). Again we denote ko = —ico„ such that 
K X = — (co„T + k-x). Now antiperiodicity, i/(0,x) — V''(j3,x), implies e"^"^ ~ 
— 1 and thus the fermionic Matsubara frequencies are 

co„ = {2n + 1)7ZT , neZ. (A.57) 

With the Fourier decomposition we find 

JX K ^ 

where the free inverse fermion propagator in momentum space is 

GqHk) = -fK^ + m . (A.59) 

Although not needed for the rest of the calculation in this appendix, let us intro- 
duce a useful form of the inverse propagator in terms of energy projectors. This 
form is convenient for more involved calculations such as done in chapters |4] and |5] 
Equivalently to Eq. ( |A.59t we can write 

Go\K) = ~Y,iko + ld- eEk) fA^ , ( A.60) 

e=± 



where E/^ — \/k^ + nfi, and where the projectors onto positive and negative energy 
states are given by 

Ai^Ul+ef'-^]. (A.61) 



2 V 

These (Hermitian) projectors are complete and orthogonal, 

K + = 1 ' ^k^k = ^.e'K ■ (A.62) 

The first property is trivial to see, the second follows with {7", / } = which follows 
from the general anticommutation property {y^ = 2g^^, and with (y-k)^ = 

From the form of the inverse propagator (IA.60b we can immediately read off the 
propagator itself. 



Go{K)^ . (A.64) 
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With the properties (IA.62I ) one easily checks that Gq 'Gq = 1. One can also rewrite 
(ESllas 

Let us now come back to the calculation of the partition function. For the functional 
integration we use 

Y\d^ld^k&^V ( -L'7/A7r?7 j = detD. (A.65) 

Note the difference of this integration over Grassmann variables 77^, 77 to the corre- 
sponding formula for bosons ( IA.27b . We obtain for the partition function 

Z.det^^detir-(^« + 'i)+'", --.-^ V (A.66) 

where the determinant is taken over Dirac space and momentum space, and where 
(7i , (72, (73 are the Pauli matrices. We can use the general formula 

det (^c D ) " det(AD - BD^'CD) , (A.67) 
for matrices A, B, C, D with D invertible, to get 

where we have used {a ■ k)^ = k^. With ko = —icOn we can write this as 
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where, in the second term of the second line, we have replaced (0„ by — ©„ which 
does not change the result since we sum over all « G Z. The third line can be easily 
checked by multiplying out all terms. 

Next we need to perform the sum over fermionic Matsubara frequencies. This is 
similar to the bosonic case and yields 
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Yin " = — +21n(l+e~'^^/^) + const. (A.70) 

Using this result to evaluate Eq. ( |A.69l l and taking the thermodynamic limit yields 
the thermodynamic potential Q. = TlnZ, 



— = -2 



V J {27t) 



E, + T\n( 1 + rin ( 1 



(A.71) 



The overall factor 2 accounts for the two spin states of the spin- 1/2 fermion. To- 
gether with the particle/antiparticle degrees of freedom we recover all four degrees 
of freedom of the Dirac spinor. Again we conclude this section by computing the 
pressure for large temperatures, 

Comparing with the bosonic pressure ( |A.37| i we see that for large T a single 
fermionic degree of freedom contributes 7/8 times as much to the thermal pressure 
as a single bosonic degree of freedom. 



A.2.1 Summation over fermionic Matsubara frequencies 

It remains to prove Eq. ( IA.70b by summing over fermionic Matsubara frequencies. 
As for the bosonic case, we write 



9 , 



j2 
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-(2„ + l)%^+^ + ?^"[^ + (^" + ^)'^']- 



(A.73) 

This time, we need to use the tanh instead of the coth when we write the sum in 
terms of a contour integral. 
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(A.74) 



(We have denoted e^. = xT .) The poles of tanh[ft)/(2r)] are given by the zeros of 
g(H/(2r) _|_ g-(u/(2r)^ j jjjgy located at ; times the fermionic Matsubara fre- 
quencies, (0 = icOn- The contour C encloses these poles and none of the poles of 
1 / (0)2 — e^). Then, with the residue theorem and with 



,a>/i2T) _ -a>/i2T) 



J_ (a>l(2T)^^-a>l(2T)\ 

do}\ J 
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(A.75a) 
(A.75b) 
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one confirms Eq. ( |A.74| i. We can now close the contour in the positive hatf-plane to 
obtain 



1 /■'■-+'' , 1 a 
aco — =r tanh - 



where 



/.(e)-^^ (A.77) 



is the Fermi distribution function. Inserting this resuh into Eq. ( IA.73I ) yields 



const 



l\v,[\+e-^^l' ) + const, (A.78) 



which proves Eq. ( |A.70| i. 



Problems 

A.l. Matsubara sum for boson loop 

Show via contour integration that 



2^ ^2 _ c2U/',-,„ _ t„A2 _ c2l ^ 2- 



(A.79) 

with kQ = —ico„, po = —icOin bosonic Matsubara frequencies, and ei , £2 > 0. 
A.2. Matsubara sum for fermion loop 

Prove via contour integration the following result for the summation over fermionic 
Matsubara frequencies, 

t(fcre?)Po + <?o)2-e|] 

1 y (£1 -ei^i)(£2-e2fc)/f (-ei£i)/F(e2£2) 



4£i£2,,,^=± ^o-ei£i+e2e2 /B(-ei£i +e2£2) 

where A:o = — /O),, with fermionic Matsubara frequencies (»„, and qo = —i(0,„ with 
bosonic Matsubara frequencies (0,„, and where ^i,^2i£ij£2 > are real numbers. 
The result of this problem is used in the calculation of the neutrino emissivity in 
chapter|5] 
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2SC phase Color superconductor in which strange quarks and quarks of one color 
remain unpaired. Because of the asymmetry induced by the strange quark mass, 
viable candidate for the ground state of quark matter at moderate chemical potential. 
In these lectures we discuss the 2SC phase in the context of neutrino emissivity, to 
illustrate the effect of both paired and unpaired quarks. 

AdS/CFT correspondence Theoretical tool not discussed in these lectures, but an 
interesting approach to tackle QCD at strong coupling. The idea is that - relatively 
simple - calculations in the gravity approximation of a certain string theory provide 
results for the - otherwise hard to access - strong coupling limit of a corresponding 
("dual") field theory. The problem is that currently no gravity dual of QCD is known. 

anomalous propagator Technically speaking, off-diagonal components of the 
propagator in Nambu-Gorkov space; nonzero in the case of a superconductor or 
a superfluid. More physically speaking, anomalous propagators describe a fermion 
which is, via the Cooper pair condensate, converted into a fermion hole. 

asymptotic freedom Important property of QCD which says that the running cou- 
pling constant of QCD becomes small for large exchanged momenta. For our con- 
text this means that quarks at large densities, where the distance between them is 
small and hence the exchanged momentum large, are weakly interacting; quarks at 
infinite density are free. In compact stars, however, the density is large, but by no 
means asymptotically large. 

axial anomaly Non-conservation of the axial current in QCD. In our context of 
(moderately) dense matter originating mainly from instantons which are certain 
semi-classical gauge field configurations. Leads to an explicit breaking of the ax- 
ial U{1)a, which is a subgroup of the chiral group, and thus gives a large mass to 
the 77'. 

bag model (MIT bag model) Simple model to take into account confinement. Via 
the bag constant, an energy penalty is introduced by hand for the deconfined phase. 
The model amounts to the picture of a hadron as a bag which confines the quarks; 
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the bag exerts an external pressure on the quarks, given by the bag constant. In our 
astrophysical context, the bag model is a simple way to compare free energies of 
dense quark matter and dense nuclear matter 

BCS theory Original theory for electronic superconductors, developed in 1957 by 
Bardeen, Cooper, and Schrieffer. Many concepts and approximations can be adopted 
for nuclear and quark matter. In color-superconducting quark matter, an important 
difference to BCS theory is the parametric dependence of the pairing gap on the 
coupling constant due to long-range interactions via magnetic gluons. 

/3 -decay Process due to the weak interaction of the form n ^ p + e + Ve in nu- 
clear matter and c/ — > m + e + v^. in quark matter. Relevant in these lectures for two 
reasons: firstly, equilibrium with respect to this process (j3 -equilibrium) yields im- 
portant constraints for the chemical potentials and secondly, this process contributes 
to the neutrino emissivity which in turn is responsible for the cooling of a compact 
star. 

Bogoliubov coefficients Momentum-dependent coefficients in the theory of su- 
perconductivity and superfluidity which characterize the mixing of fermions and 
fermion holes due to Cooper pair condensation. In these lectures, the Bogoliubov 
coefficients arise naturally in the calculation of the neutrino emissivity in color- 
superconducting quark matter. 

Cabibbo-Kobayaslii-Masliawa (CKM) matrix Matrix that characterizes the rel- 
ative strength of the weak interaction for different quark flavors. In these lectures 
relevant for the calculation of the neutrino emissivity in quark matter 

chiral symmetry For massless quarks, QCD possesses a global symmetry for 
right- and left-handed quarks separately, called chiral symmetry. This symmetry can 
be spontaneously broken, giving rise to Goldstone modes. These Goldstone modes 
(or pseudo-Goldstone modes in the case of nonzero quark masses) are for instance 
pions and kaons. In these lectures we discuss kaon condensation in nuclear and 
quark matter (in quark matter, chiral symmetry is spontaneously broken in the CFL 
phase). 

color superconductivity Cooper pair formation and condensation in cold and 
dense quark matter, analogous to electronic superconductivity in metals. If quark 
matter is present in compact stars, it can be expected to be a color superconductor 

color-flavor locking (CFL) Ground state of three-flavor quark matter at asympto- 
tically large densities. Particularly symmetric color superconductor where the order 
parameter is invariant only under simultaneous color and flavor transformations. 
May persist down to densities where the hadronic phase takes over or may be re- 
placed before this transition by a different color superconductor because of the ef- 
fects of the strange quark mass. 

compact star Very dense astrophysical object with a mass close to the sun's mass 
and a radius of about ten kilometers. The term shows our ignorance of the exact 
composition of these objects. They may be neutron stars, hybrid stars, or quark stars. 
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In a more general terminology, compact star also is used to include white dwarfs and 
black holes, neither of which are the subject of these lectures. 

constituent quark mass Quark mass including the quark's interaction energy in 
a baryon, such that the sum of the three constituent quark masses adds up to the 
baryon mass. More generally, in dense matter the density-dependent "constituent" 
quark mass includes any finite-density effects. Can be hundreds of MeV larger than 
the current quark mass. 

Cooper pairs Microscopic explanation for superfluidity and superconductivity 
within BCS theory. Arise from an instability of the Fermi surface in the presence 
of an arbitrarily small interaction. In compact stars, there are possibly Cooper pairs 
of neutrons, protons, hyperons and/or quarks. 

crust Outer, km thick, layer of a neutron star or hybrid star Composed of ordinary 
nuclei which form a crystalline structure and which, upon increasing the density and 
thus going further inside the star, become more and more neutron rich. In the inner 
crust a neutron superfluid is immersed in the lattice of nuclei. Quark stars have, if at 
all, much thinner crusts. 

current quark mass Quark mass without effects from the interactions with other 
quarks and gluons, see also constituent quark mass. Since interactions become weak 
at asymptotically large densities (much larger than densities in compact stars), cur- 
rent and constituent quark masses become identical in this limit. 

dense matter In these lectures, dense matter means matter at densities of a few 
times nuclear ground state density, as expected in the interior of compact stars. Gov- 
erned by the strong interaction, and thus very difficult to describe theoretically. We 
discuss several theoretical concepts and sometimes have to escape to lower or even 
higher densities, just to make life simpler. 

equation of state Relation between the pressure and the energy density for a given 
form of dense matter. In our context, the equation of state determines, together with 
the TOV equation, the mass-radius relation of a compact star In particular, a stiff 
(soft) equation of state allows for a large (small) maximum mass. 

Goldstone boson Massless boson arising from spontaneous symmetry breaking of 
a global symmetry. The only exact (i.e., truly massless) Goldstone boson in dense 
matter is the one associated to superfluidity, i.e., to the breaking of baryon number 
conservation. Such a mode exists in a nuclear superfluid as well as in the color-flavor 
locked phase. 

hybrid star Compact star with a quark matter core and a nuclear mantle. Most 
likely scenario to find quark matter in a compact star. 

hyperon Baryon with nonzero strangeness. Hyperons may occur in hadronic mat- 
ter at sufficiently large densities. In these lectures only discussed briefly, in the con- 
text of Walecka-like models. 
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incompressibility Thermodynamic property of nuclear matter at the saturation 
density, sometimes also called compression modulus. Can be (at least indirectly 
and approximately) determined in the experiment and thus yields a value that can, 
among other quantities, be used to fit the parameters of theoretical models, in these 
lectures the coupling constants of the Walecka model with scalar interactions. 

kaon condensation Possible example of Bose-Einstein condensation in a compact 
star May appear at sufficiently large densities. Is possible not only in nuclear matter, 
but also in quark matter, where kaons exist in the CFL phase. These kaons carry the 
same quantum numbers as the usual kaons, however are made of two quarks and 
two quark holes. 

Kepler frequency Absolute upper limit for the rotation frequency of compact stars 
beyond which mass shedding at the equator sets in. Given by the equality of the 
centrifugal and gravitational forces (more precisely, the general relativistic version 
thereof). For typical compact stars in the ms^' regime, i.e., for some pulsars actually 
observed rotation frequencies are not too far from that limit. Below that limit stars 
can suffer from other rotational instabilities, for instance the r-mode instability. 

Landau mass Effective mass of (nonrelativistic) fermions at the Fermi surface, in 
the framework of Landau's Fermi liquid theory. In these lectures, the Landau mass 
for nucleons is mentioned in the context of the Walecka model where its experimen- 
tal value serves to fit the parameters of the model. 

lattice QCD QCD on the computer Powerful method to perform calculations from 
first principles. Not discussed in these lectures, mostly because lattice QCD is cur- 
rently unable to provide results at large chemical potential and small temperature 
because of the so-called sign problem. 

Low-mass X-ray binary (LMXB) System of two stars, where a pulsar is accreting 
matter from its companion which has a mass typically smaller than one solar mass 
(as opposed to high-mass X-ray binaries where the companion has a mass larger 
than about ten solar masses). Measured rotation frequencies of pulsars in LMXBs 
are mentioned in our brief discussion of the r-mode instability of rotating compact 
stars. 

magnetar Compact star with unusually large magnetic field, up to lO'^G at the 
surface and possibly larger in the interior 

Matsubara frequency In thermal field theory, the time direction in Minkowski 
space becomes imaginary and compact, giving rise to Euclidean space with discrete 
energies, given by the Matsubara frequencies. In these lectures we mostly consider 
the zero-temperature limit, but in some instances we have to perform a sum over 
Matsubara frequencies. 

mixed phase Coexistence of two (or more) phases which occupy certain volume 
fractions - for instance bubbles of one phase immersed in the other phase - in a 
given total volume. In our context, global charge neutrality, as opposed to local 
charge neutrality, allows for mixed phases for instance of nuclei and nuclear matter 
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or quark and hadronic matter. These phases may be disfavored by large surface 
energy costs. 

Nambu-Gorkov space Contains Nambu-Gorkov spinors which arise from a dou- 
bhng of the fermionic degrees of freedom in the theoretical description of super- 
conductors and superfluids. Allows to introduce Cooper pairing in the off-diagonal 
elements of the Nambu-Gorkov propagators. See also anomalous propagators. 

Nambu-Jona-Lasinio (NJL) model Phenomenological model, not discussed in 
these lectures, where the QCD interaction between quarks is replaced by a point- 
like four-quark interaction. Since it has attraction in the same channels as QCD, this 
model is frequently used to describe color-superconducting quark matter at moder- 
ate densities. 

neutron star Compact star made of neutron-rich nuclear matter. In some literature 
the term neutron star is also used to include the possibility of a quark matter core. 
Mostly, also in these lectures, these stars are called hybrid stars. 

nuclear pasta Mixed phase of ordinary nuclei (ions) and nuclear matter, typically 
found in the inner cores of neutron stars. Because different geometries can be real- 
ized - spheres, rods, slabs, the latter two reminiscent of spaghetti or lasagna - these 
phases have been termed nuclear pasta. In these lectures we discuss the possibility 
of mixed phases of quark and hadronic matter. 

pion condensation Bose-Einstein condensation of pions in nuclear matter Al- 
though pions are lighter than kaons in the vacuum, kaon condensation seems to be 
more likely in dense nuclear matter. Therefore, in these lectures, kaon condensation, 
not pion condensation, is discussed. 

pseudo-Goldstone boson Less impressive brother of the Goldstone boson, arising 
from spontaneous breaking of a global symmetry which is broken explicitly by a 
small amount (small compared to the scale of the spontaneous breaking). Light, but 
not exactly massless. Dense matter is full of pseudo-Goldstone modes, for instance 
mesons in nuclear matter or color-flavor-locked quark matter, arising from the spon- 
taneous breaking of chiral symmetry which is explicitly broken by quark masses. 

pulsar Star whose radiation is observed in periodic pulses. Pulsars are rotating 
compact stars with large magnetic fields; their apparent pulsation is due to the align- 
ment of the radiation in a beam along the magnetic axis. When the magnetic axis is 
different from the rotation axis, the beam may point towards the earth periodically, 
just as the light of a lighthouse flashes periodically when you observe it from the 
beach. 

pulsar glitch Sudden spin-up of a rotating compact star. Not discussed in detail in 
these lectures but very interesting phenomenon since closely related to the micro- 
scopic physics, presumably to crystalline structures and vortices in the star. 

QCD phase diagram Collection of equilibrium states of QCD, typically depicted 
in the plane of quark (or baryon) chemical potential and temperature. We roughly 
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know where compact stars sit in this diagram, but we do not know the phase(s) that 
occupy this region of the diagram. These lectures are about exploring this unknown 
territory. 

quantum chromodynamics (QCD) Theory of the strong interaction. Governs the 
physics that determines the ground state of dense matter present in a compact star. 
In these lectures we perform one explicit calculation in QCD and discuss several ef- 
fective approaches to this very elegant, but for most practical purposes very difficult, 
theory. 

quarkyonic matter Form of dense matter covering a large portion of the QCD 
phase diagram for the case of asymptotically large number of colors. Not discussed 
in these lectures because there are only three colors in the real world. However, it is 
a viable option that a small region of quarkyonic matter survives and thus becomes 
also important for compact stars. 

quasiparticle Term originally used in condensed matter physics and carried over 
to dense QCD matter Absorbs interactions of the original particles into effective 
new particles. For instance, quasiparticles in a superconductor are gapped due to the 
attractive interaction between the original particles. 

r-modes Non-radial oscillations of a star with the Coriolis force as the restoring 
force. Interesting for dense matter physics because they grow unstable in a pulsar 
unless the matter inside the star is sufficiently viscous. 

rotated electromagnetism Effect in some color superconductors which is respon- 
sible for them being no electromagnetic superconductors. Therefore important for 
the physics of compact stars since magnetic fields penetrate these color supercon- 
ductors. Technically speaking, rotated electromagnetism refers to a gauge boson 
which is a mixture of a gluon and the photon. 

saturation density Density at which the binding energy is minimized, here always 
used in the context of nuclear matter for which the saturation density is approxi- 
mately 0.15 baryons per fm^ and the corresponding binding energy per nucleon is 
about 16MeV. 

sign problem Problem of QCD lattice calculations at finite values of the baryon 
chemical potential. For finite chemical potential, the action, more precisely the 
quark determinant in the functional integral of the partition function, loses its pos- 
itivity and even becomes complex. This makes the probabilistic sampling method 
("Monte Carlo method"), on which lattice QCD is based, unfeasible. In our context 
this means that currently there is no input from lattice calculations to the properties 
of dense matter. 

strange quark matter hypothesis Hypothesis that strange quark matter, not nu- 
clear matter, is the ground state at zero pressure. The hypothesis does not contradict 
our existence since, even if the hypothesis is true, the transition from nuclear matter, 
made of u and d quarks, to strange quark matter is essentially forbidden. We discuss 
that, within the bag model, the strange quark matter hypothesis is true if the bag 
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constant is between a lower bound (since we know that ordinary nuclear matter is 
stable with respect to two-flavor quark matter) and an upper limit (beyond which 
nuclear matter is absolutely stable). 

strange star (quark star) Compact star made entirely out of quark matter, thus the 
most radical scenario for quark matter in compact stars. 

strangelet Small nugget of strange quark matter. Stretching the original meaning a 
bit - well, from femtometers to kilometers - a strange star is a huge strangelet. Rel- 
evant for us in the context of the strange quark matter hypothesis; since strangelets 
would convert neutron stars into strange stars, the unambiguous observation of a 
single neutron star would invalidate the strange quark matter hypothesis, provided 
that there are enough strangelets in the cosmos to hit neutron stars. 

supernova Compact stars are expected to be born in (type II) supernova explo- 
sions, where a giant star, after burning its nuclear fuel, undergoes a gravitational 
collapse. The energy of the explosion is mostly released in the form of neutrinos. 
The theoretical description of supernovae requires very complicated hydrodynami- 
cal simulations. 

Tolman-Oppenheimer-Volkov (TOV) equation Differential equation from gen- 
eral relativity for the mass, pressure, and energy density as functions of the distance 
from the center of the star. In connection with the equation of state, which relates 
energy density and pressure, used to compute the mass-radius relation for a compact 
star. 

unpaired quark matter Term used for (dense) quark matter which does not form 
Cooper pairs and thus is no color superconductor Since dense quark matter is ex- 
pected to be some kind of color superconductor, completely unpaired dense quark 
matter is unlikely to exist. Therefore mostly used for reference calculations or when, 
for the computed quantity, it is a good approximation to paired quark matter. 

Urea process Most efficient process for neutrino emission, and thus for the cooling 
of the star In quark matter the direct Urea process is given by m + e c/ + Vg and 
variants thereof. We compute the emission rate of this process in detail in these 
lectures. In the modified Urea process, a spectator particle is added which increases 
the available phase space. 

viscosity (bulk/shear) Transport coefficients of nuclear or quark matter relevant in 
particular in the context of rotation and oscillation of the star. Requires microscopic 
calculation of processes typically governed by the weak interaction. Not discussed 
in detail in these lectures. See also r-modes. 

Walecka model Phenomenological model for interacting nuclear matter, based on 
Yukawa couplings of the nucleons with the a and o meson. Used for extrapolation 
to large densities after fitting the parameters of the model at saturation density. Dis- 
cussed in these lectures as a basic example for numerous more complicated nuclear 
models of similar kind. 
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white dwarf Dense star with a mass of about the sun's mass and radius of a few 
thousand kilometers, which makes it less dense than a neutron star. Composed of 
nuclei immersed in a degenerate electron gas. 



